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Foreword 


For  the  last  couple  of  decades,  attempts  have 
been  made  to  develop  some  general  under¬ 
standing,  and  ultimately  a  theory,  of  systems 
that  consist  of  many  interacting  components 
and  many  hierarchical  layers.  It  is  common  to 
call  these  systems  complex  because  it  is  impossi¬ 
ble  to  reduce  the  overall  behaviour  of  the 
system  to  a  set  of  properties  characterising  the 
individual  components.  Interaction  is  able  to 
produce  properties  at  the  collective  level  that 
are  simply  not  present  when  the  components 
are  considered  individually.  As  an  example,  one 
may  think  of  mutuality  and  collaboration  in 
ecology.  T  he  function  of  any  ecosystem  depends 
crucially  on  mutual  benefits  between  the  differ¬ 
ent  species  present.  0  ne  example  is  the  relation 
between  legumes,  such  as  peas  and  beans,  and 
their  associated  nitrogen-fixing  bacteria:  the 
bacteria  collects  nitrogen  for  the  legume,  which 
in  turn  produces  carbohydrates  and  other 
organic  material  for  the  bacteria.  Clearly  this 


XI 


crucial  arrangement  cannot  be  studied  by  focusing  on,  say, 
the  legume  and  neglecting  the  bacteria;  the  ecological  func¬ 
tion  emerges  first  when  the  different  components  are  brought 
together  and  interaction  is  taken  into  account. 

Another  important  feature  of  complex  systems  is  their  sensi¬ 
tivity  to  even  small  perturbations.  The  same  action  is  found 
to  lead  to  a  very  broad  range  of  responses,  making  it  exceed¬ 
ingly  difficult  to  perform  prediction  or  to  develop  any  type  of 
experience  of  a  "typical  scenario.”  T  his  must  necessarily  lead 
to  great  caution:  do  not  expect  what  worked  last  time  to  work 
this  time.  The  situation  is  exacerbated  since  real  systems 
(ecological  or  social)  undergo  adaptation.  This  implies  that 
the  response  to  a  given  strategy  most  likely  makes  the  strat¬ 
egy  redundant.  An  example  is  the  effect  of  using  the  same 
type  of  antibiotic  against  a  given  type  of  bacteria.  Evolution 
soon  ensures  that  the  bacteria  develop  resistance  and  make 
the  specific  type  of  antibiotic  useless.  T  hat  complex  systems 
adapt  and  change  their  properties  fundamentally  as  a  result 
of  the  intrinsic  dynamics  of  the  system  is  clearly  extremely 
important.  N  evertheless,  for  the  sake  of  simplicity  adaptation 
is  often  neglected  in  model  studies.  Sometimes  assuming  the 
existence  of  a  stationary  state  might  be  justified  (e.g.,  if  one  is 
interested  in  "toy”  models  of  the  flow  of  granular  material 
under  a  controlled  steady  input  of  grains).  But  if  one  is  deal¬ 
ing  with  more  complex  situations  such  as  in  ecology,  and 
even  more  when  considering  social  and  political  systems, 
ignoring  adaptation  is  very  likely  to  lead  to  erroneous 
conclusions. 

We  know  from  studies  of  Self-Organised  Critical  models, 
which  the  present  book  alludes  to  (for  more  see  P.  Bak,  H  ow 
Ndture  Works,  Oxford  U  niversity  Press,  1997  and  H  .].  Jensen, 
Sdf-Orgenized  Criticdlity,  Cambridge  University  Press,  1998), 
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that  the  correlations  and  general  behaviour  exhibited  by  these 
model  systems  are  entirely  determined  by  the  assumed 
boundary  conditions  or  the  applied  drive.  The  lesson  to  be 
learned  from  this  is  that  complex  systems  cannot  be  studied 
independently  of  their  surroundings.  Understanding  the 
behaviour  of  a  complex  system  necessitates  a  simultaneous 
understanding  of  the  environment  of  the  system.  In  model 
studies,  one  assumes  often  that  the  surroundings  can  be  repre¬ 
sented  by  one  or  the  other  type  of  "noise,”  but  this  is  just  a 
trick  that  allows  one  to  proceed  with  the  analysis  without 
understanding  the  full  system  under  consideration.  It  is  very 
important  to  appreciate  that  the  "drive"  or  the  "noise"  are 
equally  crucial  to  the  understanding,  as  is  the  analysis  of  the 
"system"  itself.  One  should  bear  in  mind  that  the  separation 
into  system,  drive,  noise,  surroundings,  etc.  is  rather  arbitrary 
and  is  far  from  representing  a  complete  analysis. 

From  these  considerations,  we  see  that  it  is  vitally  important 
to  consider  warfare  as  a  complex  system  that  is  linked  and 
interacts  (in  a  coevolving  way)  with  the  surrounding  socio- 
economical  and  political  context.  From  that  perspective,  the 
present  book  is  a  "work  in  progress"  and  a  preliminary  first 
step  along  the  road  in  helping  to  analyse  and  structure  these 
difficult  and  serious  issues.  Forgetting  that  war  and  warfare 
are  an  intimate  part  of  a  much  larger  complex  system  will 
lead  to  incomplete  and  even  dangerously  incorrect  conclu¬ 
sions.  Applying  the  approach  of  Complexity  Theory  to 
warfare  leads  one  to  the  self-consistent  realisation  that  war¬ 
fare  will  have  to  be  analysed  in  its  larger  context.  Further 
work  will  need  to  examine  how  coevolution  across  the  entire 
network  of  military,  socioeconomical,  and  political  interac¬ 
tions  leads  firstly  to  emergent  effects  at  higher  levels,  and  of 
equal  importance  how  such  effects  lead  to  coevolution  at  the 
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higher  level.  It  will  also  be  important  to  consider  the  robust¬ 
ness  of  such  networks,  and  their  vulnerability  to  damage. 


H  enrikjeldtoft  Jensen 
Professor  of  M  athematical  Physics 
Department  of  M  athematics 
Imperial  College,  London 
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Chapter  1 


Complexity  in 
Natural  and 
Economic  Systems^ 


In  this  chapter  we  consider  some  of  the  key 
ideas  of  Complexity  Theory  as  applied  to 
natural  systems.  Having  established  these  key 
ideas,  in  Chapter  2  we  begin  to  see  how  these 
ideas  map  across,  in  a  broad  conceptual  sense,  to 
the  dynamics  of  conflict  in  an  Information  Age 
environment.  In  later  chapters,  we  will  look  in 
more  detail  at  the  evidence  and  the  type  of  mod¬ 
elling  that  emerges  from  this  conceptual 
connection.  In  an  organisational  context,  we 
argue  that  complexity  provides  an  explanatory 
framework  of  interrelationships,  both  metaphor¬ 
ically  and  analogously,  of  how  individuals  and 
military  organisations  interact,  relate,  and  evolve 


he  contribution  of  D  r.  M  aurice  Passman  to  this 
chapter  is  gratefully  acknowledged. 
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within  a  larger  "ecosystem.”  Complexity  explains  why  inter¬ 
ventions  may  have  unanticipated  consequences,  but  also 
explains  how  combat  effects  follow  from  these  consequences. 
The  intricate  interrelationships  of  elements  within  a  complex 
system  give  rise  to  multiple  chains  of  dependencies.  Change 
happens  in  the  context  of  this  intricate  intertwining  at  all 
scales.  We  become  aware  of  change  only  when  a  different  pat¬ 
tern  becomes  discernible.  But  before  change  at  a  macro  level 
can  be  seen,  it  is  taking  place  at  many  micro  levels  simulta¬ 
neously.  Hence,  microcomponent  interaction  and  change 
leads  to  macrosystem  evolution. 

In  a  previous  book,^  we  considered  some  of  the  issues  to  be 
addressed  at  the  political/  military  level  as  a  consequence  of 
such  emergent  behaviour  and  the  "resultant  likelihood  of 
complex  and  unexpected  interactions,  arising  from  previ¬ 
ously  unexpected  sources."  T  he  diligent  reader  is  directed  to 
that  work  for  further  discussion  of  these  issues.  H  ere,  we  sim¬ 
ply  wish  to  add  weight  to  the  points  made  by  Professor  Jensen 
in  his  Foreword  to  this  present  work.  In  all  that  follows,  the 
recursion  of  the  process  up  to  this  political/  military  level 
must  be  kept  in  mind,  and  will  be  one  of  the  key  areas  of 
future  research. 

Our  lead  is  taken  from  current  military  doctrinal  thought 
both  in  the  UK  and  in  the  United  States  (particularly  the 
U.S.  Marine  Corps  [1]).  The  Chief  Analyst  of  the  UK 
Defence  Science  and  Technology  Laboratory  (DstI),  Roger 
Forder,  makes  the  following  point  in  his  discussion  of  the 
future  of  defence  analysis  [2]: 


hapter  1  of:  M  offat  J  (2002).  Command  and  Control  in  the  Information  A  ga 
T  he  Stationery  0  ffice.  L  ondon,  U  K . 
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0  ne  effect  of  the  human  element  in  conflict  situations  is  to  bring 
a  degree  of  complexity  into  the  situation  such  that  the  emergent 
behaviour  of  the  system  as  a  whole  is  extremely  difficult  to  pre¬ 
dict  from  the  characteristics  and  relationships  of  the  system 
elements.  Detailed  simulation,  using  agent- based  approaches,  is 
always  possible  but  the  highly  situation- specific  results  that  it 
provides  may  offer  little  general  understanding  for  carrying  for¬ 
ward  into  robust  conclusions  of  practical  significance  U  sable 
theories  ofcompl&dty,  which  would  allow  understanding  of 
emergent  behaviour  rather  than  merely  its  observation,  would 
therefore  have  a  great  deal  to  offer  to  some  of  the  central  problems 
facing  defence  analysis.  I  ndeed  they  might  well  be  thesinglemost 
desirable  theoretical  d&dopment  that  we  should  seek  over  the 
next  few  years. 

A  similar  thought  was  aired  at  a  Royal  U  nited  Services  Insti¬ 
tute  (RUSI)  conference  on  future  Intelligence,  Surveillance, 
Target  Acquisition,  and  Reconnaisance  (ISTAR)  [3],  Vice 
Admiral  Cebrowski,^  U.S.  Navy,  centred  his  keynote  address 
on  Network  Centric  Warfare  (NCW)  as  the  capstone  concept 
for  the  U.S.  N  avy  after  N  ext.  H  e  described  it  in  terms  of  the 
achievement  of  I  nformation  Superiority  with  characteristics  of 
gross  asymmetries  and  a  diversity  of  "players.”  T  hese  ideas  are 
explicitly  derived  in  his  description  from  the  new  physics  of 
nonlinearity,  complexity,  and  chaos  as  exemplified  by  the 
Santa  Fe  Institute  corpus  of  ideas.  The  relationship  between 
complexity  and  "information-based"  warfare  is  (as  he 
described  it)  less  deterministic  and  more  emergent;  less  focused 
on  the  physical,  and  more  behavioural;  less  focused  on  things, 
and  more  on  relationships.  Command  and  Control  (C2) 


ow  head  of  the  0  ffice  of  Force  T ransformation,  T  he  Pentagon,  U  .5. 
DoD. 
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emphasises  speed,  sharing,  and  decentraiisation.  in  summary, 
ADM  Cebrowski  defined  NCW  as  the  robust  networking  of 
weii-informed,  geographicaiiy  dispersed  forces. 

1/1/  here  do  these  ideas  come  from? 

I  n  looking  at  where  these  key  ideas  of  complexity  come  from, 
let  us  make  a  start  by  considering  systems  and  their  behaviour 
in  the  natural  world,  in  the  classical  view,  such  physical  or  bio¬ 
logical  processes  are  reducible  to  a  few  fundamental 
interactions.  T  his  leads  to  the  idea  that  under  well-defined  con¬ 
ditions,  a  system  governed  by  a  given  set  of  laws  will  follow  a 
unique  course  (likethe  planets  of  the  solar  system).  M  oreover,  a 
slight  change  in  the  cause  will  likewise  produce  a  slight  change 
in  the  effects  (i.e.,  the  system  is  linear  in  nature).  Recently,  an 
increasing  amount  of  experimental  data  challenging  this  idea 
has  become  available  and  this  imposes  a  new  attitude  concern¬ 
ing  the  description  of  nature,  in  natural  systems,  under 
appropriate  conditions,  a  multitude  of  self-organisation  phe¬ 
nomena  on  a  macroscopic  scale  (a  scale  of  an  order  of 
magnitude  larger  than  the  range  of  fundamental  interactions) 
in  the  form  of  spatial  or  temporal  patterns  may  be  generated. 

A  Simple  Example 

To  illustrate  this,  let  us  consider  a  simple  thermodynamic 
thought  experiment.  Imagine  a  layer  of  fluid  limited  by  two 
horizontal  parallel  plates  whose  lateral  dimensions  are  much 
longer  than  the  width  of  the  layer,  as  shown  in  Figure  1.1. 


Figure  1.1:  T  woH  orizontal  Plates  Containing  a  Layer  of  Fluid 
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Left  to  reach  equilibrium,  the  fluid  will  rapidly  tend  to  a 
homogeneous  state  that  is  statistically  identical.  T  he  homoge¬ 
neity  of  this  system  extends  to  all  of  its  properties,  particularly 
to  its  temperature,  which  will  be  the  same  at  all  parts  of  the 
fluid  and  equal  to  the  temperature  of  the  limiting  plates  or, 
alternatively,  to  the  temperature  of  the  "external”  world.  All  of 
these  properties  are  characteristic  of  a  system  in  a  particular 
state,  the  state  of  equilibrium,  for  which  there  is  neither  bulk 
motion  nor  temperature  difference  with  the  outside  world. 
What  occurs  if,  for  example,  the  temperature  on  a  small  sec¬ 
tion  on  one  of  the  plates  is  temporarily  perturbed?  At 
equilibrium,  this  temperature  perturbation  has  no  influence, 
since  the  temperature  rapidly  becomes  uniform  again  and 
equal  to  its  initial  value.  In  other  words,  the  perturbation  dies 
out;  the  system  keeps  no  track  of  it.  Such  a  state  is  said  to  be 
dsymptoticdlly  stdble. 

From  the  standpoint  of  a  very  small  observer  inside  the  sys¬ 
tem,  not  only  does  the  homogeneity  of  the  fluid  make  it 
impossible  for  the  development  of  an  intrinsic  concept  of 
space,  but  also  the  stability  of  the  state  of  equilibrium  eventu¬ 
ally  makes  all  time  instances  identical.  It  is  therefore 
impossible  for  this  observer  to  develop  an  intrinsic  concep¬ 
tion  of  correlation  or  coincidence  (e.g.,  things  happening  at 
the  same  time,  or  in  the  same  place).  We  can  increase  the 
complexity  of  the  system  by,  for  instance,  heating  the  fluid 
layer  from  below.  In  doing  this,  we  communicate  energy  to 
the  system  in  the  form  of  heat.  M  oreover,  as  the  temperature 
of  the  lower  plate  is  now  higher  than  the  upper,  the  equilib¬ 
rium  condition  is  violated.  In  other  words,  by  applying 
externdl  constreints  to  the  system,  we  do  not  permit  the  system 
to  reach  equilibrium.  T  he  presence  of  an  external  constraint 
therefore  implies  energy  flux  and  vice-versa. 
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Now  suppose  at  first  that  the  constraint  is  weak,  i.e.  the  change 
in  temperature,  Ar,  is  small.  The  system  will  again  adopt  a 
simple  and  unique  state  in  which  the  only  active  process  is  a 
transfer  of  heat  from  the  lower  to  the  upper  plate,  from  which 
heat  is  lost  to  the  external  world.  T  he  only  difference  from  the 
state  of  equilibrium  is  that  temperature,  and  consequently  den¬ 
sity  and  pressure,  are  no  longer  uniform.  T  hey  vary  from  warm 
regions  to  cold  regions  in  an  approximately  linear  fashion.  This 
phenomenon  is  known  as  thermel  conduction.  In  this  new  state 
that  the  system  has  reached  in  response  to  a  constraint,  stability 
will  prevail  again  and  the  behaviour  will  eventually  be  as  simple 
as  at  equilibrium.  H  owever,  by  removing  the  system  from  equi¬ 
librium  further  and  further,  through  an  increase  in  at,  we 
observe  that  suddenly,  at  a  value  of  at  that  we  will  call  critical, 
matter  begins  to  perform  a  bulk  movement.  Moreover,  this 
movement  is  far  from  random;  the  fluid  is  structured  in  a  series 
of  small  structures  known  asBenard  cells. 

0  wing  to  thermal  expansion,  the  fluid  closer  to  the  lower  plate 
is  characterised  by  a  lower  density  than  that  nearer  the  upper 
plate.  T  his  gives  rise  to  a  gradient  of  density  that  opposes  the 
force  of  gravity.  This  configuration  is  thus  potentially  unstable. 
Consider  a  small  volume  of  the  fluid  near  the  lower  plate. 
Imagine  that  this  volume  is  displaced  upward  by  a  perturba¬ 
tion.  This  volume,  now  in  a  colder  and  hence  denser  region, 
will  experience  an  upward  Archimedes  force,  amplifying  the 
ascending  movement  further.  If,  on  the  other  hand,  a  small 
droplet  initially  close  to  the  upper  plate  is  displaced  down¬ 
ward,  it  will  penetrate  an  environment  of  low  density  and  the 
Archimedes  force  will  tend  to  amplify  the  initial  descent.  T  he 
fluid  thus  generates  the  observed  currents.  The  stabilising 
effect  of  viscosity,  which  generates  an  internal  friction  oppos¬ 
ing  movement,  counteracts  the  destabilising  effects.  This,  and 
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thermal  conduction,  which  tends  to  average  any  temperature 
difference  between  the  displaced  droplet  and  its  environment, 
explains  why  currents  do  not  appear  as  soon  as  ^T  is  not 
strictly  zero. 

Benard  cells  also  show  the  complexity  of  movement.  T  he  cells 
unfold  along  the  horizontal  axis,  adopting  successively  right- 
handed  or  left-handed  rotation.  Our  very  small  observer  can 
now  locate  his  position  in  space  by  considering  the  rotation  of 
the  cell  he  occupies  and  by  counting  the  number  of  cells  he 
passes  through.  The  emergence  of  this  notion  of  space  is 
known  as  symmetry  bredking.  When  at  is  below  the  critical 
value,  the  homogeneity  of  the  fluid  in  the  horizontal  direction 
renders  its  different  parts  independent  of  each  other.  In  con¬ 
trast,  beyond  the  threshold,  it  is  as  if  each  volume  element  is 
watching  the  behaviour  of  its  neighbours  and  istaking  this  into 
account  in  order  to  play  its  role  adequately  and  to  participate 
in  the  overall  pattern.  T  his  suggests  the  existence  of  corrddtions 
of  statistically  reproducible  rate  relations  between  distant  parts 
of  the  system.  The  characteristic  space  dimension  of  a  Benard 
cell  is  in  the  millimetre  range,  whereas  the  characteristic  space 
scale  of  the  intermolecular  forces  is  in  the  Angstrom  range. 
That  large  numbers  of  particles  can  behave  in  a  coherent  fash¬ 
ion  at  this  long  range,  despite  random  thermal  motion,  is  one 
of  the  principal  properties  characteristic  of  such  self-organisation 
and  emergent  complex  behaviour. 

This  experiment  is  reproducible;  the  same  convection  patterns 
will  appear  at  the  same  threshold  value  and  the  process  is  sub¬ 
ject  to  a  strict  determinism.  However,  the  direction  of  the 
rotation  of  the  cells  is  unpredictable.  T  he  form  of  the  particu¬ 
lar  perturbation  that  prevails  at  the  moment  of  the  experiment 
will  decide  whether  a  given  cell  is  right-  or  left-handed.  When 
the  constraint  is  sufficiently  strong,  several  solutions  are  possi- 
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ble  for  the  same  parameter  values  and  chance  alone  will 
decide  which  of  these  solutions  is  realised.  I  n  this  way,  the  sys¬ 
tem  has  been  perturbed  from  a  state  of  equilibrium  or  near¬ 
equilibrium  to  a  state  of  sdf- or  gen  isati  on,  with  a  number  of  pos¬ 
sible  modes  of  behaviour. 

What  happens  to  the  Benard  cell  system  when  the  thermal 
constraint  is  increased  beyond  this  first  threshold?  For  some 
range  of  values  the  Benard  cells  will  be  maintained  globally 
but  some  of  their  specific  characteristics  will  be  modified.  Fur¬ 
ther  constraint  induces  the  system  to  move  beyond  another 
critical  point  and  turbulence  is  witnessed.  Note  that  all  of  these 
critical  behaviours  are  different  from  the  phase  changes  we 
normally  associate  with  closed  thermodynamic  systems.  The 
reason  behind  this  is  that  a  nonequilibrium  constraint  is  being 
applied.  For  example,  the  dendritic  structure  associated  with 
snowflakes  has  nothing  to  do  with  the  structure  of  the  underly¬ 
ing  ice-crystal  lattice.  The  scale,  size,  and  spacing  of  the 
emergent  structure  is  of  an  order  of  magnitude  larger. 

To  summarise,  nonequilibrium  has  enabled  the  system  to 
transform  part  of  the  energy  communicated  from  the  envi¬ 
ronment  into  an  ordered  behaviour  of  a  new  type:  the 
dissipative  structure.  This  regime  is  characterised  by  symmetry 
breaking,  muitipie  modes  of  behaviour,  and  correlation.  Such  a  system 
is  called  "open”  since  it  isopen  to  the  effect  of  energy  or  infor¬ 
mation  flowing  into  and  out  of  the  system.  It  is  also  called 
"dissipative”  because  of  such  energy  flows,  and  the  resultant 
dissipation  of  energy. 

Open,  Dissipative  Structures 

Consider  then  a  system  embedded  in  an  environment  with 
which  it  communicates  through  the  exchange  of  certain  prop- 
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ertiesthat  we  shall  call  fluxes.  We  now  know  to  call  this  an  open 
system  (in  contrast  to  a  closed  or  isolated  system).  Figure  1.2  is 
a  schematic  representation  of  such  an  open  system,  communi¬ 
cating  with  the  environment  through  the  exchange  of  such 
properties  as  mass,  energy,  or  information.  T  he  rate  of  amount 
transported  per  unit  surface  is  the  flux  of  the  corresponding 
property  across  the  system.  In  our  simple  example  above,  the 
amountofheatingisthefluxofenergyintothe  system . 


Figurel.2:  Schematic  of  an  Open  System 

As  a  result  of  these  exchanges,  the  variables  describing  the 
instantaneous  state  of  the  system,  {X,},  vary  in  time  and  attain 
values  typically  different  from  those  characterising  the  state  of 
the  environment  {X/e}.  Whatever  the  detailed  interpretation 
ofX/  might  be,  the  evolution  of  the  system  under  consideration 
may  be  described  in  the  following  general  form: 

Rate  of  change  of  the  system  state  =  function  of  (system  state 
and  control  variables). 

Thus  we  have: 


dt 
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where  f/  denote  the  laws  concerning  the  rate  of  change  of  the 
system,  and  are  a  set  of  parameters  present  in  the 

problem,  which  can  be  modified  by  the  external  world.  T  hese 
quantities  are  known  as  control  perem^ers.  U  nder  certain  condi¬ 
tions,  this  relation  will  have  a  single  solution  that  minimises 
some  measure  of  negative  utility  (which  we  call  a  loss  function). 
This  solution  is  then  a  unique  optlmel  control  for  the  system.  For 
each  time  t,  it  defines  an  optimal  value  for  the  settings  of  the 
control  parameters  . 

An  Exam  ple 

If  we  think  of  a  guided  missile  attempting  to  manoeuvre 
towards  a  target,  the  measure  of  loss  is  the  miss  distance  rela¬ 
tive  to  the  aim  point.  T  he  control  parameters  are  the  settings 
for  the  missile  fins  at  a  given  time  t.  For  simple  forms  of  llneer 
guidance  (e.g.,  early  forms  of  laser-guided  bombs),  this  leads  to 
what  is  called  beng-beng  control,  where  the  missile  fins  "bang” 
from  one  extreme  setting  to  another  in  order  to  keep  the  mis¬ 
sile  on  course.  The  Appendix  goes  into  this  in  more  detail  and 
shows  that  such  solutions  correspond  to  maximising  or  mini¬ 
mising  a  Flamiltonian  function.  This  is  due  to  Pontryagin's 
maximum  principle.  Applied  to  a  llneer  control  system,  this 
maximum  principle  leadsto  the  solution  of  bang-bang  control. 

A  characteristic  feature  of  many  of  the  systems  encountered  in 
nature,  however,  is  that  the  f 's  are  complicated  nonllneer  func¬ 
tions  of  the  X's.  The  equations  of  evolution  of  this  type  of 
system  should  then  admit,  under  certain  conditions,  several 
solutions  (rather  than  just  the  one  optimal  solution)  since  a 
multiplicity  of  solutions  is  the  most  typical  feature  of  a  nonlin¬ 
ear  equation.  Our  assumption  will  be  that  these  solutions 
represent  the  various  modes  ofbehevlour  of  the  underlying  system. 
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Modesof  Behaviour  in  Nonlinear  Systems 

How  do  these  different  modes  of  behaviour  arise  in  such  a 
nonlinear  system?  Thermodynamic  equilibrium  is  character¬ 
ised  by  detdiled bdlence,  i.e. 

Probabilityof  a  "direct"  process  =  Probabilityof  a  "reverse"  process 

We  can  understand  that  in  such  a  state,  any  attempt  at  diversi¬ 
fication  and  self-organisation  will  besmeared  out  immediately: 
equilibrium  is  a  state  of  full  homeostasis,  characterised  both  by 
uniqueness  and  robust  stability  properties.  Our  aim  is  to 
extend  our  ideas  of  equilibrium  to  the  nonequilibrium  dynam¬ 
ics  of  an  open,  dissipative  system. 

The  most  useful  view  of  equilibrium  is  as  follows.  We  represent 
the  evolution  of  the  system  in  a  space  spanned  by  the  state  vari¬ 
ables  (phdse  space).  An  instantaneous  state  of  the  system  is  thus 
represented  in  phase  space  by  a  point.  As  the  system  evolves 
over  time,  a  succession  of  such  states  is  produced,  giving  rise  to 
a  curve  in  phase  space,  which  iscalled  a  phase  space  trajectory.  In  a 
dissipative  dynamical  system,  as  time  progresses,  the  phase 
space  trajectory  will  tend  to  a  limit  representative  of  the  regime 
reached  by  the  system  when  all  transients  die  out.  We  call  this 
regime  the  attractor.  The  attractor  representing  an  equilibrium 
position  is  unique  and  describes  a  time-independent  situation. 
T  his  gives  a  phase  space  point  towards  which  all  possible  histo¬ 
ries  converge  monotonically.  The  state  of  equilibrium  is 
therefore  a  universal  point  attractor.  The  goal  of  self- organisation  is 
thus  the  search  for  new  attractors  that  arise  when  a  system  is 
driven  away  from  its  state  of  equilibrium.  For  example,  consider 
the  type  of  coupled  chemical  system  studied  by  Ilya  Prigogine 
[4]  (for  which  he  received  the  N  obel  Prize  in  C  hemistry): 
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A  +  B  _ 

^2 


3X 


H  ere  the  concentration  [x]  of  product  X  is  taken  to  be  the  only 
state  variable,  being  understood  that  A  and  B  are  continuously 
supplied  from  or  removed  to  the  outside  to  maintain  fixed  con¬ 
centrations.  At  equilibrium,  detailed  balance  implies  that  the 
rate  equation  is: 


A:j[a][x]^  =  k2{xf 
k2[x\  =  k^[b'\ 

These  relations  fix  the  equilibrium  value  [Xg]  of  x  uniquely  and 
impose  a  condition  on  the  concentrations  of  constituentsA  and  6: 

K  k^ 

[/)J  _  ^1^3 

[aj  k^k^ 

In  a  dissipative  stationary  state  far  from  equilibrium,  it  is  not 
necessary  for  each  individual  reaction  to  balance  in  both  direc¬ 
tions.  Cancelling  the  overall  effect  of  the  two  forward  reactions 
by  that  of  the  backward  reactions  is  sufficient,  and  this  yields  a 
cubic  equation: 

-^2  +K\.a]  [x^  ]^  +  ^3  [^j  -  ^4  [^]  =  0  ■ 


Open,  DissipativeStructures 


Chapter  1 


13 


T  his  equation  may  have  up  to  three  soiutions  (hence  three  pos- 
sibie  modes  of  behaviour)  for  certain  vaiues  of  [a]  and  [()].  it 
may  therefore  be  said  that  nonequiiibrium  reveaisthe  potenti- 
aiities  hidden  in  the  noniinearities,  which  remain  "dormant" 
at  or  near  equiiibrium. 

T  he  monotonic  character  of  the  approach  to  the  state  of  equi¬ 
iibrium  impiies  that  the  evoiution  iaws  of  the  system  shouid 
obey,  in  the  neighbourhood  of  eguiiibrium,  very  particuiar  condi¬ 
tions.  introducing  the  deviations  of  X/  from  the  equiiibrium 

vaiues,X/e- 

the  evoiution  of  {x,}  near  equiiibrium  can  then  be  written: 

_  V'  p 

dt  j 

o  is  a  thermodynamic  potentiai  taking  its  minimum  at  equi¬ 
iibrium  and  {r.J  is  a  symmetric  matrix.  This  symmetry  can  be 
traced  back  to  the  property  of  detaiied  baiance  or  the  property 
of  the  invariance  of  the  equiiibrium  state  to  time  reversai.  We 
can  see  from  this  that  the  system  response  is  essentiaiiy  i inear 
near  to  equiiibrium  (i.e.,  smaii  changes  iead  to  smaii  effects). 

The  Far-from-Equi librium  State 

The  search  for  a  generaiised  thermodynamic  potentiai  o  in 
the  noniinear  range  weii  away  from  the  equiiibrium  state  has 
attracted  a  great  deai  of  attention,  but  these  efforts  have,  so 
far,  not  made  much  progress.  Typicaiiy,  therefore,  beyond  the 
i inear  domain  for  such  irreversibie  processes,  the  above  con- 
troi  equations  are  expected  to  breakdown.  A  first  consequence 
is  that  the  steady  state  point  ettrector  of  the  system  (i.e.,  the 
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point  to  which  the  system  state  moves  as  the  system  evolves 
over  time),  extrapolating  the  state  of  equilibrium  as  the  dis¬ 
tance  from  equilibrium  is  increased,  can  now  be  approached 
through  damped  oscillations.  This  behaviour  heralds  a  still 
more  interesting  possibility  in  which  the  oscillation  eventually 
becomes  sustained.  Topologically,  this  implies  the  emergence 
of  a  new  one-dimensional  attractor  in  phase  space.  T  he  point 
attractor  for  the  state  of  the  system  is  essentially  replaced  by  a 
circle.  I  n  the  limit,  the  system  moves  endlessly  around  this  cir¬ 
cle,  which  is  thus  known  as  a  limit  cycle. 

By  allowing  the  intrinsic  nonlinearity  to  be  manifested  in  the 
regime  of  detailed  balance,  nonequilibrium  can  also  lead  to 
the  coexistence  of  multiple  attractors  in  state  space.  T  he  state 
space  can  then  be  carved  up  into  a  set  of  basins.  Each  of  these 
corresponds  to  the  set  of  states  that,  if  the  system  were  to  start 
from  there,  would  evolve  to  a  particular  attractor.  These  are 
known  as  the  basins  of  attraction.  The  ridges  separating  these 
basins  of  attraction  are  called  separatrices.  The  coexistence  of 
multiple  attractors  constitutes  the  natural  mode  of  systems 
capable  of  showing  adapted  behaviour  and  of  performing 
regulatory  tasks.  We  would  thus  expect  to  see  the  system  stay¬ 
ing  within  one  basin  of  attraction  (corresponding  to  resistance 
to  change)  and  then  at  some  point  switching  between  differ¬ 
ent  attractors  (corresponding  to  a  change  in  the  long-term 
mode  of  behaviour)  as  we  further  vary  the  initial  state  of  the 
system.  The  existence  of  one-dimensional  attractors  (points 
and  circles)  suggests  the  possibility  of  higher  dimensional 
attracting  objects  in  phase  space.  These  model  multiperiodic 
and  chaotic  behaviour,  which  is  observed  under  appropriate 
experimental  conditions. 

Nonequilibrium  phenomena  show  a  variety  of  behaviours  and 
therefore  correspond  to  the  movements  of  the  system  towards 
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different  attractors.  The  simplest  mechanism  to  depict  this  is 
the  bifurcation  diagram,  where  a  single  control  parameter  (the 
thermal  gradient  in  the  case  of  Benard  cells)  affects  the  dynam¬ 
ics  of  the  system. 


Bifurcations  AND  Universality  Limits 

Consider  a  system  described  by  a  set  of  evolution  laws  of  the 
form  of  the  equation: 


dt 


(/  =  !,. ...n) 


where,  as  before,  f  /  are  the  rate  laws,  and  X,  are  the  control 
parameters.  In  a  typical  natural  phenomenon,  the  number  of 
variables  n  is  expected  to  be  very  high.  T  his  will  considerably 
complicate  the  search  for  all  possible  solutions.  However,  sup¬ 
pose  that  by  experiment  we  know  one  of  the  solutions.  By  a 
standard  method,  known  as LinedrStebilityAndlysis,  we  can  then 
determine  the  parameter  values  for  which  this  solution 
regarded  as  the  reference  state  switches  from  asymptotic  sta¬ 
bility  to  instability. 


Robustness  AND  Linear  Stability  Analysis 

stability  or  "robustness  to  change”  is  essentially  determined  by 
the  response  of  the  system  to  perturbations.  It  is  therefore  nat¬ 
ural  to  transform  the  dynamical  laws  into  a  form  in  which  the 
perturbations  appear  explicitly: 

x,(0  =  x,,+x,(0 

and 

dx 

dt 
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T  hese  equations  are  homogeneous  in  the  sense  that  the  right 
hand  side  vanishes  if  all  X/  =  0.  Expanding,  we  may  write: 

dX:  ^ 

=  S  (  {^y  }  ’  ^)  '  =  1^2, n 

where  L/j  are  the  coefficients  of  the  linear  part  and  hj  are  the 
nonlinear  part.  It  is  assumed  that  the  asymptotic  stability  of 
the  reference  state  (i.e.,X  =X5orx  =  0)  of  the  system  is  iden¬ 
tical  to  that  of  the  linearised  part: 

dX:  ^ 

This  is  reasonable  provided  that  the  perturbation  is  not  too 
large  and  the  system  is  "well  behaved.” 

I  n  general,  a  multivariate  system  gives  rise  to  a  wide  spectrum 
of  values  for  the  rate  of  change  of  perturbations.  Fora  unique 
control  variable  'k  =  'k,,  two  cases  can  be  distinguished.  First, 
the  perturbations  are  nonoscillatory  and  the  bifurcations  (i.e., 
alternative  modes  of  behaviour)  will  correspond  to  steady- 
state  point  attractors,  or  secondly,  the  perturbations  are  oscil¬ 
latory  and  the  bifurcations  will  correspond  to  time-periodic 
solutions  in  the  form  of  limit  cycles.  M  ore  intricate  solutions 
can  also  be  envisaged  leading  to  secondary,  tertiary,  or  higher 
order  interactions;  however,  the  complete  stable  unfolding  of 
the  problem  remains  an  open  question.  It  is  left  to  other 
methods,  such  as  explicit  simulation  using  cellular  automata 
(i.e.,  agent-based  simulation)  to  attempt  to  solve  the  problem 
of  describing  the  dynamic  behaviour  of  nonequilibrium  sys¬ 
tems.  The  use  of  such  a  simulation  approach  can  be  described 
as  "experimental  mathematics." 
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Self-Organisation  in  Nature-An  Example 

Consider  an  ecosystem  consisting  of  a  large  number  of  inter¬ 
acting  species,  each  evolving  in  response  to  the  environment 
created  by  the  rest  of  the  ecosystem  (i.e.,  each  species  iscoei/o/i/- 
ing).  Such  a  system  consists  of  many  components  that  interact 
through  some  kind  of  exchange  of  forces  or  information.  In 
addition  to  the  internal  interactions,  the  system  may  be  driven 
by  some  external  force- natural  selection  in  this  case.  T  he  sys¬ 
tem  will  now  evolve  over  time  under  the  influence  of  the 
external  driving  forces  and  the  internal  interactions.  What 
happens  when  we  observe  such  a  system?  Is  there  some  simpli¬ 
fying  mechanism  that  produces  a  typical  behaviour  shared  by 
large  classes  of  such  systems?  T  he  mechanism,  it  turns  out,  is 
that  of  clustering. 

A  simple  cellular  automaton  model  of  such  an  ecosystem  is 
given  by  theBak-Sneppen  evolution  model.  It  is  an  example  of 
the  experimental  mathematics  we  described  earlier,  as  a  way 
of  analysing  the  coevolution  of  such  a  complex  system  over 
time.  A  description  of  this  model  of  coevolution  within  an  eco¬ 
system  was  given  in  the  last  chapter  of  [5]  as  part  of  an 
introduction  to  ideas  from  complexity  mathematics  and  the 
development  of  mathematical  "metamodels”  of  future  Infor¬ 
mation  Age  conflict.  For  ease  of  reference,  some  of  that 
description  is  repeated  here. 

The  Bak-Sneppen  Evolution  Model 

In  this  automaton  model,  we  have  a  d-dimensional  lattice  (Fig¬ 
ure  1.3)  and  random  numbers  f/  drawn  without  replacement 
from  the  interval  [0,1]  occupy  the  lattice  sites.  At  each  update 
step,  the  extremal  site  (that  is,  the  one  with  the  smallest  value  of 
fj)  is  chosen,  and  then  it  and  its  2d  immediate  neighbours  are 
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assigned  new  random  numbers.  As  a  model  of  evolution,  the 
values  f/  correspond  to  "fitness”  values.  Changing  both  the  site 
and  neighbouring  sites  captures  the  process  of  local  coevolution. 
It  follows  from  [6]  that  the  set  of  such  active  sites  is  a  fractal  in 
space-time  (see  particularly  Figures  1  and  28  of  that  reference). 


Figure  1.3:  T  heLatticeof  Species  Interactions  in  a  M  odd  Ecosystem 

As  described  by  [5,  Chapter  6],  the  approach  to  the  critical 
attractor  of  the  process  (at  which  avalanches/  clusters  of  all 
sizes  are  possible)  iscontrolled  by  the  "gap  equation": 

ciG(s)_  l-G(s) 

* 

where  Gfsj  is  the  maximum  extremal  value  fj(s)  at  times,  cor¬ 
responding  to  the  gap  opened  up  between  the  existing  state  of 
the  system  and  the  zero  or  ground  state,  L  is  the  linear  size  of 
the  lattice,  and  (s  is  the  average  avalanche  size  at  time- 
step  s.  (An  avalanche  or  cluster  consistsof  a  set  of  extremal  val¬ 
ues  f/,  each  of  which  is  a  neighbour  of  the  previous  extremal 
value.  T  he  size  of  an  avalanche  is  the  number  of  timestepsfor 
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which  this  process  continues.)  Schematicaiiy,  we  thus  have  the 
foiiowing  picture  of  the  process  (Figure  1.4): 


Figurel.4:  M  ovement  of  tfie E  cosystem  towards  a  5elf-0 rganised  C ritical  Point 


Figure  1.4  shows  how  the  equation  (the  hatched  iine)  approxi¬ 
mates  the  seif-organised  movement  of  the  system,  via  a  series 
of  avaianches/  ciusters  towards  the  criticai  attractor  of  the  sys¬ 
tem  at  which  the  system  has  optimai  fiexibiiity  (in  the  sense 
that  ciusters  of  aii  sizes  can  be  created).  T  his  criticai  point  cor¬ 
responds  to  a  fitness  vaiue  f^.  At  this  point,  there  are  no  fitness 
vaiuesbeiow  this  criticai  vaiue  and  a  fiat  distribution  of  fitness 
vaiues  in  the  range  from  to  1.0.  T  his  is  in  compiete  contrast 
to  the  behaviour  of  a  cioski  system  such  as  an  ideai  gas  in  an 
isoiated  container,  where  the  gas  evoives  from  (for  exampie) 
being  partitioned  in  part  of  the  container  to  the  equiiibrium 
state  where  it  is  spread  equaiiy  throughout. 

Such  criticai  systems  are  of  particuiar  scientific  interest.  Sys¬ 
tems  in  criticai  states  do  not  have  any  characteristic  scaie  and 
may  therefore  exhibit  the  fuii  range  of  behaviourai  characteris¬ 
tics  within  the  particuiar  system  restraints.  This  means  that 
systems  at  the  point  of  criticaiity  are  in  a  position  of  optimai 
fiexibiiity  in  some  sense,  as  we  have  noted,  it  couid  thus  be 
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argued  [5]  that  one  of  the  requirements  of  military  command 
is  to  so  arrange  things  that  the  forces  collaborate  locally  and 
thusself-organise  into  this  optimal  state. 

From  the  previous  equation,  the  rate  of  change  of  the  gap  is 
inversely  proportional  to  the  average  avalanche  size: 


G(5)oc 


1 


T hus at  the  critical  value  °°  and  near  to  critical¬ 

ity,  the  average  avalanche  size  satisfies  the  scaling  law: 


(5)  u  ii-fy 

for  some  exponent  y  (as  has  been  confirmed  by  experiment). 

Before  the  critical  point  is  reached,  at  some  time  f,  if  fo  is  the 
smallest  random  number  on  the  lattice  in  the  evolution  model, 
then  random  numbers  created  at  the  next  time-step  will  only 
continue  the  avalanche  process  if  they  are  smaller  than  fg. 
T  hus  the  value  fo  can  be  viewed  as  the  branching  probability 
of  a  random  process  over  time.  T  his  will  give  information  on 
the  avalanche  size.  If  fo  =  branching  probability,  then  for 
larger  fo  we  have  larger  avalanches.  We  thus  assume  [6]  a  scal¬ 
ing  relation  of  the  form: 

p(s,fj  =  s-^g(s(f^-f,)>"n 

to  describe  the  probability  distribution  of  avalanches/  clusters 
of  sizes  corresponding  to  an  extremal  value  fo.  Such  a  relation 
has  also  been  confirmed  by  simulation  experiments.  In  partic¬ 
ular,  when  fo  equals  f^,  the  critical  value,  the  probability  of  a 
cluster  of  sizes  is  given  by  P{S,fj  u  5  ^  We  can  think  of  the 
function  g  as  a  'cutoff'  corresponding  to  the  fact  that  the  ava- 
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lanches  created  will  have  a  finite  size.  As  we  move  to  the 
critical  state,  this  cutoff  dies  away.  We  shall  call  r  the  clustering 
exponent  for  this  coevolving  ecosystem.  In  the  context  of 
manoeuvre  warfare,  this  describes  the  statistics  of  local  cluster¬ 
ing/  collaboration  at  a  transient  point  fo  heading  towards  the 
critical  attractor  value  f^.  T  he  parameters  r  and  cr  are  model 
dependent  and  g  is  our  "scaling  function.”  Note  again  that  the 
average  size  of  the  fo  avalanche  diverges  as  f^,  i.e. 

(5)  u  a-y;r. 

If  we  mark  each  of  the  minimal  sites  on  the  lattice  as  it  is  iden¬ 
tified  as  an  extremal  value  fg,  then  the  set  of  marks  generated 
over  time  forms  a  fractal  in  space-time  [6]  as  we  have  already 
noted.  Cuts  of  this  fractal  in  the  space  direction  at  a  given  time 
identify  the  site  that  is  "active"  (i.e.,  chosen  as  the  minimal  site) 
at  that  time.  Cuts  in  the  time  direction  produce  a  fractal  time 
series.  I  n  C  hapter  3,  we  show  that  the  time  series  of  casualties 
in  conflict  has  the  character! sties,  in  some  cases,  of  such  a  frac¬ 
tal  time  series.  0  ne  interpretation  of  such  effects  isthat  it  isthe 
dynamics  of  local  clustering  (by  one  side)  that  is  leading  to 
casualties  to  the  other  side.  This  can  also  be  shown  to  occur  in 
cellular  automata  modelsof  conflict  such  as  the  ISAAC  model 
produced  by  the  U.S.  Marine  Corps  Combat  Development 
Centre,  asdiscussed  in  Chapter  6  of  reference  [5]. 

Self-Organised  Criticality 

In  a  paper  published  in  1987,  Bak,  Tang  and  Weisenfeld  [7] 
first  proposed  the  hypothesis  that  a  system  consisting  of  many 
interacting  constituents  may  exhibit,  in  certain  cases,  a  specific 
general  emergent  behaviour  characteristic  of  the  system.  Bak 
described  the  behaviour  of  this  type  of  system  by  the  term  seff- 
orgenised  criticdlity  (SOC).  Self-organisation  has  for  many  years 


Self-Organisation  in  Nature- An  Exampie 


22 


C  omplexity  T  heory  and  N  etwork  C  entric  W  arfare 


been  used  to  describe  the  ability  of  certain  nonequilibrium  sys¬ 
tems  to  develop  structure  and  patterns.  T  he  word  criticdlity  has 
a  very  precise  meaning  in  thermodynamics.  It  is  used  in  con¬ 
nection  with  phase  transitions.  At  all  temperatures  other  than 
the  transition  temperature,  perturbations  of  the  system  will 
only  locally  influence  system  components  At  the  critical  tem¬ 
perature,  the  perturbation  affects  the  whole  system,  even 
though  only  the  nearest  neighbour  system  components  inter¬ 
act  directly.  The  system  becomes  critical  in  the  sense  that  all  of 
the  members  of  the  entire  system  influence  each  other.  For  the 
example  ecosystem  above,  the  system  self-organises  itself  into 
the  critical  state  corresponding  to  this  ability  of  the  entire  sys¬ 
tem  to  be  influenced  through  the  propagation  of  local 
coevolution  influences  and  the  resultant  clusters/  avalanches  of 
species  that  coevolution  created. 

Separation  of  Internal  and  External 
T IMESCALES 

Such  a  self-organised  dynamical  state  requires  the  separation 
of  the  external  and  internal  timescales.  For  example,  the  stress 
in  the  earth's  crust,  built  up  over  a  period  of  time,  isof  a  differ¬ 
ent  scale  to  the  subsequent  earthquake  lasting  merely  minutes 
or  seconds.  The  force  applied  to  an  individual  tectonic  plate 
must  overcome  a  threshold  in  order  to  produce  an  earthquake. 
T  his  means  that  the  plate  may  exist  in  a  multitude  of  interme¬ 
diate  states- mefasfaWesfafes-on  the  way  to  criticality.  Among  all 
of  the  metastable  states,  some  are  of  particular  importance. 
T  hese  states  are  merginelly  stable;  a  slight  disturbance  may  lead 
to  a  wide  variety  of  responses.  Bak  envisaged  that  these  mar¬ 
ginally  stable  states  are  characterised  by  the  lack  of  any  typical 
time  or  length.  T  his  configuration  is  of  a  similar  type  to  that 
seen  in  thethermodynamicsof  phase  changes.  The  lack  oftyp- 
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ical  scale  leads  to  algebraic  correlation  functions  associated 
with  power-laws'^  Hence,  the  distribution  functions  describing 
the  frequency  with  which  SOC  events  occur  exhibit  power- 
law  characteristics.  For  our  earthquake  example,  if  E  is  the 
energy  released  during  an  earthquake,  then  the  probability  of 
an  earthquake  of  that  size  is  given  by  the  power-law  relation¬ 
ship  p{E)-e  '^ .  As  noted  in  [8],  such  a  simple  law  should 
have  an  elegant  explanation.  We  shall  see  later,  in  Chapter  3 
when  examining  the  evidence  for  such  emergent  behaviour  in 
warfare,  that  the  number  of  casualties  in  war  also  has  such  a 
power-law  distribution.  However,  as  yet  no  such  elegant 
explanation  is  available,  either  for  earthquakes  or  wars. 

Clustering  in  Space  and  Time 

The  formation  of  clusters  that  are  frectel  in  both  space  and 
time  is  common  in  natural  systems,  as  we  have  already  seen.  It 
wasthistype  of  behaviour  that  Bak  wanted  to  explain.  We  will 
see  that  these  relate  to  ideas  of  corrddtion  in  space  or  time  (in 
contrast  to  coincidence  in  space  or  time).  Correlation  in  space 
or  time  is  a  signal  of  local  clustering  and  collaboration  spatially 
(e.g.,  across  a  battlespace)  or  in  time  (e.g.,  across  an  informa¬ 
tion  grid-reading  e-mail  creates  a  correlation  in  time  between 
individuals,  taking  a  phone  call  creates  a  coincidence  in  time). 
T  he  properties  of  fractals  and  their  link  to  chaotic  behaviour 
have  been  examined  intensively  over  the  last  two  decades.^ 


power-law  / (x)  =  has  the  property  that  the  relative  change 

f(^)  ,  a 

- =  k 

m 

is  independent  ofx.  Power-laws  in  this  sense  lack  characteristic  scale. 

^See  for  example  the  standard  text  C haos  and  Fractals  by  H  einz-0  tto 
Peitgen,  H  artmut  Jurgens  and  Dietmar  Sauape.  Springer.  1992. 
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Despite  this*  very  little  (until  now)  has  been  known  about  why 
fractals  form.  Fractal  structures  are  not  the  lowest  energy  con¬ 
figuration  that  can  be  selected  in,  for  example,  thermo¬ 
dynamic  systems,  therefore  some  kind  of  dynamic  selection  of 
configuration  must  betaking  place. 

Bak  explains  the  connection  in  the  following  way.  A  signal  will 
be  able  to  evolve  through  the  system  as  long  as  it  is  able  to  find 
a  connected  path  of  above-threshold  regions.  When  the  system 
is  either  driven  at  random  or  started  from  a  random  initial 
state,  regions  that  are  able  to  transmit  a  signal  will  form  some 
kind  of  random  network.  This  network  is  correlated  by  the 
interaction  of  the  internal  dynamics  with  the  external  field. 
The  complicated  interrelation  between  the  two  driving 
dynamics  means  that  a  complex,  finely-balanced  system  is  pro¬ 
duced.  As  the  system  is  driven,  after  this  marginally  stable  self- 
organised  state  has  been  reached,  we  will  see  flashes  of  activity 
as  external  perturbations  interact  with  internal  drivers  to  spark 
off  avalanches  (i.e.,  clusters)  of  activity  through  different  routes 
in  the  system.  Bak's  assertion  is  that  the  structure  of  this 
dynamic  network  is  fractal.  If  the  activated  clusters  consist  of 
fractals  of  different  sizes,  then  the  duration  of  the  induced  pro¬ 
cesses  travelling  through  these  fractals  will  vary  greatly. 
Different  timescales  of  this  type  lead  to  what  is  termed  1/  f  noise. 

1/  f  noise  is  a  label  used  to  describe  a  particular  form  of  time 
correldtion  in  nature.  If  a  time  signal  fluctuates  in  a  seemingly 
erratic  way,  the  question  is  whether  the  value  of  the  signal 
N(Tg)  at  time  has  any  correlation  to  the  signal  measured  at 
time  To -rr  (w(To -rx)).  The  amount  of  causation  is  character¬ 
ised  by  a  temporal  correlation  function: 

G(t)  =  (N(x,  )N(x,  +  -  {N(x„))l . 
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The  correlation  function,  for  a  stationary  process,  is  linked  to 
the  power  spectrum  as  follows: 


S{f)  =  2j  JTG(x)cos(27t/c) 

0 

where  the  power  spectrum  is  defined  in  terms  of  the  square 
amplitude  of  the  Fourier  transform  of  the  time  signal: 


1 

T 


2 


1/f  noise  corresponds  to  the  case  where  s{f)  =  y  and  corre¬ 
sponds  to  a  fractal  clustering  of  the  signal  amplitude  in  time. 


Fractal  Clustering  in  Space 

For  a  system  that  is  not  poised  in  a  critical  state  and  thus  not 
about  to  change  its  mode  of  behaviour,  the  reaction  of  the  sys¬ 
tem  is  described  by  a  characteristic  response  time  and 
characteristic  length  of  scale  over  which  the  perturbation  is  felt. 
Flowever,  for  a  critical  system,  the  same  perturbation  applied 
at  different  positions  or  the  same  position  at  different  times  can 
lead  to  a  response  of  any  size.  T  he  average  is  not  therefore  a 
useful  measure  of  response.  T  he  amount  of  the  system  involved 
is  a  cluster  in  the  spatial  dimensions  of  the  system. 

An  Example-A  Pile  of  Dry  Sand  on  a  Beach 

T  he  nature  of  this  critical  state  can  be  illustrated  using  a  sim¬ 
ple  example-a  pile  of  sand.  Anyone  can  try  this  experiment 
on  a  sunny  day  by  the  seaside.  Normally,  the  grains  of  sand 
only  interact  locally  and  nothing  surprising  happens.  How- 
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ever,  if  the  critical  state  (defined  by  the  slope  of  the  pile)  is 
perturbed  by  the  addition  of  one  grain  of  sand  to  a  randomly 
chosen  position  on  top  of  the  pile,  the  extra  grain  will  induce 
an  avalanche  (a  cluster)  characterised  by  spatial  and  temporal 
characteristics  such  as  the  total  number  of  sand  grains  s 
involved  in  the  avalanche  and  the  lifetime  of  the  avalanche  t 
The  statistical  distributions  describing  the  response  are 
denoted  byPfsj  and  P(t).  In  the  critical  state,  we  expect  broad 
power-law  distributions  of  the  form  P{s)-s^^  and  P(0~r“. 
The  particular  values  of  aandjB  are  then  characteristic  of 
how  the  system  can  create  such  correlations  in  time  and 
space.  T  hese  distributions  will  typically  be  bounded  by  spetidl 
end  temporel  lower  end  upper  cut-offs.  For  example,  an  avalanche 
cannot  involve  the  displacement  of  less  than  onegrain  of  sand, 
and  the  duration  of  an  avalanche  cannot  be  shorter  than  the 
time  it  takes  one  grain  to  move  a  distance  equal  to  the  size  of  a 
single  grain. 

In  the  thermodynamic  theory  of  phenomena,  the  spatial  corre¬ 
lation  function  plays  a  fundamental  role.  If  the  system  is 
described  by  a  space-time  field  n{r,t),  the  spatial  correlation 
function  is  defined  as: 


G(r)  =  (n(ro)n(r,  +r))^^  -(<^o))\ 

where  a  thermal  average  and  an  average  over  the  position 
is  assumed. 

Away  from  the  critical  temperature,  the  correlations  decay 
exponentially,  G(r)~exp(-r/^),  beyond  the  correlation 
length  The  correlation  length  diverges,  as  the 

critical  temperature  is  approached.  At  the  critical  point, 
T  =  T^,  the  correlation  function  changes  functional  behaviour 
from  exponential  to  algebraic  dependence  upon  r,  i.e. 
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Gir)-r-\J  he  divergence  is  considered  to  bethesignai  of  the 
iack  of  a  characteristic  iength  scaie.  Seif-organising  systems 
may  be  characterised  and  examined  in  this  way. 

Cellular  Automata 

The  introduction  in  1987  of  the  Seif-0  rganised  Criticaiity  con¬ 
cept  empioyed  the  ianguage  of  avaianches  (ciusters)  of  sand 
grains  it  proposed  a  simuiation  modei  (a  cdluler  eutometd  modei 
to  be  precise)  of  the  most  essentiai  features  of  sand  dynamics. 
This  ceiiuiar  automata  modei  is  indeed  characterised  by 
power-iaws  and  exhibits  criticai  behaviour.  As  an  exampie,  iet 
us  iook  again  at  the  modeiiing  of  an  ecosystem  consisting  of  a 
number  of  coevoiving  species  (the  Bak-Sneppen  evoiution 
modei)  that  we  mentioned  eariier.  This  automaton  process 
considers  the  points  of  a  grid  and  has  a  simpie  set  of  ruies 
determining  how  the  system  changes  from  onetime-step  to  the 
next,  to  represent  the  changing  fitness  of  the  species  itseif,  and 
the  coevoiutionary  impact  of  that  change  on  the  fitness  of 
cioseiy  iinked  species  in  the  ecosystem.  N  ote  that  this  iinkage  is 
one  of  iocai  species  influence  and  coevoiution.  it  does  not  nec- 
essariiy  assume  physicai  cioseness.  Aithough  these  ruies  are 
simpie,  the  emergent  behaviour  of  the  system  is  compiex  and 
surprising-a  characteristic  of  such  noniinear  interactions. 

An  Example:  C lust ERiNG  and  CoEVOLUTiON  on 
A  Global  Inform  at  iON  Grid 

We  can  use  the  evolution  model  just  described  to  gain  insight 
into  the  effect  of  a  G  lobal  I  nformation  G  rid.  I  magine  such  a 
grid  in  two  dimensions.  At  each  grid  point  is  positioned  an 
element  of  our  force.  Each  such  force  element  has  a  "fitness" 
value  corresponding  to  its  ability  to  evolve  and  adapt  to  local 
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circumstances  as  a  function  of  the  information  available  on 
the  grid.  We  assume  these  fitness  values  are  random  at  first. 
At  each  step  of  the  process,  we  assume  that  the  force  element 
with  the  smallest  fitness  is  likely  to  have  to  adapt  fastest  to  its 
local  environment.  In  so  doing,  it  will  change  the  fitness  val¬ 
ues  of  the  units  closest  to  it  on  the  information  grid  (i.e.,  there 
is  local  coevolution).  Note  that  these  force  elements  may  be 
separated  by  large  and  varying  distances  in  space.  With  these 
assumptions,  over  time  the  force  elements  will  form  clusters  of 
coevolution  oftheform  predicted  by  the  Bak-Sneppen  model. 

I  n  particular,  the  statistics  of  emergent  cluster  size  can  be  pre¬ 
dicted  mathematically  to  converge  to  a  power-law.  Thus,  the 
probability  of  a  cluster  of  size  S,  P(S),  is  of  the  form 
where  r  is  a  characteristic  cluster  exponent  for  this  particular 
information  grid. 

T  hus  far,  we  have  looked  at  complex  effects  and  correlations 
across  space  and  time.  What  happens  if  we  restrict  ourselves  to 
looking  at  the  boundary  between  two  different  regimes  (such 
as  two  different  nationalities  or  two  opposing  armed  forces), 
and  how  this  would  move  over  time  depending  on  the  local 
coevolution  of  the  elements  involved? 

Movement  of  a  Boundary 

I  n  natural  systems,  we  can  consider  the  movement  of  a  bound¬ 
ary  through  a  medium  (for  example,  the  boundary  of  an 
atomic  surface,  the  boundary  of  a  growing  cluster  of  bacteria, 
or  the  front  of  advance  of  a  fluid  "invasion”  of  a  medium  such 
as  a  crystalline  rock).  T  his  has  been  studied  extensively  in  rela¬ 
tion  to  the  laying  down  of  single  atom  surfaces  using  molecular 
beam  epitaxy  [9].  The  most  relevant  case  from  our  point  of 
view  is  the  front  of  advance  of  fluid  "invasion”  of  a  medium. 
As  described  in  [9],  we  can  represent  the  medium  itself  as  con- 
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sisting  of  a  lattice  of  cells,  each  with  either  a  1  or  0  in  it.  A  "1" 
represents  the  fact  that  that  cell  can  be  wetted.  T  he  proportion 
of  cells  containing  a  "1”  is  defined  as  p.  For  large  configura¬ 
tions,  we  can  also  interpret  p  as  the  probability  that  a 
particular  cell  containsa  "1.”  A  "0”  represents  the  fact  that  the 
cell  cannot  be  wetted-it  thus  "pins”  the  advance  of  the  fluid 
through  the  medium,  at  least  locally.  I  n  the  standard  D  irection 
Percolation  Depinning  (DPD)  model,  we  start  off  with  (for 
example)  a  two-dimensional  square  lattice  of  cells,  some 
labelled  with  Is  and  the  rest  with  Os  to  represent  the  distribu¬ 
tion  of  these  "pinning”  forces  throughout  the  medium. 
Initially,  we  wet  one  edge  of  the  lattice.  In  Figure  1.5,  we  show 
the  wetted  right-hand  edge.  In  the  standard  DPD  model,  one 
of  the  unpinned  cells  in  the  next  column  is  chosen  at  random, 
and  the  fluid  invades  that  cell. 
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Figurel.5:  Local  Pinning  of  a  Fluid  Boundary 
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Itturnsout  that  for  this  case,  when  the  pinning  probabiiityp  is 
greater  than  a  criticai  vaiue  p^,  the  growth  of  the  interface  is 
haited  by  a  spanning  path  of  pinning  ceiis.  Such  modeis  of 
interface  or  boundary  movement  exhibit  fractai  properties  of 
the  interface,  as  discussed  in  detaii  in  [9],  Weshaii  seesimiiar 
effects  iater  in  our  discussion  in  Chapter  4  of  the  controi  of  the 
battiespace  using  ideas  based  on  preventing  the  fiow  of  oppos¬ 
ing  forces  and/ or  third  parties  through  the  space.  Rather  than 
choosing  the  next  ceii  to  invade  at  random,  as  in  the  DPD 
modei,  we  can  use  a  modei  of  the  process  that  is  more  akin  to 
the  manoeuverist  principie  of  appiying  your  strength  where 
the  opponent  is  weak-in  other  words,  the  ceii  next  to  be  wet¬ 
ted  is  the  one  where  the  iocai  pinning  force  of  the  medium  is 
weakest.  Such  a  modei  of  the  boundary  movement  isthe/ni/a- 
sion  Percoldtion  model.  We  can  create  a  modei  of  this  process  in  a 
way  that  is  consistent  with  our  description  of  the  Bak-Sneppen 
evoiution  modei  of  iocai  coevoiution  [6].  We  start  by  assign¬ 
ing,  as  with  the  Bak-Sneppen  evoiution  modei,  random 
numbers  f/  between  0  and  1  to  the  points  of  a  d-dimensionai 
iattice.  initiaiiy,  one  side  of  the  iattice  is  the  wetted  ciuster. 
The  random  numbers  at  the  boundary  of  the  wetted  ciuster 
are  examined.  At  each  update  step  s,  the  site  with  the  smaiiest 
random  number  f/5  on  the  boundary  of  the  wetted  area  is 
iocated  and  added  to  the  ciuster.  i  n  this  case,  we  can  interpret 
the  vaiues  f/  as  the  vaiues  of  the  iocai  pinning  force,  and  the 
ciuster  advances  at  those  points  where  the  pinning  force  is 
smaiiest.  As  noted  in  [6],  an  important  physicai  reaiisation  of 
invasion  percoiation  is  the  dispiacement  of  one  fiuid  by 
another  in  a  porous  medium.  The  boundary  of  the  ciuster  cre¬ 
ated  by  this  process  isfractai  and  has  a  fractai  dimension  in  the 
range  1.33-1.89  dependent  on  the  exact  definition  of  the 
boundary  [6,  Appendix]. 
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An  Analytical  Modelofthis  Process 

In  natural  systems,  the  boundary  of  such  an  interface  that  is 
moving  through  a  medium  can  be  characterised  by  its  "rough¬ 
ness.”  This  is  defined  as  follows. 

We  assume  that  the  interface  is  defined  across  a  linear  space  of 
length  L  (as  shown  in  Figure  1.6).  T  he  roughness  of  the  interface 
at  a  given  timet  and  for  a  span  of  length!,  is  defined  [9]  as: 


where  the  interface  length  is  divided  into  a  number  of  natu¬ 
ral  cells  of  unit  length,  h(i,t)  is  the  height  of  the/th  column  at 
time  t  (as  indicated  in  Figure  1.6)  and  k(t)  is  the  average  of 
these  heights  at  time  t.  T  he  roughness  1/1/  (L,t)  is  thus  the  stan¬ 
dard  deviation  of  the  height  at  a  given  time.  For  many 
natural  systems,  the  roughness  first  goes  through  a  transition 
period  before  stabilising  at  an  equilibrium  value.  0  nee  it  has 
stabilised,  the  expected  behaviour  is  that  this  "saturation" 
value  1/1/ saf  scales  with  the  interface  width  L,  i.e.  we  expect 
the  relationship  1/1/ saf  u  l"  . 
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When  this  occurs,  the  exponent  a  is  referred  to  as  the  roughness 
exponent.  It  is  then  a  characteristic  exponent  of  the  invasion 
process  under  study.  A  typical  value  of  a  is  about  0.6  (as 
shown  in  [9]).  This  idea  can  also  be  used  to  characterise 
Brownian  motion,  as  we  discuss  later  in  this  chapter  in  the 
context  of  stock  price  dynamics. 

By  using  symmetry  arguments,  we  can  derive  [9]  an  analytical 
expression  of  the  rate  at  which  such  an  interface  moves 
through  a  medium.  (In  terms  of  conflict,  this  corresponds  to 
the  rate  of  advance  of  a  combat  front.)  This  is  given  by  the 
K  ardar-Parisi-Zhang  (K  PZ)  equation: 

+ -  ryhf + riix,  t) . 
dt  2 

T  he  first  term  in  this  equation  represents  linear  effects  of  the 
interface  growth,  the  second  captures  nonlinear  effects,  and 
thethird  isa  noise  term.  T  histhus  represents  the  starting  point 
for  an  analytical  expression  (i.e.,  a  metamodel)  of  the  advance 
of  a  conflict  front  through  a  locally  controlled  area,  as  dis¬ 
cussed  below. 

An  Exam ple:  Force  Control  along  a 
Boundary 

T  he  concept  of  pinning  a  fluid  locally  is  similar  to  the  idea  of 
trying  to  exert  local  control  over  a  boundary  to  prevent  the 
flow  of  other  forces  or  third  parties  across  that  boundary.  We 
will  see  later  (in  C  hapter  4)  that  the  idea  of  control  as  the  pre¬ 
vention  of  such  flows  through  an  area  has  important 
implicationsfor  the  emergent  behaviour  of  a  force  (or  two  com¬ 
peting  forces)  attempting  to  exert  control  over  a  battlespace. 
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I  n  the  case  of  control  along  a  boundary,  C  omplexity  T  heory 
in  terms  of  the  invasion  percolation  model,  can  be  used  to 
analyse  the  effect  of  two  forces  (an  attack  and  a  defence  force) 
interacting  across  a  boundary,  when  the  boundary  moves  at 
the  point  where  the  defending  (pinning)  force  is  weakest.  If  the 
defence  pinning  force  is  coevolving  locally,  then  the  boundary 
should  form  a  fractal  with  a  fractal  dimension  in  the  range 
1.33-1.89  [6,  Appendix],  as  we  have  seen.  In  Chapter  3,  we 
show  that  there  is  historical  evidence  in  warfare  for  such  an 
effect,  and  for  values  in  this  range. 

An  Example  of  Complex  Behaviour  and 
Fractal  Time  Series  in  Economics 

Brownian  Motion,  Fractals,  and  Similarity 

M  andelbrot  [10]  has  considered  the  movement  of  stock  prices 
using  fractal  ideas  and  we  start  by  using  a  simple  example 
based  on  Brownian  motion  (since  this  is  the  basis  of  most  cur¬ 
rent  predictors  of  stock  price  volatility).  Consider  then 
Brownian  motion  in  one  space  variable,  thus  the  motion  of 
particles  is  restricted  to  a  line.  T  he  impacts  affect  the  particle 
only  from  the  left  and  right,  causing  a  displacement  of  length  / 
in  either  direction.  Can  any  prediction  be  made  about  the 
total  displacement  after  a  number  of  time-steps  n?  First,  the 
total  expected  displacement  is  zero,  as  all  displacements  are  -i-/ 
or  -/,  both  with  equal  probability  0.5.  Consider,  instead,  the 
square  of  the  displacement.  T  he  average  of  these  square  dis¬ 
placements,  called  the  man  square  dispidcement,  indicates  how 
much  the  particles  spread  in  a  given  number  of  time-stepson 
the  average.  Its  value  is n/^.  In  terms  of  Brownian  motion,  the 
number  of  steps  corresponds  to  the  number  of  impacts  on  a 
particle  and  cannot  be  directly  measured  in  an  experiment. 


AnE  xample of  C omplex  B  ehaviour  and  F racial  T  ime Series  in  E  conomics 


34 


C  omplexity  T  heory  and  N  etwork  C  entric  W  arfare 


Consider,  therefore,  time  duration  t  Assuming  an  average 
number  ofn  impulses  during  a  timespan  f,  the  particle  travels  a 
total  length  n/.  If  1/  is  the  average  speed  of  the  particle,  then  we 
have  i/f  =  n/.  T  he  mean  square  displacement  nl^  thus  equals  vlt 
Therefore,  the  mean  square  displacement  is  proportional  to 
the  timespan.  Is  there  anything  else  that  can  be  said  regarding 
the  distribution  of  the  displacement  after  time  f?  Experimen¬ 
tally,  it  is  found  that  the  statistics  for  a  simulated  one¬ 
dimensional  Brownian  motion  follow  a  Gaussian  or  normal 
distribution.  This  is  not  surprising  as  a  normal  distribution 
arises  where  independent  and  identically  distributed  random 
events  are  averaged  and  is  an  implication  of  the  Central  Limit 
Theorem  of  Statistics. 

Brownian  motion  is  also  associated  with  scaling  rddtionships  and 
fractal  behaviour.  Asa  nonfractal  object  is  magnified,  no  new  fea¬ 
tures  are  revealed.  As  a  fractal  object  is  magnified,  finer  details 
are  revealed.  The  size  of  the  smallest  feature  of  a  nonfractal 
object  is  called  the  characteristic  scale.  A  measurement  made  at 
finer  resolution  will  include  more  of  these  smaller  pieces.  Thus 
the  value  measured  of  a  property  will  depend  upon  the  resolu¬ 
tion  used  to  make  the  measurement.  How  a  measured 
property  depends  on  the  resolution  used  to  make  the  measure¬ 
ment  is  called  the  scaling  relationship.  A  fractal  object  has 
features  over  a  broad  range  of  sizes.  Fractal  phenomenological 
characteristics  are: 

1.  Self-sim  ilarity:  behavioural  characteristics  are  "simi¬ 
lar"  at  different  resolutions. 

2.  Scaling:  the  value  measured  for  a  property  depends 
upon  the  resolution  at  which  it  is  measured. 

3.  Dimension:  thedimension  of  an  object  gives  a  quanti¬ 
tative  measure  of  self-similarity  and  scaling.  It  tells  us 
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how  many  new  pieces  of  an  object  are  revealed  as  it  is 
viewed  at  higher  magnification. 

4.  Nonstatistical  properties  may  be  observed. 

M  oments  may  be  zero  or  nonfinite  (e.g.,  the  mean  tends 
towards  zero  and  variance  tends  towards  infinity). 

T  here  are  two  types  of  self-similarity: 

1.  Geometrical:  pieces  of  the  object  are  exact  smaller 
copies  of  the  whole  object. 

2.  Statistical:  the  value  of  the  statistical  property  Ofrj 
measured  at  resolution  r  is  proportional  to  the  value  of 
Q(dr)  measured  at  a  resolution  ar  such  that  0  far)  =kQ(r). 
For  statistical  probability  distribution  functions  (pdf), 
this  implies  that:  pdflC^farj]  =  pdf[/cC)('rj]. 

Statistical  self-similarity  is  possible  in  both  space  (spatial)  and 
time  (temporal). 

In  summary  then,  such  self-similarity  implies  a  scaling  rela¬ 
tionship.  T  he  simplest  form  of  the  scaling  relationship  is  that 
the  measured  value  of  a  property  Q  (r)  depends  on  the  resolu¬ 
tion  used  to  make  the  measurement,  as  in  the  equation 
Q(r)  =Br‘’,  i.e.  LogQ(r)  =  LogB  +  bLog  r,  where  6  and  b  are  con¬ 
stants.  Hence  for  self-similarity,  we  observe  log-log  linear 
behaviour.  Experimentally,  this  implies  that  for  Brownian 
motion,  there  should  be  a  scaling  factor,  r,  that  yields  curves 
that  are  visibly  identical,  i.e.  when  scaling  Brownian  motion 
in  time,  by  a  factor  say  of  t  and  in  amplitude  by  a  factor  of  r, 
we  should  see  no  difference.  This  transformation  is  called  a 
scdling  collapse  since  it  has  the  effect  of  collapsing  the  curves  at 
different  scales  of  measurement  on  top  of  each  other  into  one 
normalised  relationship. 
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The  scaling  relationship  in  time  may  be  determined  theoreti¬ 
cally  for  Brownian  motion.  It  follows  from  the  mean  squared 
analysis  that  the  mean  squared  displacement,  of  the 
Brownian  motion  X(t),  isdescribed  by  oc^.  Now  consider 
the  rescaled  random  function: 


7(0  =  rX 


\aJ 


i.e.  the  graph  ofX  is  stretched  in  the  time  direction  by  a  fac¬ 
tor  a  and  in  the  amplitude  by  r.  The  displacements  in  y  for 
time  differences  t  are  the  same  as  those  in  X  multiplied  by  r 
for  corresponding  time  differences  f/a.  Thus,  the  squared 
displacements  are  proportional  to  r^t/a.  In  order  to  ensure 
the  same  constant  of  proportionality  as  the  original  Brown¬ 
ian  motion,  we  require  r^/a  =  \  or  r  =  ^fa.  For  example, 
when  replacing  t  by  t/  2,  i.e.  stretching  the  graph  by  a  factor 
of  2,  we  havea=2  i.e.,  r  =  V2  . 


In  this  broader  context  of  fractal  processes*  ordinary  Brownian 
motion  is  a  random  process/ ft)  with  Gaussian  increments  and 
var(X(^2)-X(^i))oc  1^2 where  H  (the  Hurst  expo¬ 
nent)  =  V2.  We  can  consider  (as  discussed  in  [9])  the  Brownian 
motion  time  series  as  describing  an  interface  (between  the 
parts  above  and  below  the  series)  stretching  between  the  time 
points  ti  and  (2.  In  terms  of  our  previous  discussion  of  interface 
roughness,  we  now  have  L  =  T  he  standard  deviation  (i.e., 
the  roughness  of  the  interface  generated  by  the  Brownian 
motion  time  series)  over  this  timespan  L  has  the  form  l" 
where  in  this  case  a  equals  the  H  urst  exponent  H  (as  we  can 
see  from  the  expression  for  the  variance  above  in  terms  of  H  ). 
Thus,  standard  Brownian  motion  corresponds  to  a  roughness 
exponent  of  V2.  0  ther  values  of  H  are  possible,  corresponding 
to  rougher  or  smoother  forms  of  time  series. 
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Stock  Pricesand  Brownian  Motion 

M  andelbrot  considered  the  movement  of  stock  prices  using  a 
Brownian  random  walk  process  as  his  starting  point  [10],  T  he 
behaviour  of  such  a  variable  z  can  be  understood  by  consider¬ 
ing  the  changes  in  its  value  in  small  intervalsof  time.  Consider 
such  a  small  interval  of  time  and  define  Az  as  the  change 
in  z  during  M.  tsz  must  have  two  basic  properties:  first, 
Az  =  £-Va7  where  £•  is  a  random  draw  from  a  standardised 
normal  distribution  (i.e.,  a  normal  distribution  with  mean  zero 
and  standard  deviation  of  1.0);  and  second,  the  values  of  Az 
for  any  two  different  short  intervals  of  time  A^  are  indepen¬ 
dent.  It  follows  that  Az  has  a  normal  distribution  with  mean 
zero,  standard  deviation  of  and  variance  A^ .  T  he  second 
property  implies  that  z  follows  a  M  arkov  process.  Now  con¬ 
sider  the  change  in  the  value  of  z  during  a  relatively  long 
period  of  time  T.  This  can  be  denoted  by  z(T)-z(0).  This 
change  can  be  regarded  as  the  sum  of  the  changes  in  z  in  A/ 

small  time  intervalsof  length  A^,  where  7V  =  — .Thus: 

At 

z(r)-z(0)= 

i=\ 

where  s.  (i  =  1,2,.. .A/)  are  random  drawings  from  a  standard¬ 
ised  normal  distribution.  It  follows  that  z(T)-z(0)  is  normally 
distributed  with  mean  zero,  variance  of  A/  At  =  T  and  stan¬ 
dard  deviation  ^/T.  The  process  described  so  far  has  a  drift 
rate  of  zero  and  a  variance  rate  of  1.  T  he  expected  value  of  z 
at  any  future  time  is  equal  to  its  current  value  and  the  variance 
ofthe  change  in  z  in  a  time  interval  T  equals!. 

A  generdlised  1/1/  iener  process  for  a  variable  x  generalises  the  con¬ 
cept  of  the  drift  rate  and  variance  of  such  a  process,  and  may 
be  defined  in  terms  of  dz  as  dx  =  adt  +  bdz  where  a  and  b  are 
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constants  T  headfterm  impliesthatx  hasan  expected  drift  rate 
of  a  per  unit  time.  Without  the  Mz  term,  the  equation  becomes: 

,  ,  ;  ^  dx 

ax  -  adt ,  I.e.  —  =  a  ■ 
dt 

TheMz  term  may  be  regarded  as  adding  noise  or  variability  to 
the  path  followed  by  x.  For  a  small  time  interval,  A^,  the 
change  in  x.  Ax,  is  given  by  /sx  =  aAt+be4^t  where  £•  is  a 
random  draw  from  a  standardised  normal  distribution.  Ax  has 
a  normal  distribution  with  mean  ofaA^,  standard  deviation  of 
b4^t ,  and  a  variance  of  b^b.t.  Similarly,  the  mean  change  in  x 
for  any  time  interval  T  is  normally  distributed  with  mean 
change  in  x  given  by  al,  standard  deviation  of  change  in  x 
given  by  by/r ,  and  variance  of  change  in  x  as  b^T. 

Even  more  generally,  a  stochastic  process  is  defined  as  an  I  to 

process  \t  dx  =  a{x,t)dt  +  b{x,t)dz  . 

A  basic  simulation  of  stock  price  movement  would  be  to  use  a 
generalised  Wiener  process.  T  his  is  clearly  inadequate  as  this 
assumes  that  both  a  constant  drift  rate  and  constant  variance 
rate  occur,  i.e.  the  percentage  stock  return  is  dependent  upon 
stock  price.  The  constant  expected  drift  rate  assumption  is 
inappropriate  and  is  replaced  by  the  assumption  that  the 
expected  drift,  expressed  as  a  proportion  of  the  stock  price,  is 
constant.  Thus,  if  S  is  the  stock  price,  the  expected  drift  rate  in 
S  is  //A  for  some  constant  parameter  ju  and  for  a  small  time 
interval.  At,  the  expected  change  in  S  is  juSAt.  If  the  variance 
rate  of  the  stock  price  is  always  zero,  then: 

ds 

dS  =  juSdt  or  —  =  Mdt 

O 

i.e.,  S  =  S,e^‘ 
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where  is  the  stock  price  at  time  zero.  This  indicates  that 
when  the  variance  rate  is  zero,  the  stock  price  grows  (or 
declines)  at  a  continuously  compounded  rate  of  //  per  unit 
time.  In  practice,  as  stock  prices  exhibit  volatility,  a  reasonable 
assumption  is  that  the  variance  of  the  percentage  return  in  a 
short  period  of  time.  At,  is  the  same  regardless  of  stock  price. 
Define  as  the  variance  rate  of  the  proportional  change  in 
the  stock  price.  T  hus  a^At  is  the  variance  of  the  proportional 
change  in  stock  price,  S,  during  time  At.  The  instantaneous 
variance  rate  for  S  is  .  T  his  implies  that  S  can  be  repre¬ 
sented  by  an  Ito  process  that  has  an  instantaneous  drift  rate  juS 
and  instantaneous  variance  rate  This  can  be  written: 

dS  =  //Sdt  +  crSdz 
or 


—  =  judt  +  CJdz 
S 


and 


(J 

~1 


2  ^ 


dt  +  adz . 


The  change  in  In  S  between  times  t  and  T  is  thus  normally 
distributed: 


V 


{T  -t),<7-sjT  -t 


where  is  the  stock  price  at  a  future  time  T,  S  is  the  stock 
price  at  the  current  time,  and  ^(m,s)  denotes  a  normal  distri¬ 
bution  with  mean  m  and  standard  deviation  s. 
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It  follows  that: 


ln;S  + 

r 

1^ — ^ 

{T  ~t),(J^T  -t 

l  2  ; 

thus  In  St  is  normally  distributed  so  that  Sj  has  a  log-normal 
distribution. 

U  nder  closer  examination,  stock  prices  in  fact  depart  from  log¬ 
normal  behaviour.  Examination  of  movements  in  stock  prices 
presents  changes  greater  than  our  model  predicts.  Stock 
returns  exhibit  leptokurtosis,  i.e.  the  likelihood  of  returns  near 
the  mean  and  of  large  returns  is  greater  than  our  Brownian 
motion  model  predicts,  whilst  other  returns  tend  to  be  less 
likely  [11,  12].  T able  1.1  gives  a  precis  of  the  statistics  of  this 
non-normal  behaviour. 


D  ecade 

M  ean 
% 

Standard 

Deviation 

% 

Skewness 

Kurtosis 

M  ax 
% 

M  in 
% 

1920s 

0.0058 

2.1453 

0.1321 

11.9272 

11.6396 

-13.7203 

1930s 

-0.0216 

1.8233 

0.2856 

4.5422 

13.8635 

-8.7776 

1940s 

0.0098 

0.7618 

-1.0890 

12.8516 

6.5275 

-7.0431 

1950s 

0.0470 

0.6578 

-0.9201 

7.3343 

4.0476 

-6.7660 

1960s 

0.0063 

0.6558 

0.0437 

5.4873 

4.5787 

-5.8815 

1970s 

0.0025 

0.9262 

0.2693 

1.7889 

4.9517 

-3.5660 

1980s 

0.0489 

1.1607 

-4.3115 

99.3933 

9.6661 

-25.6315 

Overall 

0.0142 

1.1329 

-0.7919 

28.8256 

13.8635 

-25.6315 

Table  1.1:  Decade-by-Decade  Behaviour  of  Daily  Returns 
from  the  D  ow  J  ones  I  ndex 


Clustering  in  Tim e 

W  hat  is  of  particular  interest  is  that  T urner  and  Weigel's  data 
strongly  suggest  the  occurrence  of  temporal  clustering.  Over 
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the  1928-1989  period,  12.5  and  37.5  percent  of  all  extreme 
positive  jumps  in  the  S&P  500  occurred  within  one  and  five 
days  respectively  of  another  positive  jump  in  equity  prices. 
Positive  jumps  in  the  Dow  Jones  were  similarly  clustered  with 
11.3  percent  of  the  positive  jumps  taking  place  within  one  day, 
and  36.2  percent  transpiring  within  5  days  of  each  other. 

The  second  defect  with  the  Brownian  motion  model  is  that  if 
the  model  holds,  then  stock  returns  should  be  proportional  to 
elapsed  time  and  the  standard  deviation  of  returns  should  be 
proportional  to  the  square  root  of  elapsed  time.  T  his  is  based 
upon  the  scaling  properties  of  Brownian  motion.  Turner  and 
Weigel's  data  also  demonstrated  that  monthly  and  quarterly 
volatilities  are  higher  than  annual  volatilities  and  conversely, 
that  daily  volatilities  are  lower  than  annual  volatilities,  i.e. 
their  research  shows  that  the  stock  returns  do  not  scale  in  a 
Brownian  motion  sense. 

T  he  answer  to  these  problems  can  perhaps  be  found  by  con¬ 
sidering  the  mathematical  form  of  these  distributions  [13].  We 
have  shown  that  scaling  is  observed  within  Brownian  motion 
and  that  this  may  be  characterised  by  a  Gaussian  probability 
distribution  plot.  A  log-normal  distribution  may  be  crudely 
assumed  to  be  a  Gaussian  plot  with  an  increased  tail.  The  tails 
of  Pareto  distributions  also  die  off  much  more  slowly  than 
Gaussian  or  log-normal  tails.  Such  fetteil  probability  distribu¬ 
tions  are  thus  increasingly  describing  the  greater  and  greater 
volatility  of  the  system.  Current  financial  models  rely  on  an 
extension  of  the  basic  Brownian  motion  model,  either  assum¬ 
ing  a  different  volatility  distribution  (as  in  "jump  volatility 
calculations")  or  attempting  to  empirically  analyse  and  fit  to 
the  observed  volatility.  For  example,  for  a  Pareto  probability 
distribution  of  the  random  variabley,  with  y>l  then  we  obtain: 
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i.e.  a  power-law.  For  further  discussion  of  these  and  related 
ideas  from  Complexity  Theory  and  dynamical  systems  used  in 
financial  mathematics,  see  references  10  and  13. 

Summary 

In  summary,  we  have  looked  in  some  depth  at  the  complex 
behaviour  of  natural  biological  and  physical  systems.  From 
our  analysis  of  these  open  and  dissipative  systems,  it  is  clear 
that  there  are  a  number  of  key  properties  of  complexity  that 
are  important  to  our  consideration  of  the  nature  of  future  war¬ 
fare.  Such  futures,  involving  the  exploitation  of  loosely 
coupled  command  systems  such  as  Network  Centric  Warfare, 
will  have  to  take  account  of  these  key  properties.  A  list  of  these 
is  given  here,  and  then  discussed  further  in  Chapter  2  in  the 
context  of  N  etwork  C  entric  Warfare. 

1.  Nonlinear  Interaction:  this  can  give  rise  to  surpris¬ 
ing  and  non-intuitive  behaviour,  on  the  basis  of  simple 
local  coevolution. 

2.  Decentralised  Control:  the  natural  systems  we 
have  considered,  such  as  the  coevolution  of  an  ecosys¬ 
tem  or  the  movement  of  a  fluid  front  through  a 
crystalline  structure,  are  not  controlled  centrally.  T  he 
emergent  behaviour  is  generated  through  local 
coevolution. 

3.  Self-Organisation:  we  have  seen  how  such  natural 
systems  can  evolve  over  time  to  an  attractor  correspond¬ 
ing  to  a  special  state  of  the  system,  without  the  need  for 
guidance  from  outside  the  system. 
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4.  Nonequilibrium  Order:  the  order  (for  example,  the 
space  and  time  correlations)  inherent  in  an  open,  dissi¬ 
pative  system  far  from  equilibrium. 

5.  Adaptation:  we  have  seen  how  such  systems  are  con¬ 
stantly  adapting- clusters  or  avalanches  of  local 
interaction  are  constantly  being  created  and  dissolved 
across  the  system.  T  hese  correspond  to  correlation 
effects  in  space  and  time,  rather  a  top-down  imposition 
of  large-scale  coincidences  in  space  and  time. 

6.  Collectivist  Dynamics:  the  ability  of  elementsto 
locally  influence  each  other,  and  for  these  effects  to  rip¬ 
ple  through  the  system,  allows  continual  feedback 
between  the  evolving  states  of  the  elements  of  the 
system. 
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Concepts  for 

VIPRFfiRE  FROM 

Complexity  Theory 


As  a  starting  point  for  our  journey, 
Chapter  1  established  in  some  depth  and 
detail  the  key  ideas  and  methods  from  Complex¬ 
ity  T  heory  that  we  can  bring  to  bear  in  thinking 
about  and  modelling  future  warfare.  I  n  his  RU  SI 
keynote  talk  [1],  ADM  Cebrowski  (H  ead  of  the 
Office  of  Force  Transformation,  U.S.  DoD) 
indicated  that  Network  Centric  Warfare  is  an 
emerging  theory  of  war  based  on  the  concepts  of 
nonlinearity,  complexity,  and  chaos.  It  is  less 
deterministic  and  more  emergent;  it  has  less 
focus  on  the  physical  than  the  behavioural;  and 
it  has  less  focus  on  things  than  on  relationships. 
It  is  clear  from  the  discussion  and  modelling  in 
Chapter  1  that  Complexity  Theory  is  the 
essence  of  these  ideas. 
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I  n  a  previous  book  [2]  we  showed  how  it  is  possible  to  capture 
the  effects  of  command  and  control  in  agent-based  simulations 
of  I  nformation  Age  warfare.  T  his  is  done  by  representing  the 
process  as  the  interaction  of  top-down  and  bottom-up  effects. 
These  are  described  as  Deliberate  Planning  and  Rapid  Plan¬ 
ning.  Deliberate  Planning  is  appropriate  when  ample  time  is 
available  for  the  consideration  of  a  number  of  alternative 
courses  of  action  by  either  "side”  and  a  course  of  action  can  be 
chosen  that  is  considered  to  be,  in  some  sense,  optimal.  Rapid 
Planning  is  appropriate  when  time  is  short  and  expert  deci¬ 
sionmaking  under  stress  leads  to  a  pattern-matching  approach. 
To  quote  from  [2]: 

C ombdt  is,  by  its  neture,  a  compiex  ectivity.  Ashby's  L  aw  of 
R  equisite  Veriety... which  emerged  from  thetheoreticei  consider- 
dtion  of  genera/  systems  as  pert  of  Cybernetics,  indicetes  thet  to 
property  controi  such  a  system,  the  veriety  of  the  controiier  (the 
number  ofeccessibiestetes  which  it  cen  occupy)  must  metch  the 
vdridty  of  the  combdt  system  itself  T  he  control  system  itself,  in 
other  words,  hes  to  be  complex.  Somepr&ious  ettempts  at  repre¬ 
senting  C 2  in  combdt  models  hove  tdken  the  View  thet  this  must 
in&itdbly  ledd  to  extremely  complex  models.  H  owei/er,  recent 
d&dopments  in  C omplexity  T  heory... indicdtednother  way  for¬ 
ward.  T  heessentidi  idea  is  thdtd  number  of  interdcting  units, 
behdving  under  smell  numbers  of  simple  rules  or  digorithms,  cen 
generdte  extremely  complex  beheviour,  corresponding  to  an 
extremely  lerge  number  of  dccessiblestetes,  or  a  high  veriety  con- 
figurdtion,  in  Cybernetic  terms.  It  follows  thdt,  if  we  choose  these 
simple  interdctions  cdrefully,  the  resul  tent  represented  on  ofC2 
will  be  sufficient  to  control,  in  an  eccepteblewey,  the  underlying 
combdt  model.  As  pert  of  this  cereful  choice,  we  need  to  ensure 
thdt  the potentielly  cheotic  beheviour  genereted  by  the  interection 
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of  these  simple  rules  is  'demped'  bye  top-down  C2  structure 
which  remains  focused  on  the  overall,  high  level,  campaign 
objectives.... 

...It  follows,  from  whatwehave  just  said,  that  the  representation 
of  the  C 2  process  must  reflect  two  different  mechanisms.  T  he 
first  is  thelower-l&el  interaction  of  simple  rules  or  algorithms, 
which  generate  the  required  system  variety.  T  he  second  is  the 
need  to  damp  these  by  a  top-down  C 2  process  focused  on  cam¬ 
paign  objectives.  E  ach  of  these  has  to  be  capable  of  being 
represented  using  the  same  6  eneric  HQ/C  ommand  A  gent  object 
architecture  1/1/  e  have  chosen  to  do  this  by  following  the  general 
psychological  structure  of  R  asmussen's  L  adder,  as  a  schema  for 
the  decisionmaking  process.  A  t  the  lower  levels  of  command 
(below  about  Corps,  and  equivalent  in  other  environments),  this 
will  consist  of  a  stimulus/  response  mechanism.  In  cybernetic 
terms,  this  is  feedback  control.  At  the  higher  level,  a  broader 
(cognitive based)  review  of  the  options  available  to  change  the 
current  campaign  plan  (if  necessary)  will  be  carried  out  In 
cybernetic  terms,  this  is  feedforward  control  since  it  involves  the 
use  of  a  'model'  (i.e,  a  model  within  our  model)  to  predict  the 
effects  of  a  particular  system  change 

I  n  the  last  chapter  of  [2]  (C  hapter  6:  "Paths to  the  Future”)  the 
following  point  is  made,  which  is  the  foundation  for  all  of  the 
work  and  ideas  presented  here: 

M  odd  ling  and  analysis  to  ddiermine  the  effect  of  such  phenomena 
underpin  our  thinking  about  such  future  conflict,  the  representa¬ 
tion  of  information  and  command  bang  at  thdr  heart  A  new 
approach  to  capturing  these  effects  has  been  put  forward  in  this 
book,  and  is  having  a  significant  influence  on  the  approach  to 
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modelling  these  phenomend.  How&er,  capturing  the  process  of 
intelligent  agents  in  conflict,  sef  within  a  widely  divergent  sd:  of 
possible  futures,  leads  to  a  rich  set  of  possible  trajectories  of  sys¬ 
tem  evolution  for  analysis  to  consider.  1/1/  e  thus  need  to 
complement  this  effort  with  other  work  to  categorise  and  under¬ 
stand  the  classes  of  behaviours  which  might  emerge  from  such  a 
complex  situation.  T  his  is  the  domain  of  C  omplexity  T  heory. 

The  overall  aim  is  thus  to  develop  an  "Operational  Synthesis” 
(as  discussed  in  [3])  of  both  agent-based  modelling  approaches 
(as  described  in  [2])  and  higher  level  mathematical  metamod¬ 
els  based  on  Complexity  T  heory.  Reference  [2]  lays  out  some 
initial  ideas  on  how  to  develop  such  an  understanding  based 
on  a  theoretical  approach  to  the  development  of  a  higher  level 
"metamodel”  of  a  cellular  automaton  model  of  conflict  such  as 
the  ISAAC  model  developed  by  the  U.S.  M  arine  Corps  Com¬ 
bat  Development  Centre.  We  will  develop  these  ideas  further 
in  Chapters  4  and  5.  In  Chapter  5  in  particular,  we  will  con¬ 
sider  again  the  ISAAC  model  and  demonstrate  how  the  ideas 
of  Complexity  Theory  lead  to  an  understanding  of  the  clusters 
forming  and  reforming  in  such  a  model,  and  how  they  relate 
to  the  emergent  behaviour  of  the  model. 

For  the  moment,  let  us  consider  the  list  of  key  concepts  from 
C  omplexity  T  heory  at  the  end  of  C  hapter  1.  T  hese  are  on  the 
left-hand  side  of  Table  2.1.  On  the  right-hand  side  is  an  inter¬ 
pretation  in  terms  of  military  behaviour  and  doctrine  (which 
we  have  termed  "an  Information  Age  force  structure"). 

T  he  nature  of  N  etwork  Centric  Warfare  for  such  future  I  nfor- 
mation  Age  forces  can  be  outlined  as:  within  a  broad  intent 
and  constraints  available  to  all  the  forces,  the  local  force  units 
self-synchronise  under  mission  command  in  order  to  achieve 
the  overall  intent  [4]. 
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Com PLExiTY  Concept 

iNFORMATION  AGE  FORCE 

Nonllneer  Interaction 

C  ombat  forces  composed  of  a  iarge 
number  of  noniineariy  interacting  parts 

D  ecentrallsed  C  ontrol 

T  here  is  no  master  "oracie"  dictating  the 
actions  of  each  and  every  combatant. 

Self-Organization 

Locai  action,  which  often  appears 
"chaotic,"  induces  iong-range order. 

Nonequilibrium  Order 

M  iiitary  confiicts,  by  their  nature,  pro¬ 
ceed  far  from  equiiibrium.  Correiation 
of  iocai  effects  is  key. 

Adaptation 

Combat  forces  must  continuaiiy  adapt 
and  coevoive  in  a  changing 
environment. 

Collectivist  Dynamics 

T  here  isa  continuai  feedback  between 
the  behaviour  of  combatants  and  the 
command  structure. 

Table  2.1:  Relation  Between  Complexity  and  Information  Age  Warfare 


T  his  process  is  enabied  by  the  abiiity  of  the  forces  invoived  to 
robustiy  network.  We  can  describe  such  a  system  as  loosely  cou¬ 
pled  to  capture  the  iocai  freedom  avaiiabie  to  the  units  to 
prosecute  their  mission  within  an  awareness  of  the  overaii 
intent  and  constraints  imposed  by  high-ievei  command.  This 
aiso  emphasises  the  iooser  correiation  and  nonsynchronous  reia- 
tionship  between  inputs  to  the  system  (e.g.,  sensor  reports)  and 
outputs  from  the  system  (e.g.,  orders),  in  this  process,  informa¬ 
tion  is  transformed  into  "shared  awareness,”  which  is  avaiiabie 
to  aii.  T  his  ieadsto  units  iinking  up  with  other  units,  which  are 
either  iocai  in  a  physicai  sense  or  iocai  through  (for  exampie) 
an  information  grid  or  intranet  (seif-synchronisation).  This  in 
turn  ieadsto  emergent  behaviour  in  the  battiespace,  as  shown 
in  Figure2.1. 
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information 


shared  self-  emergent 

awareness  synchronization  behavdoiir 


Figure2.1:  Information  L  eading  to  E mergent  B ehaviour 

Compare  these  ideas  with  the  broad  conceptual  framework  of 
Complexity  Theory,  as  summarised  below. 

The  Conceptual  Framework  of  Complexity 

Prof.  M  urray  Gell-M  ann  [5]  traces  the  meaning  to  the  root  of 
the  word.  Plexus  means  braided  or  entwined,  from  which  is 
derived  compl&us,  meaning  braided  together,  and  the  English 
word  complex  is  derived  from  the  Latin.  Complexity  is  therefore 
associated  with  the  intricate  intertwining  or  inter-connectivity 
of  elements  within  a  system  and  between  a  system  and  its  envi¬ 
ronment.  In  a  human  system,  connectivity  means  that  a 
decision  or  action  by  any  individual  (group,  organisation,  insti¬ 
tution,  or  human  system)  will  affect  all  other  related 
individuals  and  systems.  T  hat  effect  will  not  have  equal  or  uni¬ 
form  impact,  and  will  vary  with  the  state  of  each  related 
individual  and  system  at  that  time.  T  he  state  of  an  individual 
and  system  will  include  its  history  and  its  constitution,  which  in 
turn  will  include  its  organisation  and  structure.  Connectivity 
applies  to  the  interrelatedness  of  individualswithin  a  system,  as 
well  as  to  the  related  ness  between  human  social  systems,  which 
include  systems  of  artifacts  such  as  information  systems  and 
intellectual  systems  of  ideas. 

T  he  term  complexity  is  used  to  refer  to  the  theories  of  complex¬ 
ity  as  applied  to  complex  adaptive  systems  (CAS).  These  are 
dynamic  systems  able  to  adapt  and  change  within,  or  as  part 
of,  a  changing  environment.  It  is  important  to  note,  however. 
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that  there  is  no  dichotomy  between  a  system  and  its  environ¬ 
ment  in  the  sense  that  a  system  always  adapts  to  a  changing 
environment.  T  he  notion  to  be  explored  is  rather  that  of  a  sys¬ 
tem  closely  linked  with  all  other  related  systems  making  up  an 
"ecosystem.”  Within  such  a  context,  change  needs  to  be  seen 
in  terms  of  coevolution  with  all  other  related  systems  as  we  saw 
in  Chapter  1,  rather  than  as  adaptation  to  a  separate  and  dis¬ 
tinct  environment. 

Collectivist  Dynamicsand  Feedback 

The  phenomenological  definition  of  a  complex  system  isthatit 
exhibits  nonlinear,  emergent,  adaptive  behaviour.  Nonlinear 
behaviour  is  associated  with  far-from-equilibrium,  open  sys¬ 
tems,  in  that  cause  and  effect  are  no  longer  linearly  connected. 
This  is  ultimately  due  to  the  type  of  internal-external  system 
interactions  (feedback)  affecting  our  system. 

Self-Organisation  and  Clustering 

Self-organisation  in  this  context  is  taken  to  mean  the  coming 
together  of  a  group  of  individuals  to  perform  a  particular  task. 
T  hey  are  not  directed  by  anyone  outside  the  group.  T  his  is  not 
the  same  as  "self-management,"  as  no  manager  outside  the 
group  dictates  that  those  individuals  should  belong  to  that 
group,  what  they  should  do,  or  how  it  should  be  done.  1 1  is  the 
group  members  themselves  who  choose  to  come  together,  who 
decide  what  they  will  do  and  how  it  will  be  done.  A  feature  of 
these  groups  is  that  they  are  informal  and  often  temporary. 
Enabling  self-organisation  can  often  be  a  source  of  innovation. 
M  ilitary  commanders  who  understand  the  nature  of  euftreg- 
stdktik^  have  always  understood  this:  a  commander  must 
regard  his  superior's  intention  as  sacrosanct  and  make  its 
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attainment  the  underlying  purpose  of  everything  he  does.  H  e 
is  given  a  task  and  resources  and  any  constraints,  and  within 
thisframework  he  is  left  to  make  his  plan. 

We  already  know  from  C  hapter  1  that  some  complexity  mod¬ 
els  of  natural  ecosystems  use  extremal  dynamics  as  a 
behavioural  driver.  In  systems  where  the  dynamical  evolution 
is  a  struggle  against  various  types  of  thresholds  or  barriers,  the 
action  will  predominately  occur  where  the  net  barrier  to 
change  is  the  smallest.  The  Bak-Sneppen  evolution  model 
described  in  Chapter  1  is  an  example  of  such  an  extremal 
model.  The  species  with  the  lowest  fitness  coevolves  first.  Sim¬ 
ilarly,  in  considering  the  movement  of  a  fluid  through  a 
medium,  the  boundary  moves  where  the  pinning  force  is 
smallest.  T  he  avalanches  of  active  sites  are  the  clusters  created 
in  such  models.  I  n  the  model  ecosystem,  the  system  self-organ- 
ises  towards  a  critical  point  where  it  has  the  greatest 
dynamical  freedom-clustersof  all  sizes  can  potentially  be  cre¬ 
ated,  and  the  statistics  of  the  emergent  distribution  of  such 
cluster  sizes  can  be  predicted-the  distribution  is  of  power  law 
form.  Such  extremal  dynamics  echo  one  of  the  key  tenets  of 
manoeuvre  warfare- namely,  to  focus  your  strength  against 
your  enemy's  weakness. 

Forest  Fires,  Clusters  of  Trees,  and 
Casualties  IN  VdPR 

The  "forest  fire"  model  is  another  example  of  a  system  that 
evolves  to  a  critical  point  through  a  process  of  local  interac¬ 
tion.  In  this  case,  we  start  with  an  empty  two-dimensional  grid. 
At  each  iteration  of  the  process,  with  a  certain  probability  p 


^Auftragstaktik  is  directed  control,  as  opposed  to  befehlstaktik  (detailed 
order  tactics). 
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(normally  close  to  one),  we  drop  a  tree  onto  a  random  grid 
point.  If  the  grid  point  is  blank,  the  tree  is  planted.  If  there  is 
already  a  tree  there,  the  new  tree  is  discarded.  With  probabil¬ 
ity  (I -p)  we  drop  a  spark  onto  a  random  site  instead  of  a  tree.  If 
the  site  is  bare  the  spark  goes  out.  If  there  is  a  tree  on  the  site, 
the  tree  and  all  of  its  immediate  neighbours  burn.  All  of  the 
neighbours  of  these  trees  then  burn  and  so  on  until  the  com¬ 
plete  cluster  of  linked  trees  is  burned  (this  is  termed  a  forest  fire). 
The  rate  (which  is  [I-p]  if  each  iteration  of  the  process  is 
counted  as  a  unit  of  time)  of  sparks  dropping  onto  the  grid  is 
termed  the  sperking  fregumcy.  This  sparking  frequency  sp  is  a  key 
driver  of  the  dynamics  of  the  forest  ecosystem.  If  sp  is  small, 
very  large  clusters  of  trees  are  allowed  to  form,  which  span  the 
entire  grid.  When  a  spark  is  then  dropped,  the  forest  fire  wipes 
out  an  entire  forest  stretching  from  one  side  of  the  grid  to  the 
other.  In  Complexity  Theory,  this  is  known  as  snapping 
noise [6].  This  name  comes  from  looking  at  the  behaviour  of 
the  system  over  time-large  spikes  of  tree  extinction  (forest  fires) 
are  created  at  isolated  points  in  time.  If  the  sparking  frequency 
sp  is  very  large,  then  tree  clusters  do  not  have  the  chance  to 
grow.  T  hus,  over  time,  the  system  produces  a  large  number  of 
small  spikes  of  activity,  which  are  called  popping  noise.  When  sp 
is  in  the  intermediate  regime,  the  system  self-organises  to  a 
critical  state  where  the  clusters  of  burnt  trees  have  a  distribu¬ 
tion  represented  by  a  power  law,  and  clusters  of  all  sizes  can  be 
created.  Over  time,  the  spikes  produced  by  this  process 
(i.e.,  the  time  evolution  of  forest  fires  of  various  sizes)  have  a 
similar  dynamic  to  that  produced  by  the  acoustic  dynamics  of 
crumpling  paper  [6],  and  so  this  regime  istermed  creckiing  noise. 

It  is  possible  [7]  to  relate  such  self-organised  behaviour  of  a 
forest  fire  model  to  the  statistics  of  the  scale  and  intensity  of 
conflicts.  This  is  the  beginning  of  an  explanation  as  to  why 
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casualties  in  war  follow  a  power  law  distribution.  As  noted 
in  [7],  "...the  behaviour  of  the  forest  fire  model  can  be 
explained  in  terms  of  a  cascade  model.  If  trees  are  randomly 
planted  on  a  grid,  the  distribution  of  cluster  sizes  is  exponential 
(Poissonian),  not  power  law  (fractal).  The  distribution  of  cluster 
sizes  in  the  forest  fire  model  is  power  law  (fractal).  This  is 
because  clusters  of  trees  continuously  grow  and  combine  to 
form  larger  clusters.  Small  fires  sample  this  population  of  clus¬ 
ters,  but  the  loss  of  trees  in  fires  is  dominated  by  the  largest 
fires.  There  is  a  self-similar  cascade  of  trees  from  small  to  large 
clusters.  I  n  terms  of  the  forest  fire  model,  a  spark  ignites  a  tree 
and  the  model  fire  consumes  the  entire  cluster  to  which  this 
tree  belongs.  T  his  is  also  the  case  for  real  forest  fires.  Ignition 
of  the  forest  must  take  place  for  a  fire  to  take  place,  and  the  fire 
will  then  spread  through  the  contiguous  flammable  material. 
A  war  must  begin  in  a  manner  similar  to  the  ignition  of  a  for¬ 
est.  0  ne  country  may  invade  another  country,  or  a  prominent 
politician  may  be  assassinated.  The  war  will  then  spread  over 
the  contiguous  region  of  metastable  countries...."  We  will  look 
at  this  fractal  and  power  law  behaviour  of  casualties  in  more 
detail  in  Chapter  3. 

Tuning  and  Goal  Seeking 

T  his  leads  us  to  consider  the  concept  of  tuning.  Self-organising 
systems  (such  as  the  Bak-Sneppen  evolution  model  of  an  eco¬ 
system  described  in  Chapter  1)  can,  as  their  name  implies, 
develop  local  organisation  within  the  system  in  order  to  evolve 
towards  an  attractor.  Tuning  can  be  seen  as  a  directive  way  for 
the  macrosystem  to  attempt  to  influence  the  behaviour  of  the 
microsystem.  A  controlling  intelligence  is  deemed  to  be  neces¬ 
sary  in  order  to  guide  the  system  towards  a  particular  goal.  An 
example  of  this  is  the  tuning  of  the  sparking  frequency  model 
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parameter  sp  in  the  forest  fire  model  in  order  to  obtain  critical¬ 
ity  (and  hence  power  law  statistics  for  the  size  of  clusters  of 
burnt  trees).  Varying  the  tuning  parameter  (the  sparking  fre¬ 
quency)  of  the  forest  fire  model  represents  intervention  from 
outside  the  system  in  order  to  ensure  that  it  heads  towards  a 
particular  goal.  T  his  question  of  tuning  makes  us  consider  the 
bounddries  of  the  systems  we  are  examining,  and  the  flux  of 
energy  and/  or  information  across  the  system  boundary. 

Information  Flux  across  the  Boundary  of  an 
Open  System 

I  n  such  open  dissipative  systems,  there  will  always  be  fluxes  of 
information  and/  or  energy  across  the  system  boundary.  As  an 
example  of  this,  we  show  in  Chapter  4  how  it  is  possible  to 
quantify  the  flux  of  information  across  the  boundary  of  an 
open  system  (a  wargame)  using  the  concept  of  Information 
Entropy.  We  can  then  use  this  to  relate  such  a  flux  of  informa¬ 
tion  to  emergent  properties  of  the  wargame  (such  as  the 
number  of  casualties  suffered).  T  he  idea  of  using  entropy  as  a 
measure  of  information  (and  hence  knowledge)  is  then  applied 
in  Chapter  4  to  show  how  the  benefits  of  network-centric  (as 
opposed  to  platform-centric)  approaches  to  specific  task  prose¬ 
cution  can  be  quantified  (this  latter  based  on  work  by  the 
RAND  Corporation). 

REFERENCES 

1  CEBROWSKI  A  (2000).  "Network  Centric  Warfare  and  Information 
Superiority."  Key  note  ad  dress  from  proceedings,  Royal  U  nited  Services 
I  nstitute  (RU  SI )  conference.  "C  41 STA R ;  Achieving  I  nformation 
Superiority."  RU  SI .  W  hitehall,  London,  U  K . 

2  M  0  F  FAT  J  (2002).  C  ommand  and  C  ontrot  in  the  I  nformation  h  ge:  R  epresenting  its 
Impact.  T  he  stationery  Office.  London,  U  K  . 


R  eferences 


56 


C  omplexity  T  heory  and  N  etwork  C  entric  W  arfare 


3  HOFFMAN  F  G  and  HORNE  G  E  (1998).  M  anmerWarfareSciencel998. 
Deptof  the  Navy,  HO  U.S.  M  arine  Corps.  Washington,  DC . 

4  ALBERTSD  S,GARSTKAJJ  and  STEIN  F  P  (1999).  N^ork  Centric 
W arfare:  D evdoping and  L evaaging I nformation  Superiority.  CCRP,  DoD. 
Washington,  DC,  USA. 

5  GELL-MANN  M  (1994).  T  he  Quark  and  the  J  aguar.W  H  Freeman.  New  York, 
USA. 

6  SETHNAJ  RDAHMEN  K  A,andMYERSC  R  (2001).  "Crackiing 
N Oise."  Nature.  410  8.  M  arch  2001.  pp.  242-250. 

7  RO  BERTS  D  C  andTURCOTTE  D  L  (1998).  "Fractaiity  and  Seif- 
Organised  Criticaiity  of  Wars."  Fractais.  6  No  4.  pp.  351-357. 


R  eferences 


Chapter  3 


Evidence  for 
Complex  Emergent 
Behaviour  in 
Historical  Data 


Introduction 

The  first  point  to  make  is  that  the  exploita¬ 
tion  of  manoeuvre  warfare  is  not  new. 
Commanders  have  exploited  such  an  approach 
in  previous  generations.  For  example,  LTG  Sir 
Francis  Tuker  [1]  indicated  that  at  a  three- 
dimensional  spatial  level,  manoeuvre  warfare  is 
determined  by  three  conditions: 

1.  Flanksshall  be  tactically  open  or  it  shall  be 
possible  to  create  a  flank  by  break-in  and 
breakthrough. 

2.  T  he  mobile  arm  shall  be  predominant. 
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3.  It  shall  be  possible  to  administer  the  mobile  arm  to  the 
point  where  it  will  decide  the  battle  and  gain  decisive 
victory. 

Examination  of  historical  data  should  thusgive  us  some  insight 
into  these  effects. 

Secondly  the  idea  that  execution  of  military  network-centric 
enterprises  within  a  battlespace  is  of  a  self-organising  nature 
has  interesting  consequences  for  the  future  of  both  warfighting 
and  other  operations.  A  system  that  comprises  of  clusters 
grouping  and  regrouping,  the  increasing  tempos  of  battlespace 
awareness,  operations,  and  responsiveness  etc.  require  that  the 
command  decision  process  (or  better,  the  rate  of  change  of 
decision  process)  is  undertaken  with  a  similar  tempo.  T  he  ratio 
of  "battle  speed”  to  "C2  speed'”  is  thus  a  critical  issue,  as  we 
discussed  in  [2].  W  hen  this  ratio  is  high,  the  system  reverts  to 
being  self-organising  in  nature.  In  considering  such  self-organ¬ 
ising  systems  earlier,  we  indicated  that  under  certain 
conditions,  the  time  series  of  outputs  from  such  a  system 
should  have  a  certain  fractal  character.  In  the  case  of  a  model 
ecosystem  such  as  the  Bak-Sneppen  evolution  model  discussed 
in  Chapter  1,  in  fact,  this  fractal  process  is  driven  by  the  extre¬ 
mal  (smallest)  fitness  values  on  the  lattice  of  interacting  species. 

Time  Series  Behaviour 

Applications  of  these  ideas  have  been  made  in  order  to  predict 
combat  casualty  rate  patterns  from  WWII  data  [3].  A  time 
series  model  based  on  thisfractal  self-organising  approach  has 
been  contrasted  with  three  other  prediction  methodologies:  a 
neural  network;  use  of  nonlinear  prediction  (a  prediction  is 
made  by  searching  the  data  for  the  nearest  N  points  in  a  d- 
dimensional  embedding  and  estimating  the  behaviour  charac- 
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teristics  of  the  data  for  the  next  timestep  [4]);  and  lastly,  use  of 
a  maximum  entropy  method  to  calculate  a  power  spectrum 
from  which  a  linear  prediction  may  be  made  [5],  All  of  these 
approaches  are  available  in  a  package  of  time  series  analysis 
procedures  (the  Chaos  Data  Analyser)  produced  by  the  Ameri¬ 
can  Institute  of  Physics  for  the  analysis  of  experimental  data  in 
natural  systems,  and  that  is  what  we  have  used  here.  A  range 
of  the  data  points  from  the  original  data  were  deleted  (those  at 
the  end  of  the  time  series)  and  a  prediction  made  of  these  data 
points,  which  is  then  compared  with  the  original  data  for  the 
2nd  Armoured  Division.  The  plots  (Figure  3.1  to  3.5)  are  of 
casualties  per  1000  on  they-axisand  time  in  days  on  thex-axis. 


Figures. 1: 2 nd  Armoured  Division  Ddtd 


Figures. 2: 2nd  Armoured  Di vision- Power  Spectrum  Prediction 
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Figures. 3: 2ncl  Armoured  Division- Nonlin^r  Prediction 


Figures. 4:  2nd  Armoured  D  ivision-  N  eurai  N  et  Prediction 


Figures. 5: 2nd  Armoured  Division- SO  C  Prediction 

The  first  part  of  this  time  series  (up  to  day  38)  was  in  fact  used 
to  train  a  number  of  different  time  series  prediction  methods, 
and  these  have  been  compared  with  the  predictions  for  the 
days 39  onward,  in  fitting  a  prediction  based  on  a  seif-organis- 
ing  criticaiity  (SOC)  fractai  series,  we  have  assumed  that  the 
circumstances  remain  sufficientiy  constant  that  we  can  fit  a 
singie  SO  C  process  (this  corresponds  to  a  power  spectrum  that 
isiinear  when  piotted  on  a  iog-iog  scaie).  Comparing  the  "jerk- 
iness”  of  the  SOC  prediction  and  the  reai  data,  the  generai 
pattern  of  the  process  is  very  simiiar. 
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We  repeated  the  process  with  another  data  set  drawn  from 
Kuhn  for  the  9th  Armoured  Division,  with  similar  results  as 
summarised  below. 
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Figure3.9:  9th  Armoured  Division- 50 C  Prediction 


Figures. 8:  9th  Armoured  Division- Power  Spectrum  Prediction 


Figures.?:  9th  Armoured  Division- Neurdi  NetPrediction 


Figures. 6:  9th  Armoured  Division  Ddtd 
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These  plots  indicate  that  the  assumption  of  self-organisation 
appears  to  give  casualty  behaviour  that  is  of  the  same  form  as 
for  the  real  data,  at  least  for  these  data  points. 

Kuhn  himself  suggested  that  the  casualty  rate  data  pattern, 
taken  from  his  collection  of  combat  rates  from  World  War  II 
to  the  present,  displayed  a  move  away  from  linearity.  First,  the 
number  of  "hotspots”  in  combat  did  not  increase  with  larger 
force  concentrations,  and  second,  combat  was  characterised 
by  high  rates  of  casualties  lasting  for  short  periods  of  time, 
interspersed  with  low  casualty  rates.  T  hese  results  indicated,  in 
turn,  a  move  away  from  the  type  of  modelling  associated  with 
attrition  warfare.  Q  uantitative  patterns  were  also  found  to  be 
based  upon  two  types  of  operational  forms:  that  of  a  continudi 
front  (characterised  by  combat  on  the  Western  Front  in  World 
War  II)  and  that  of  a  disrupted  front,  (characterised  by  combat 
on  the  Eastern  Front  in  World  War  II).  Combat  casualties  dra¬ 
matically  increased  as  a  consequence  of  breakthrough 
operations  in  disruptive  types  of  operations.  Combat,  there¬ 
fore,  could  be  seen  as  a  process  wherein  quiet  states  are 
interspersed  with  "critical"  irruptions. 

Lauren  [6]  has  also  considered  casualty  data  from  World 
War  II  and  compared  it  with  outputs  from  the  MAN  A  model 
(an  agent-based  simulation  similar  in  character  to  the  ISAAC 
model  we  will  discuss  in  detail  in  C  hapters  4  and  5).  FI  e  has 
shown  that  such  data  display  fractal  properties  and  power 
spectra  that  confirm  the  analysis  of  casualty  data  discussed  in 
this  chapter. 


T ime Series  Behaviour 


Chapters  63 


Further  Historical  Data  on  the 
Processes  of  "I  rruption"  and 
Breakthrough 

In  an  historical  analysis  study  [7]  of  the  operational  level  of 
combat,  it  was  found  by  Rowland  that  the  occurrence  of 
breakthrough,  defined  as  the  destruction  of  cohesiveness  of  the 
defence,  was  an  important  event  in  the  eventual  success  of  an 
offence.  Following  breakthrough,  86%  of  operations  were  suc¬ 
cessful,  whereas  if  no  breakthrough  was  achieved  only  15% 
were  eventually  successful.  Once  breakthrough  has  been 
achieved,  it  becomes  possible  for  the  attack  force  to  conduct  a 
type  of  operation  more  in  the  nature  of  exploitation  than  com¬ 
bat.  M  oreover,  variations  in  the  time  to  breakthrough  also  led 
to  differences  in  the  nature  of  campaigns;  the  timings  of 
"immediate”  (less  than  half  a  day),  "quick"  (less  than  2  days), 
and  "prolonged"  (over  2  days)  were  used  for  study  purposes. 

I  n  studying  the  course  of  operations,  one  of  the  aspects  exam¬ 
ined  was  the  nature  of  the  movements  of  attack  forces.  The 
study  of  the  rate  of  terrain  capture  with  operational  success 
showed  patterns  that  could  be  related  to  breakthrough  time. 
Several  different  measures  were  examined,  the  first  one 
derived  having  the  dimensions  of  advance  rate  (distance/ 
time).  Whilst  the  measures  change  through  time,  the  following 
results  relate  to  measurement  (such  as  the  length  of  attack 
frontage)  at  the  time  of  breakthrough. 

Area  Taken  AND  M  ean  Advance  at 
Breakthrough 

This  measure  may  betaken  as  a  possible  index  of  how  badlythe 
defence  has  fared  and  of  its  chances  of  recovery.  H  owever,  as 
the  simple  area  measure  does  not  allow  for  variations  in  size  of 
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operation,  the  function  (area  at  breakthrough)/  (attack  frontage) 
has  been  used  here.  T  his  has  a  dimension  of  length  and  would 
represent  "mean  advance”  on  the  original  attack  frontage. 

T  he  distributions  of  the  measures  in  T ables  3.1,  3.2,  and  3.3 
are  log-normal  within  a  given  breakthrough  category,  as  we 
will  show  explicitly  later.  In  each  case,  the  table  shows  the  geo¬ 
metric  mean  of  the  corresponding  log-normal  distribution. 
The  discussion  of  stock  price  returns  and  Brownian  motion  in 
Chapter  1  gives  us  some  insight  into  why  log-normal  distribu¬ 
tions  might  arise  in  such  cases. 


Operations  Leading  To 

Cam paign  Success 

Campaign  Failure 

Prolonged 

Breakthrough 

29 

8 

Quick 

Breakthrough 

5 

- 

Immediate 

Breakthrough 

11 

2.2 

Table  3.1:  Geometric  M  ean  Area/ Attack  Front  at  Breakthrough  (M  lies) 


0 PERATioNS  Leading  To 

Campaign  Success 

Campaign  Failure 

Prolonged 

Breakthrough 

4.5 

1.5 

Quick 

Breakthrough 

4.5 

1.5 

Immediate 

Breakthrough 

11 

2.2 

Table  3.2:  Geomei 

trie  M  ean  Area  Per  Day/  At 

tack  Front  at  Breakthrough 

(M  lies/  Day) 
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Operations  Leading  Tc 

Cam paign  Success 

1 : 

Campaign  Failure 

Prolonged 

Breakthrough 

4.5 

2.1 

Quick 

Breakthrough 

8.8 

- 

Immediate 

Breakthrough 

18 

5.5 

Table  3.3:  Geometric  M  ean  VArea  Per  Day  at  Breakthrough  (M  lies/  Day) 


A  first  analysis  of  the  results  in  Table  3.1  shows  that  these 
"mean  advances”  at  breakthrough  are  greater  for  those  achiev¬ 
ing  subsequent  success  than  those  failing  by  factors  of  5  for 
immediate  breakthrough  and  approximately  3.6  for  prolonged 
breakthrough.  This  larger  factor  for  immediate  breakthrough 
is  indicative  of  the  extra  "brittleness"  that  pertains  to  these  very 
quick  breakthrough  cases.  They  also  show  that  these  "mean 
advances"  at  breakthrough  are  less  for  immediate  than  for  pro¬ 
longed  breakthrough  by  mean  factors  of  0.38  for  those 
achieving  subsequent  success  and  by  0.28  for  those  failing. 

Moving  to  the  more  complicated  measure  of  irruption  in 
T able  3.2,  we  move  from  mean  advance  to  mean  advance  per 
day,  giving  a  first  order  measure  of  the  mean  rate  of  advance. 
For  this  case,  there  is  somewhat  less  variability  in  results. 
There  is,  in  fact,  no  significant  difference  between  quick  and 
prolonged  breakthrough  effects  in  this  case.  T  hese  two  groups 
can  therefore  be  pooled.  Combined  quick  and  prolonged 
breakthrough  results  in  terms  of  subsequent  success  and  failure 
are  thus  shown  in  Table  3.2. 

We  again  find  differences  on  this  measure  (mean  rate  of  advance) 
between  breakthroughs  leading  to  eventual  success  and  those 
leading  to  failure.  T  hese  differences  are  by  factors  of  5  for  imme- 
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diate  breakthrough,  and  3  for  prolonged/ quick  breakthrough. 
The  mean  rate  of  advance  for  immediate  breakthrough  is  greater 
than  for  prolonged/  quick  breakthrough. 

The  caricature  of  movement  to  breakthrough  implied  is: 


H  owever,  the  movement  post-bredkthrough  is  better  represented 
by  the  linear  model  than  the  radial  one  for  success  after  break¬ 
through  and  without  breakthrough.  The  exception  to  this  is 
success  after  immediate  breakthrough,  which  discriminates 
success  most  significantly  on  the  basis  of  the  radial  model. 

An  alternative  measure  of  the  same  dimensions  is  thus  Area/ 
Time,  again  a  representation  of  advance  rates  but  representing 
propagation  from  a  point  rather  than  a  linear  advance. 
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The  equivalent  mean  values  for  this  measure  are  shown  in 
T able  3.3.  T  he  differences  between  groups  show  similar  pat¬ 
terns  to  the  previous  measures,  but  are  greater,  and 
significant  at  the  1%  level  (f  test)  between  immediate  and 
prolonged  breakthrough. 

It  can  thus  be  observed  that,  despite  the  great  variations  in  size 
of  the  operations  studied,  there  are  patterns  to  be  deduced  and 
these  can  offer  lessons  on  the  nature  of  this  little-studied  aspect 
of  manoeuvre  warfare. 

This  process  of  irruption  has  been  identified  as  one  of  the  key 
emergent  effects  of  manoeuvre  warfare  [8].  We  consider  now 
whether  such  a  process  has  scaling  properties  of  the  type  dis¬ 
cussed  in  our  general  consideration  of  complexity.  The 
historical  data  indicates  (as  we  have  discussed)  that  for  a 
given  type  of  breakthrough  (immediate,  quick,  or  prolonged- 
I,  Q  ,  or  P),  and  subsequent  effect  on  the  campaign  (Subse¬ 
quent  Success  [SS]  or  Subsequent  Failure[SF]),  the  mean 
advance  at  breakthrough  turns  out  to  be  a  log-normal  distri¬ 
bution.  Of  even  more  interest  to  us  is  the  fact  that  if  these 
distributions  are  plotted  for  each  of  the  breakthrough/  cam¬ 
paign  effect  categories,  then  they  have  a  certain  scaling 
character,  which  we  now  define. 

Consider  Figure  3.10  (page  71).  FI  ere  we  have  plotted  the  log¬ 
normal  distribution  of  mean  advance  at  breakthrough  for  each 
of  a  number  of  breakthrough  categories.  T  he  x-axis  is  the  log 
of  the  mean  advance  at  breakthrough  (in  miles),  and  they-axis 
is  such  that  a  cumulative  normal  distribution  of  the  x  variable 
will  give  a  straight  line.  We  can  see  that  the  various  categories 
of  breakthrough  produce  cumulative  curves  that  are  parallel  to 
each  other,  with  one  exception.  FI  owever,  this  case  has  signifi¬ 
cantly  less  data  than  the  other  cases  and  we  take  this  to  be  the 
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cause  of  this  deviation.  T  he  fact  that  these  iines  iie  paraiiei  to 
each  other  means  the  foiiowing:  given  two  such  curves,  corre¬ 
sponding  to  X  variabiesx(l)  and  x(2),  there  is  a  scaiing  variabie 
X  (which  depends  on  the  two  categories  of  breakthrough  being 
considered)  such  that 

Log  x(l)  =  X  Log  x(2) 
i.e.  x(i)  =  x(2y 

and  the  distribution  of  mean  advance  at  breakthrough  coin¬ 
cides  for  the  variabies  x(l)  and  x{2f-.  i  n  this  sense,  we  can  say 
that  x(2)  can  be  scaied  by  a  power  transformation  so  that  its 
distribution  coiiapsesonto  that  ofx(l). 

An  aiternative  radiai  measure  of  irruption,  as  we  have  dis¬ 
cussed,  is  V(area  at  breakthrough)/ (days  to  breakthrough), 
with  dimensions  of  mi ies  per  day.  if  thisispiotted  on  the  same 
basis  as  the  previous  figure,  we  again  have  evidence  for  a  form 
of  "scaiing  coiiapse"  of  the  type  discussed  above  (Figure  3.11). 
M  oreover,  the  stabiiity  of  the  two  sets  of  data  indicate  that 
there  are  (at  ieast)  two  categories  of  emergent  behaviour  for 
irruption  and  subsequent  campaign  outcome- i inear  advance 
and  radiai  propagation  from  a  point. 

Each  data  point  in  Figures  3.10  and  3.11  is  a  campaign  out¬ 
come,  ciassified  in  terms  of  immediate  (i),  quick  (Q)  or 
proionged  (P)  irruption,  and  Subsequent  Success  (SS)  or  Sub¬ 
sequent  Faiiure  (SF).  T  he  key  to  the  data  points  is  given  beiow: 

■  I(SF)  •  I(SS)  □  Q(SS) 

A  P(SF)  X  P(SS) 
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The  Fractal  Front  of  Combat 

In  [9],  Lauren  discusses  the  fractal  nature  of  a  combat  front 
between  two  opponents  T  he  idea  is  that  an  essentially  straight 
line  frontage  between  two  tactical-level  opponents  will  buckle 
into  a  fractal  shape,  whose  fractal  dimension  can  be  calculated 
asafunction  of  the  force  ratio  of  the  forces  involved  (thenumber 
of  attackers  to  the  number  of  defenders)  as  derived  from  H  istori- 
cal  Analysisof  infantry  battles  carried  out  by  theU  K  Dstl. 

Lauren  [9]  uses  the  H  istorical  Analysis  result  as  a  basis  for  his 
analysis: 

F  =  (number  of  attack  infantry/  number  of  defence  infantry)o685 

where  Fisa  multiplier  for  the  base  number  of  casualties  of  the 
attacking  force  per  defence  weapon.  As  a  consequence  of  this, 
Lauren  was  able  to  show  that  the  combat  front  will  buckle 
over  time  and  in  the  limit  will  have  a  fractal  dimension 
D  =  1.685. 

From  C  hapter  1,  if  we  assume  that  this  process  is  akin  to  inva¬ 
sion  percolation  of  one  fluid  by  another  in  a  porous  medium, 
the  fractal  dimension  of  the  boundary  of  the  resulting  interface 
should  lie  in  the  range  1.33-1.89,  which  is  what  we  find  from 
historical  data. 

It  is  possible  in  this  case  to  derive  the  underlying  dynamics 
producing  this  statistical  effect.  It  turnsout  that  this  fractal  fac¬ 
tor  is  due  fundamentally  to  detection  of  targets  [10],  and 
comes  from  a  model  of  the  engagement  process  that  leads  to 
the  following  relationship: 


1 

R 


h 

T 
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where  are  constants,  R  is  the  defender  rate  of  fire, 

and  T  isthe  number  of  targets  in  view  [11], 

It  reflects  the  asymmetry  of  the  infantry  battle  in  the  following 
sense  [12],  T  he  attack  force  aim  isto  close  on  the  defence  posi¬ 
tion,  and  fire  is  used  in  a  general  suppressive  mode-actual 
casualties  caused  to  the  defence  are  only  a  small  part  of  the 
process  at  this  point.  H  owever,  from  the  defence  perspective, 
the  aim  isto  deter  the  attack,  and  casualties  to  the  attack  force 
are  very  important.  Such  casualties  to  the  attack  force  are  a 
direct  reflection  of  the  intervisibility  of  targets  to  the  defence 
force  as  discussed  above. 

As  with  most  applications  of  fractal  processes,  the  process 
breaks  down  at  some  point  due  to  the  granularity  of  the  resolu¬ 
tion.  In  this  case,  the  process  remains  valid  up  to  about  30m 
closing  distance  between  the  attack  and  defence.  At  that  point, 
a  different  mechanism  comes  into  play,  leading  to  local  defence 
surrender  and  attack  overrun  of  defence  positions  [12]. 

M  ore  generally,  the  figure  of  0.685  relates  to  open  terrain.  I  n 
urban  areas  it  is  about  0.5  [13].  The  closing  to  overrun 
appears  to  occur  differently  in  urban  and  wooded  terrain  as 
compared  with  open  terrain  [12].  For  example,  in  open  condi¬ 
tions,  the  closing  part  of  the  battle  occurs  across  the  front.  By 
contrast,  in  urban  conditions,  the  attack  force  is  split  into  small 
subunits  who  individually  close  on  defence  locations  leading  to 
local  surrender  and  overrun. 
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Power  Law  Relationships  in  Combat  Data 

The  Hartley  Model 

In  [14],  H  artley  has  analysed  eight  separate  databases  of  his¬ 
torical  combat  data.  Four  were  developed  by  Helmbold;  the 
fifth  dataset,  "Inchon,"  was  developed  by  Busse.  T  he  last  three 
datasets  were  developed  by  Dupuy.  T  hese  eight  datasets  span 
several  centuries  in  time,  include  both  air  and  land  conflicts, 
and  span  the  range  from  small  to  large  interactions.  Details  of 
the  databases  are  given  in  [14].  On  the  basis  of  this  extremely 
extensive  set  of  data,  H  artley  was  able  to  develop  a  stable  anal¬ 
ysis  of  the  relationship  between  casualties  in  conflict  and  the 
initial  force  ratio  based  on  earlier  ideas  of  Helmbold.  He 
defines  the  following  two  dimensionless  variables: 

HELMRAT  =  ^ ^ 

yo-y 

FOR  RAT  =  — 

To 

where  in  each  case,  Xq  isthe  starting  value  of  force  size,  and  x  is 
the  final  value  (similarly  for  y).  H  artley  has  established  a  power 
law  relationship  between  these  two  variables,  HELM  RAT  and 
FORRAT,  on  the  basis  of  the  comprehensive  data  sets 
described  above.  H  e  has  shown  that  (in  logarithmic  terms): 

Ln  (HELMRAT)  =  aLn  (FORRAT)+(3 

where  the  expected  value  of  a  is  approximately  1.35  and  the 
value  of  (3  is  approximately  normally  distributed  about  the 
value  -0.22  with  standard  deviation  of  0.7.  H  artley  shows  that 
the  value  of  a  has  the  characteristics  of  a  universal  constant, 
being  stable  over  four  centuries  of  time  [14,  Figure  17],  and 
stable  when  considering  conflicts  of  different  sizes,  ranging 
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from  force  sizes  of  less  that  5,000  to  more  than  100,000  [14, 
Tables4,  5,  and  6], 

If  it  is  assumed  that  the  mechanism  that  produces  this  remark¬ 
ably  stable  relationship  between  casualty  effects  and  force  ratio 
isof  Lanchester  type,  then  H  artley  shows  that  it  must  be  of  lin¬ 
ear-logarithmic  form.  However,  the  relationship  is  based  on 
the  empirical  data  alone,  and  other  explanations  are  possible. 
For  example,  in  [15]  an  analysis  based  on  self-organisation  (in 
particular  the  forest  fire  model  of  C  hapter  2)  is  put  forward  as 
the  basis  for  the  equally  remarkable  scaling  of  conflict  size.  It  is 
thus  persuasive  that  such  complexity-based  effects,  rather  than 
a  Lanchester  process,  lie  at  the  base  of  the  scaling  relationship 
established  by  H  artley.  T  his  analysis  by  T urcotte  and  Roberts 
is  next  discussed. 

In  their  paper  [15],  they  begin  by  comparing  the  predictions 
of  the  theoretical  self-organising  forest  fire  model  with  the  sta¬ 
tistics  of  the  relative  sizes  of  real  forest  fires.  Four  data  sets  are 
considered:  4,284  forest  fires  in  the  USA  Fish  and  Wildlife 
Service  Lands  during  the  period  1986-1995;  120  of  the  largest 
fire  areas  in  the  western  USA  from  tree  ring  data,  spanning 
the  period  1155-1960  (800  years);  164  fires  in  the  Alaskan 
boreal  forests  during  1990-1991,  and  298  fires  in  the  Austra¬ 
lian  Capital  Territory  during  1926-1991.  The  results  are  in 
good  agreement  with  a  power  law  statistical  distribution  of  size 
of  fire  versus  frequency,  with  a  power  law  exponent  of 
between  1.3  and  1.5.  T  he  remarkable  thing  is  the  stability  of 
the  trend  across  such  a  long  period  of  time,  during  which  tech¬ 
nology  has  changed,  as  have  ways  of  fighting  such  fires.  The 
authors  then  show  that  a  similar  power  law  relationship  (also 
with  an  exponent  in  the  same  range)  holds  for  the  intensity  of 
conflict  versus  its  frequency.  This  work  extends  the  research  of 
Richardson  [16],  who  also  showed  a  power  law  relationship 
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between  the  intensity  of  war  and  its  frequency.  T urcotte  and 
Roberts  base  their  results  on  two  data  sets:  that  of  Levy,  which 
tabulates  the  intensities  of  119  wars  from  1495  to  1973;  and 
that  of  Small  and  Singer,  who  considered  118  wars  over  the 
period  1816-1980. 

The  similarity  of  the  power  law  exponent  for  both  forest  fire 
statistics  and  war  intensity  leads  Turcotte  and  Roberts  to 
hypothesise  that  war  deaths  are  caused  by  a  self-organising 
mechanism  akin  to  that  of  the  forest  fire  model.  This  is  at  least 
the  beginning  of  an  explanation  of  why  casualties  in  war 
should  give  rise  to  such  a  simple  power  law  relationship,  stable 
over  centuries  of  time. 
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Mathematical 
Modelling  of 
Complexity, 
Knowledge, 
AND  Conflict 


Introduction 

Understanding  the  behaviour  of  agent-based 
simulation  models  of  conflict  is  now 
becoming  more  important,  especially  as  (with 
improved  representation  of  Command  and  Con¬ 
trol  [1])  the  agents  gain  intelligence  and  try  to 
outsmart  each  other,  producing  potentially  very 
complex  behaviour.  In  the  modelling  of  natural 
systems  (such  as  fluid  dynamics*  or  heat  flow),  the 
principal  variables  in  these  models  can  often  be 
separated  out  from  the  rest  of  the  model  to  pro¬ 
duce  a  mathematical  metamodel  that  is  aimed  at 
relating  the  outputs  of  the  model  to  these  driving 
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inputs  in  a  more  transparent  and  explicit  way.  If  this  can  be 
achieved,  it  improves  our  understanding  of  the  system  and  its 
likely  emergent  behaviour,  as  well  as  complementing  the  use  of 
detailed  simulation.  Such  an  approach  is  consistent  with  the 
idea  of  'Operational  Synthesis'  as  espoused  by  Dr.  Alfred 
Brandstein,  then  Chief  Scientist,  U.S.  M  arineCorps[2]. 

As  an  example,  in  developing  a  metamodel  we  consider  the 
relationship  between  a  key  outcome  of  the  model,  a,  and  a  set 
of  input  variables  as  follows: 

a  =  f{a„...a^,b„b^) 

(T  his  is  easily  generalised  to  an  arbitrary  number  of  bs.)  T  he 
arguments  ay,...a^  have  independent  dimensions.  That  is,  the 
dimension  of  any  a  cannot  be  expressed  as  a  combination  of 
the  dimensions  of  the  other  as.  In  contrast,  the  dimension  of 
each  b  variable  can  be  expressed  as  such  a  combination.  The 
arguments  can  be  transformed  using  a  gauge  transformation 
so  that: 


a[  =  A^a^,...,a[  =Ai^a^- 

T  hese  correspond  to  a  change  in  the  "gauge"  (e.g.,  from  centi¬ 
metres  to  metres  or  kilometres)  in  the  measurement  of  a 
variable.  Physically,  if  a  gauge  change  makes  no  difference  to 
observed  behaviour  for  all  observers  using  different  gauges, 
the  variable  is  said  to  be  sdf-simildr. 
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The  metamodel  function  can  then  be  shown  [1]  to  have  the 
property: 


a  =  /(«! , ,  7^2 )  =  af .. 


h 


P\  71  ’  r- 

V<— <  al\..ai  j 


=  af...a;o(n„n2) 


(1) 


where  Hj 


Py 

aj 


and  n. 


T  hese  "  n  ”  variables  are  sometimes  called  simildrity  variables. 
T  his  is  because  two  natural  systems,  with  different  values  of  as 
and  bs  but  the  same  value  of  n,  will  tend  to  have  similar 
emergent  behaviour.  An  example  is  the  flow  of  air  past  an  air¬ 
craft  in  the  atmosphere  or  past  a  model  of  the  aircraft  in  a 
wind  tunnel.  The  n  variable  in  this  case  is  the  Reynolds 
number.  If  this  is  the  same  in  both  cases,  then  the  model  in 
the  wind  tunnel  will  give  results  relevant  to  the  full-scale  air¬ 
craft  in  the  atmosphere. 

Self-similar  solutions  correspond  to  problems  where  the  values 
of  the  variable  (hj  for  example)  tend  to  zero  or  infinity.  T  hree 
possibilities  are  available  (see  Reference  [1]  for  further 
discussion): 

Type  1  metamodel,  o  tends  to  a  non-zero  finite  limit  as  one 
of  its  arguments  tends  to  zero  or  infinity.  T  his  means  that  in 
most  practical  cases,  this  argument  can  be  eliminated  from  the 
relationship,  giving  a  simplified  form  in  equation  (1). 

Type  2  metamodel,  o  has  power-law  asymptotics  of  the 
form: 


(2) 
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as  the  argument  tends  to  zero  or  infinity.  In  theoretical 
physics,  these  systems  are  examined  from  a  gauge  theory  point 
of  view  using  a  "renormalisation  group”  approach  in  which 
the  parameter  is  considered  (through  the  repeated  applica¬ 
tion  of  a  renormalisation  group)  at  larger  and  larger  (or 
smaller  and  smaller)  gauges,  giving  an  asymptotic  expression 
of  the  form  required  for  a  type  1  or  2  metamodel.  Later,  we 
will  discuss  the  renormalisation  group  in  more  depth  and 
relate  it  to  concepts  of  control  of  the  battlespace. 

Type  3  metamodel.  Neither  1  nor  2  holds  and  self-similarity 
is  not  observed;  o  has  no  finite  limit  different  from  zero  and 
no  power- 1  aw  asymptotics. 

T  his  approach  directs  us  (for  evidence  of  metamodels  of  types 
1  and  2)  to  search  for  evidence  of  power-law  relationships  of 
the  form  y  =  x" ,  which,  if  plotted  on  a  log-log  scale,  give  a 
straight  line  whose  slope  is  the  power-law  exponent.  Such 
expressions  arise  naturally  in  certain  types  of  complex  systems, 
particularly  where  fractal  structures  are  involved,  and  are 
referred  to  as  seeling  rdetionships  since  they  have  no  preferred 
gauge.  E  vidence  of  such  scaling  relations  is  thus  &idence  in  support  of  the 
assumption  of  relative  gauge.  Chapter  3  shows  that  there  is  clear 
evidence  for  such  an  assumption  in  historical  conflict  data,  and 
that  metamodels  of  types  1  and  2  should  be  expected  for 
agent-based  models  of  conflict.  In  addition,  we  should  expect 
to  see  evidence  for  normalised  "scaling  collapse"  as  exempli¬ 
fied  by  the  function  0(0;, Hj)  in  equation  (2),  and  we  should 
expect  to  see  the  effect  of  renormalisation  groups. 

If  we  consider  an  agent-based  "distillation"  such  as  the  ISAAC 
model  developed  under  Project  Albert  by  the  U.S.  Marine 
Corps  Combat  Development  Centre,  we  can  consider  the 
emergent  behaviour  of  such  a  model  in  terms  of  both  the  spa- 
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tial  clustering  of  the  agents,  and  the  attrition  that  they  inflict 
on  the  opponent.  For  such  a  distillation,  a  metamodel  of  type  2 
applies[l]  that  allows  us  to  relate  the  attrition  rate  for  one  side 
to  the  clustering  dynamics  of  the  opposing  side,  as  measured 
by  the  mean  fractal  dimension  of  these  clustering  agents.  As  a 
simple  example,  (given  in  [1]),  assume  that  the  command  pro¬ 
cess,  say  for  Red,  is  represented  by  the  following  effects: 

1.  The  number  of  discrete  clusters  of  Red  agents  at  timet, 

A/ ft),  is  specified  ahead  of  the  simulation. 

2.  A/ ft)  is  a  decreasing  function  oft. 

T  hese  assumptions  are  meant  to  suggest  that  the  number  of 
Red  clusters  decreases  in  time,  reflecting  the  desire  of  Red  to 
concentrate  force.  With  these  assumptions  let  us  further 
assume  that  the  smallest  cluster  of  Red  agents,  X  ft),  at  timet,  is 
taken  and  added  to  another  randomly  chosen  cluster  of  Red 
agents.  T  his  process  thus  represents  both  the  concentration  of 
Red  force  and  the  reconstitution  of  force  elements. 

Let  us  now  define  y?(x,0  =  (expected  number  of  clusters  of  Red 
agents>sizex  at  timet)/  (initial  total  number  of  clusters  of  Red 
agents)  and  A/ft)=(the  total  number  of  remaining  clusters  of 
Red  agents  at  time  t)/  (initial  total  number  of  clusters).  Given 
the  assumptions  and  definitions  above,  it  can  then  be  shown 
that  (p{x,t) ,  the  cumulative  distribution  of  cluster  sizes  at  time 
f,  approaches  a  self-similar  distribution  as  time  progresses 
(i.e.,  a  scaling  collapse  takes  place).  Thus  the  cluster  size  distri¬ 
bution  evolves  over  time  by  a  scaling  relation.  (p{x,t)  can  then 
be  represented  in  the  self-similar  form: 


(p{x,t) 


g{xl  X{t)) 


X{t) 
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where  g  is  some  positive  continuous  function  and  where 
g(l)=A/ ft)X  ft).  T  his  self-similar  form  means  that  we  can  define 
the  distribution  of  relative  cluster  size  in  a  way  that  is  time 
invariant  (although  the  actual  cluster  sizes  will  change). 

Now  assume  that  the  evolution  of  the  distribution  of  ^ix,t)  is 
smooth  (a  small  change  in  time  t  leads  to  a  small  change  in 
(p{x,t)-,  this  is  equivalent  to  saying  that  the  renormalisation 
group  issmooth  [1]).  If  ^{x,i\-S)t)  isthe expected  cluster  size 
X  at  time  and  ^{x,t)  isthe  same  expectation  at  timet, 
then  this  assumption  means  we  can  find  a  constant  b  to  first 
order  such  that: 


(p{x,  ^)  =  (1  -I-  bS)(p(x,  (1  -  S)t) . 

It  then  follows  that: 


log  (p{x,  t)  =  b  log  t  +  c 

for  some  constant  c  and  the  normalised  expected  cluster  size  at 
timet,  varies  as  a  power-law  with  increasing  timet  and 
scaling  constants. 

If  AB  isthe  change  in  the  number  of  Blue  agents,  Lauren  [3] 
has  shown  that: 


At 

is  proportional  to  the  product  of  (Red  unit  effect! veness)x (the 
probability  of  meeting  a  Red  cluster)x(the  expected  number  of 
Red  units  per  cluster).  It  is  assumed  that  unit  effectiveness  is 
constant.  Keeping  the  cluster  size  constant  for  the  moment, 
this  indicates  [3]  that  the  rate  of  change  of  Blue  agents  is  given 
by  an  expression  of  the  form: 
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where  D  is  the  average  fractal  dimension  of  Red  (and  therefore 
an  indication  of  how  Red  clusters/ collaborates  locally)  and 
both  r  and  q  are  exponents.  This  equation  is  a  form  of 
Lanchester  law  where  the  rate  constant  is  dependent  upon  the 
clustering  of  Red  agents.  If  Red  cluster  size  varies  according  to: 


^(x,0 


g(x/X(t)} 

x(0 


wherex  is  the  cluster  size  and  X  (t)  the  smallest  cluster  at  time  t, 
then  we  can  write: 


where  A/ ft),  inversely  related  to  X  ft),  is  the  normalised  number 
of  clusters  of  Red  at  timet  and  gfyftjj  is  the  mean  of  the  distri¬ 
bution  of  cluster  size,  which  evolves  as  a  power  law  (as  we  have 
shown  in  certain  cases). 

More  General  Distribution  of  Cluster  Sizes 

For  self-organising  groups  of  agents  that  are  approaching  a 
critical  point  in  the  form  of  a  Bak-Sneppen  evolution  model, 
but  are  not  necessarily  attempting  to  concentrate  force  or  to 
reconstitute  force  in  the  sense  described  above,  our  previous 
analysis  in  Chapter  1  indicated  that  the  distribution  of  cluster 
sizes  at  some  intermediate  point  s  is  given  by: 
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where  is  the  minimal  site  probability  at  timestep  s.  As 
noted  before,  this  converges  to  a  power-law  distribution  ass 
tends  to  infinity. 

In  [1],  we  also  relate  the  cluster  fractal  dimension  predicted  by 
the  Bak-Sneppen  evolution  model  of  local  species  coevolution, 
to  that  which  emerges  from  initial  experiments  with  the 
ISAAC  cellular  automaton  model.  (For  a  description  of 
I SAAC  see  the  Web  site  at  reference  [4].) 

The  Renormalisation  Group 

The  gauge  invariant  approach  to  metamodelling  outlined 
above  also  leads  us  to  consider  the  role  of  the  renormalisa¬ 
tion  group,  which  explicitly  appears  in  terms  of  its  effect  on 
the  distribution  function  0(ni,n2)  in  our  characterisation  of 
metamodels.  Let  us  explore  this  a  little  more  here.  Suppose 
that  u{x,t)  is  a  function  of  the  two  variables x  and  t  Using 
the  "static  scaling"  assumption  of  renormalisation  [5],  we 
assume  that  we  have  a  group  of  renormalisations  of  the  form 
so  that: 


Z(o) 

From  the  group  property,  =  R,,^^ 

It  follows  that  Z{a)Z{b)  =  Z{ab) 

Thus  Z{b)  =  b“  for  some  exponent  a 
If  u\x,t)  is  a  fixed  point  of  the  renormalisation,  then; 
u\x,t)  =  b~“u\b‘^x,bt)  Vh 

Choose  b  =  -,  then; 
t 

u  {x,t)  =  fu 
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The  repeated  application  of  the  renormalisation  process  will 
thus  (in  the  limit)  produce  a  functional  relation  of  the  form  u. 
This  explains  why  the  type  2  metamodel  has  the  form 
assumed,  and  also  why  the  (renormalisation  invariant)  cluster 
size  distribution  of  Carr  and  Pego  discussed  earlier  has  the 
form  derived  by  them.  Note  the  simplification  that  has  been 
achieved  in  going  from  an  unknown  function  of  two  variables 
to  a  normalised  unknown  function  of  only  one  variable. 

T  he  ability  of  a  force  to  control  an  area  of  operations  can  be 
related  to  the  fractal  dimension  of  the  force  through  the  use  of 
such  a  renormalisation  process  as  we  now  show. 

Control  of  the  Battlespace 

Looking  now  at  the  phenomenon  of  control  of  the  battlespace, 
we  can  consider  the  problem  using  the  renormalisation 
approach  as  in  a  type  2  metamodel.  Consider,  as  shown  in  Fig¬ 
ure  4.1,  the  Area  of  Operations  (AO)  of  a  military 
commander.  For  simplicity  we  assume  this  is  a  square  of  side  L . 


F igure  4.1:  A  red  of  0 peretions 
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We  assume  that  the  commander  aims  to  establish  control  in 
this  area.  Firstly  we  have  to  define  what  this  means.  Each  unit, 
shown  by  a  dot  in  Figure  4.1,  has  an  area  surrounding  it  that  it 
can  control.  The  size  of  this  area  is  defined  by  the  nature  of  the 
force  and  its  associated  sensors  [6],  giving  rise  to  a  "bubble”  of 
control  around  the  unit.  I  n  two  dimensions,  let  this  area  corre¬ 
spond  to  a  square  of  side  /.  We  assume  that  /  is  significantly 
smaller  than  the  dimensions  of  the  AO.  (Note:  we  have 
restricted  the  battlespace  here  to  two  dimensions  to  simplify 
the  discussion.  Flowever,  the  same  approach  should  work  in 
three  dimensions,  corresponding  to  the  complete  battlespace.) 

Now  letD  be  the  fractal  dimension  of  the  force  under  the  com¬ 
mander's  control  within  the  AO .  Suppose  we  partition  the  AO 
into  square  cells  of  width  /.  Let  A/  be  the  total  number  of  such 
cells,  so  that  =  ^ .  Let  A/  (I)  equal  the  number  of  cells  in  the 
AO  that  are  occupied  by  one  of  the  units  making  up  the  force. 
By  definition  of  the  fractal  dimension,  we  have  that  N(l)  =  r 
(normalising  the  constant  of  proportionality  to  1).  If  p  is  the 
probability  that  a  cell  chosen  at  random  in  the  AO  is  under 
control,  then: 


_  N(i)  _ 

N  ~  N  ~  A  ' 

Note  that  D  always  lies  between  0  and  2,  so  that  p  is  well 
defined. 

In  discussion  with  senior  UK  commanders  who  have  had 
recent  operational  experience  at  a  high  level,  the  concept  of 
control  of  an  area  as  corresponding  to  the  prevention  of  flow 
through  an  area  (flow  in  terms  of  an  opposing  force,  or  per¬ 
haps  some  third  party)  has  been  endorsed  as  a  good  analogy. 
We  thus  define  the  commander  as  having  "weak  control"  of 
his  area  of  operations  if  he  can  to  some  extent  control  move- 
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merit  through  the  AO.  This  is  similar  to  the  problem  of 
determining  whether  fluid  can  seep  through  a  block  of  semi- 
porous  rock,  as  discussed  in  [7]  where  a  renormalisation  group 
approach  was  used.  We  thus  define  weak  control  as  corre¬ 
sponding  to  a  span  of  controlled  areas  that  stretch  either  from 
side  to  side  or  top  to  bottom  of  the  AO.  Following  on  from 
this,  we  define  the  commander  to  have  "strong  control"  of  the 
AO  if  there  is  a  span  of  controlled  areas  stretching  both  from 
side  to  side  end  top  to  bottom  of  the  AO,  resulting  in  a  strong 
constraint  on  the  flow  of  people  through  the  area. 

T  he  question  at  issue  is  then:  how  do  these  concepts  of  control 
relate  to  the  ability  of  the  force  to  collaborate  locally  (the  frac¬ 
tal  dimension)? 

Consider  first  a  cell  of  four  elements  where  each  cell  is  a 
square  of  side  /,  which  a  single  unit  can  control.  We  now  con¬ 
sider  the  probability  Pi  of  weak  or  strong  control  ofthissquare 
cell  of  side  2/  in  terms  of  the  probability: 

^2-D 


of  a  unit  controlling  each  of  the  squares  of  side  I.  We  consider 
each  of  the  five  different  classes  of  configuration  for  this  cell,  as 
shown  in  Figure 4.2. 

I  n  Figure  4.2,  we  show  the  five  classes  a  to  e  of  configuration, 
and  mark  beside  each  case  whether  this  gives  weak  or  strong 
control,  by  considering  the  span  of  controlled  areas. 
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«  I  I  I  0  weak.  0  strong 

''  ffl  ffl  ffl  E  0»,,,k.0«»ng 

‘  EB  3  ffl  B  E  ffl  4w...,k.()„r„ng 

4  weak,  4  strong 

e  1  weak.  1  strong 

Figure4.2:  FiveConfigurdtion  Classes 

T  he  probability  of  each  configuration  can  be  derived  in  terms 
ofp.  For  example,  the  probability  of  any  of  the  cases  in  config¬ 
uration  d  is  p^(l-p).  By  adding  up  the  configurations 
corresponding  to  weak  control,  and  taking  into  account  the 
probability  of  each  such  configuration,  we  have  the  relation: 

Pi(weak)  =  4p^(l-py  +4p^(l-p)  +  p'' . 

We  can  do  the  same  thing  for  strong  control,  leading  to  the 
relation: 

p^ (strong)  =  4p^(\- p)  +  p‘^ . 

Using  the  renormalisation  group  approach  [5],  we  iterate  at 
increasing  levels  of  cell  size,  leading  to  the  relations: 

•  Weak  control: 

•  Strong  control:  ;7„^(4-3;7„) 

T  hese  give  rise  to  the  recursive  schemes  shown  in  Figures  4.3 
and  4.4.  T  he  relations  for  weak  and  strong  control  above  cor¬ 
respond  to  the  relationships  respectively: 

f{x)  =  (4-4x4- x^) 

g(x)  =  x^(4-3x) 
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y 


y 


Figure4.4:  Plot  ofy  =  g(x)  dncly  =  x  ■  strong  control 

T  he  fixed  points  in  the  recursive  relation  of  weak  control  cor¬ 
respond  to  the  values X  shown  in  Figure  4.3  such  that  y  =  f(x) 
intersects  y=  x.  Similarly  for  strong  control,  the  fixed  points 
correspond  to  the  values  x  such  that  y  =  g(x)  intersects  y  =  x. 
For  both  weak  and  strong  control,  there  are  stable  fixed  points 
at  X  =  0  and  X  =  1.  FI  owever,  there  is  also  an  unstable  fixed 
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point  between  these  that  is  different  for  strong  and  weak  con¬ 
trol.  This  was  calculated  to  be  0.382  for  weak  control,  and 
0.768  for  strong  control. 

Control  and  Fractal  Dimension 

In  either  case,  starting  with  a  given  fractal  dimension  for  the 
force,  and  thedimensionsof  the  AO,  we  can  calculate  a  corre¬ 
sponding  starting  probability: 

l2-0 


For  side  length  L  of  the  AO  there  will  be  a  corresponding 
value  of  iteration  order  n  such  that  ti  =  l.  By  using  the  recur¬ 
sive  scheme  above,  we  can  calculate  for  this  value  of  n  the 
corresponding  probability  of  weak  or  strong  control  of  the 
AO.  Consideration  of  Figures 4.3  and  4.4  indicates  that  there 
is  a  critical  value  of  the  probability: 

12-D 


such  that  values  above  Po  polarise  towards  very  good  control, 
whereas  values  below  Po  polarise  towards  very  poor  control.  In 
fact,  the  point  Po  corresponds  to  a  phase  change  in  the  behav¬ 
iour  of  such  a  system. 

Examination  of  Figures  4.3  and  4.4  shows  that  it  is  easier  to 
iterate  towards  good  weak  control  than  towards  good  strong 
control,  as  we  would  expect  (since  weak  control  is  easier  to 
achieve  than  strong  control).  Figure  4.5  shows  how  this  itera¬ 
tion  works  for  a  starting  probability  of  0.65  and  the 
requirement  of  weak  control. 
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Lockout 

From  a  game  theoretic  perspective,  we  can  see  that  each  side  is 
trying  to  drive  its  own  value  of  control  up,  and  the  other  side's 
down.  T  he  analysis  above  indicates  that  there  should  be  rapid 
lockout,  i.e.  one  side  should  rapidly  gain  control  and  lock  the 
other  side  out. 

Percolation  Theory  and  the 
Renormalisation  Group 

In  relating  these  ideas  to  the  behaviour  of  natural  systems,  the 
various  configurations  of  2x2  controlled  areas  are  identical  to 
the  porous  and  nonporous  regions  in  semipermeable  rock 
structures  [7].  T  he  study  of  such  processes  is  referred  to  as  Per- 
coldtion  T  heory.  A  good  step-by-step  introduction  to  the  theory 
and  someworking  examples  of  how  such  processes  work  in  two 
dimensions  are  provided  at  the  Web  site  reference  [8].  If  p  is 
the  probability  of  an  individual  rock  (or  crystal)  domain  being 
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porous,  then  p  is  the  driving  parameter  of  the  process.  The 
behaviour  can  be  tuned  in  the  sense  that  we  discussed  in 
Chapter  2.  Ifp  is  below  the  critical  value Po,  then  the  clusters 
are  not  large  enough  to  form  a  path  of  percolation  from  one 
side  of  the  structure  to  the  other  (popping  noise).  When  the 
probability  p  is  above  the  critical  value  Po,  then  suddenly  all 
clusters  tend  to  spread  from  one  side  to  the  other,  allowing 
percolation  throughout  the  structure  (these  are  thus  called  per- 
coldtion  clusters)  and  corresponding  to  a  phase  change  in  the 
dynamic  of  the  system  (snapping  noise).  Near  to  the  critical 
point,  it  can  be  shown  [8]  that  the  distribution  in  size  of  perco¬ 
lating  clusters  is  a  power-law,  corresponding  to  a  fractal 
distribution  of  cluster  size  (crackling  noise).  In  fact,  for  p 
slightly  greater  than  Po,  the  fraction  F  of  individual  domains 
that  form  part  of  a  percolating  cluster  takes  the  form 
F  =  Fo(p-poy .  From  this,  we  can  see  that  as  p  approaches 
the  critical  value  Po  from  above,  the  fraction  of  domains  form¬ 
ing  part  of  the  percolating  cluster  tends  to  zero  (for  a  very  large 
initial  configuration).  Thus  such  clusters  can  become  very  ten¬ 
uous  close  to  the  critical  point. 

These  configurations  come  originally  from  attempts  to  model 
lattices  of  magnetic  spins  in  more  than  one  dimension  [5]. 
Such  2x2  configurations  are  then  referred  to  as  block  spins 
since  they  are  composed  of  four  individual  spins,  each  of 
which  may  be  either  up  or  down  (the  block  spin  is  defined  as 
the  sum  of  the  signs  of  the  individual  spins,  so  it  still  has  the 
value  -1-1  or  -1).  By  developing  the  idea  of  the  renornidllsdtlon 
group,  iterating  to  larger  and  larger  domains,  Wilson  (building 
on  work  by  K  adanoff),  was  able  to  show  that  such  arrays  of 
magnetic  spins  do  indeed  exhibit  phase  change  effects  (as  we 
have  shown  above  for  the  phase  change  between  being  out  of 
control  and  being  in  control  of  a  region),  and  that  the  param- 
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eters  involved  can  be  calculated  explicitly.  For  this  work, 
Wilson  was  awarded  the  Nobel  Prize  for  Physics. 

Implications  FOR  Self-Organising 
Information  Networks 

If  a  Self-0  rganising  Information  Network  is  thought  of  as  a 
grid  of  connections,  and  we  have  a  probability  p  of  creating  a 
link  between  one  element  of  the  network  and  another,  then 
percolation  across  the  network  corresponds  to  being  able  to 
send  a  signal  from  one  end  of  the  network  to  the  other.  From 
Percolation  Theory,  we  can  thus  see  that  we  would  expect  a 
phase  change  in  the  dynamic  of  such  a  network.  If  p  is  small, 
then  only  local  connections  can  be  made.  Flowever,  at  some 
critical  value  of  p,  there  will  be  a  phase  change  such  that  con¬ 
nections  across  the  full  network  can  suddenly  be  made.  Near 
to  this  critical  value  of  p,  the  clusters  formed  by  those  con¬ 
nected  on  the  network  will  form  a  fractal  set,  and  the 
distribution  of  such  cluster  sizes  is  described  by  a  power-law 
relationship  between  size  and  frequency  of  that  size.  Two 
questionsthat  arise  are: 

1.  What  is  the  benefit  (and  cost)  of  being  above  the  critical 
threshold  so  that  the  connections  are  robust? 

2.  Flow  can  we  measurethe  benefit  of  using  the  knowledge 
obtained  by  such  networked  interconnection? 

Recent  work  by  Perry  [9]  has  exploited  the  idea  of  information 
entropy  to  address  the  second  question  (with  a  reduction  in 
entropy  across  the  network  corresponding  to  an  increase  in 
knowledge,  and  this  then  being  equivalent  to  a  reduction  in 
delay  in  prosecuting  an  action).  To  understand  where  this  idea 
originates,  we  first  look  at  the  influence  of  knowledge  in  war- 
games  from  an  open  systems  perspective. 
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VUexrgamesasOpen  Systems  Sustained  by 
Knowledge  Flowing  Across  the 
Boundary 

As  noted  recently  by  Roske  [10],  a  wargame  isan  open  system 
of  the  type  introduced  in  Chapter  1.  As  he  notes, 

"In  a  cl  dssic  commend  post  exercise,  we  Inject  humen  decislon- 
meking  Into  a  structured  aivironmmt  to  generete  open  system 
behdviours.  H  umen  decisionmaking  represents  energy  crossing 
the  structured  system  boundary...  " 

In  [11],  Perry  and  M  offat  developed  the  idea  of  using  Infor¬ 
mation  E  ntropy  (from  Shannon  I  nformation  T  heory)  as  a  way 
of  capturing  the  knowledge  available  to  a  military  com¬ 
mander's  decisionmaking  during  a  wargame.  T  his  was  initially 
applied  to  looking  at  the  benefit  to  be  gained  from  advanced 
airborne  standoff  radars  (such  as  J STARS  in  the  U.S.  or 
ASTOR  in  the  UK). 

A  series  of  wargames  was  carried  out  in  the  UK  to  quantify  the 
difference  in  combat  effectiveness  of  a  force  without  airborne 
standoff  radar  (ASTOR)  in  comparison  with  a  force  with 
ASTOR,  or  a  force  with  other  weapon  systems  whose  life- 
cycle  costs  approximately  equalled  those  of  ASTOR.  The 
common  thread  through  all  of  the  cases  to  be  examined  was 
the  stream  of  decisions  made  by  the  friendly  force  com¬ 
mander.  We  describe  here  how  we  were  able  to  capture  the 
flow  of  information  to  the  Blue  commander  during  the  war- 
games  using  the  concept  of  Information  Entropy,  and  then 
turning  that  into  a  measure  of  "knowledge."  This  approach 
allowed  us  to  measure  the  quantity  of  knowledge  flowing 
across  the  boundary  of  the  system  in  order  to  influence  the 
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decisions  made  by  the  Blue  commander  and  the  onward  evo¬ 
lution  of  the  wargame. 

Wargame  Structure 

T  here  are  three  major  problems  with  the  use  of  wargames  to 
support  military  studies:  (1)  too  little  output  data,  (2)  the  likeli¬ 
hood  of  atypical  results,  and  (3)  oversimplification.  The  first 
problem  stems  from  the  fact  that  wargames  are  generally  slow, 
cumbersome,  and  resource  intensive.  Consequently,  most  ana¬ 
lysts  who  use  them  to  support  studies  plan  only  a  small  number 
of  games,  thus  precluding  significant  statistical  results.  The 
second  problem  recognises  the  possibility  that  the  sequence  of 
decisions  taken  by  the  players  in  these  games  represents sfaf/sf/- 
cal  outliers.  Players  may  adopt  extreme  strategies  that  exist 
"outside”  of  what  is  considered  to  be  a  typical  military 
response.  The  third  problem  reflects  the  fact  that  human  play¬ 
ers  can  only  approximate  the  results  of  combat  operations.  I  n 
our  studies,  we  addressed  these  problems  in  three  ways:  by 
arguing  that  our  wargames  are  quasi -memoryless  processes  for  tac¬ 
tical  situation  assessment;  by  introducing  the  epitomising  strategy 
principle  in  wargames;  and  by  embedding  computer  models  to 
adjudicate  engagements  in  the  manual  games.  We  discuss  the 
first  two  of  these  concepts  next. 

The  M  emoryless  Property  of  Wargames 

The  wargames  used  to  support  this  study  were  two-sided, 
zero-sum  games  played  over  several  discrete  time  periods  or 
cycles  (Bowen  [12]).  Consequently,  the  entire  campaign  can 
be  viewed  as  a  dynamical  process  in  which  the  state  variables 
are  the  force  levels  on  both  sides.  In  each  of  the  wargames, 
both  the  friendly  and  enemy  commanders  formulated  opera- 
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tions  plans  based  on  the  stated  campaign  objectives.  During 
the  play  of  each  game,  both  commanders  assessed  the  overall 
situation  at  the  strategic  and  tactical  levels.  Assessment  at  the 
strategic  level  consisted  of  examining  the  need  to  alter  the 
campaign  plan.  At  the  tactical  level,  it  supported  force  alloca¬ 
tion  decisions  consistent  with  the  implementation  of  the 
campaign  plan. 

A  dynamic  process  is  said  to  be  memoryless  or  M  erkovien  if  at 
each  cycle,  the  state  of  the  system  is  influenced  only  by  the 
state  of  the  system  in  the  previous  cycle,  and  not  by  the  specific 
history  of  the  system  (Stark  and  Woods  [13]).  The  wargames 
conducted  to  support  this  study  epproximetely  satisfied  these 
conditions  as  follows: 

•  Tactical  Decisions:  In  all  of  the  wargames  played, 
the  strategic  situation  was  such  that  neither  commander 
found  it  necessary  to  alter  their  original  plan.  Conse¬ 
quently,  the  commanders'  decisions  centred  on  the 
allocation  of  their  forces  T  his  forced  them  to  focus 
exclusively  on  the  tactical  situation  in  the  current  game 
cycle  and  their  assessmentsof  the  likely  situation  in  sub¬ 
sequent  cycles.  T ransitions  in  the  state  variables 
therefore  depended  upon  the  status  of  the  forces  at  the 
end  of  the  previous  cycle,  the  decisions  taken  in  the  cur¬ 
rent  cycle,  and  the  combat  attrition  experienced  in  the 
current  cycle. 

•  Cycle  Independence:  Both  commanders  were 
assumed  to  act  outside  the  opponent's  decision  cycle. 

T  hat  is,  the  opponent  was  assumed  to  have  sufficient 
time  within  each  cycle  to  redeploy  his  forces  so  that  intel¬ 
ligence  on  unit  identity,  type,  and  location  gained  in  the 
previous  cycle  was  no  longer  valid  in  the  current  cycle; 
past  history  had  no  effect  on  the  commanders'  current 
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decisions.  In  actual  practice,  we  found  this  assumption 
only  partially  valid,  as  will  be  made  more  apparent  below 
in  the  discussion  of  the  UK  FAST  HEX  gaming  model 
used  in  the  wargames  for  this  study. 

U  nder  these  conditions,  the  campaign  can  be  viewed  as  a 
sequence  of  decisions  taken  in  a  discrete  dynamical  system. 
Each  commander  attempts  to  select  a  set  of  decisions  (strategy) 
that  will  maximise  the  likelihood  that  he  will  achieve  his  mis¬ 
sion  and  that  is  consistent  with  the  overall  campaign  plan. 
Because  the  decisions  are  made  under  a  connecting  campaign 
plan,  we  assert  that  the  process  is  only  quasi- memoryless.  This 
assumption  is  most  relevant  when  the  situation  on  the  ground 
is  in  a  state  of  rapid  flux-the  most  interesting  case. 

The  Epitomising  Strategy  Principle 

In  the  play  of  the  games,  we  strove  to  ensure  that  the  com¬ 
manders  took  actions  that  epitomised  the  side's  historic 
conduct  in  battle.  We  attempted  to  avoid  bold,  daring,  bril¬ 
liant  manoeuvres  as  well  as  bungling  incompetence.  For  the 
purposes  of  analysis,  more  cautious  conservative  strategies  that 
are  consistent  with  accepted  doctrine  are  preferred.  T  here  is  a 
real  danger  that  uncontrolled  play  would  have  resulted  in  pro¬ 
ducing  only  outlier  results,  when  what  is  wanted  are  typical 
results.  We  achieved  this  by  playing  the  wargames  "open."^ 
Red  and  Blue  players  were  able  to  discuss  strategies  and  deci¬ 
sions  in  the  presence  of  a  game  controller  and  thus  we  ensured 
that  the  actions  contemplated  by  either  did  not  constitute 
extraordinary  tactics. 


wargameis"open  to  Blue"  if  the  physical  state  of  Red  isfully  known  to 
the  Blue  player  (Bowen  [12]). 
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In  keeping  with  this  principle,  great  pains  were  taken  to  ensure 
that  each  side  "knew”  only  those  things  that  would  normally 
be  known  through  the  available  surveillance  and  intelligence 
assets.  Players  were  forced  to  examine  the  information 
received  through  sensors  and  surveillance  assets  carefully  to 
ensure  that: 

1.  Sound  militaryjudgement  was  used  in  considering  the 
decision  options  available. 

2.  Players  used  what  they  learned  and  interpreted  from 
sensor  reports,  and  not  what  they  saw  on  the  "game 
board." 

3.  The  information  received  was  consistent  with  the  limita¬ 
tions  and  capabilitiesof  the  equipment  being  used,  the 
employment  of  the  surveillance  assets,  and  other  envi¬ 
ronmental  conditions. 

VUexrgaming  with  fast  hex 

I  n  the  play  of  the  FAST  H  EX  wargame,  the  continuity  of  the 
battle  was  modelled  as  a  series  of  discrete  timesteps  referred  to 
as  "game  cycles"  (Figure  4.6).  T  he  length  of  a  cycle  can  be  set 
by  the  players,  but  is  usually  chosen  to  be  2  hours.  Within  each 
cycle,  a  linear  sequence  of  actions  is  taken  and  the  conse¬ 
quences  of  each  is  evaluated  to  simulate  events  within  the 
cycle.  Blue  and  Red  alternately  begin  the  sequence  in  order  to 
smooth  out  the  advantages  or  disadvantages  of  "going  first." 

T  he  flow  chart  in  Figure  4.6  can  be  thought  of  as  a  series  of 
modules  in  a  fully  automated  simulation  of  combat,  less  the 
decision  module.  The  modules  are  rather  simple  representa¬ 
tions  requiring  extensive  player  input.  For  example,  the  player 
selects  the  type  of  reconnaissance  system  to  be  used,  states  the 
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current  environmental  conditions,  and  the  hexagon^  on  the 
game  board  to  be  searched.  With  this  information,  the  model 
applies  the  appropriate  probability  of  detection  or  recognition 
and  reports  the  results.  Environmental  conditions  and  terrain 
features  are  not  known  by  the  model  and  are  inserted  through 
rulesconstraining  the  play. 


Ad\  ance  to 
tlie  next  cycle 


Figure  4. 6:  FAST  H  EX  Game  Cycle  Sequence 

The  Decision  Problem 

With  the  epitomising  principle  in  mind  and  the  constraints  on 
gaming  imposed  by  the  FASTH  EX  model,  we  examine  the 
decisions  open  to  each  commander  with  particular  emphasis 
on  the  Blue  commander.  Great  care  has  been  taken  to  ensure 
that  decisions  required  in  actual  combat  had  their  equivalent 


^FASTH  EX  uses  a  hexagonal  game  board  much  like  I  DA  H  EX  .  For  these 
games,  each  hexagon  is  7.5  km  from  face  to  face. 
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during  the  play  of  the  wargame.  T  his  is  made  possible  by  the 
fact  that  players  can  override  almost  all  computer-generated 
outcomes  in  the  game.  Therefore,  realism  can  be  imposed  in 
those  cases  where  the  model  obviously  strays. 

The  commander  makes  operational  and  tactical  decisions  at 
each  combat  phase  in  the  FAST  HEX  wargame  cycle  in  keep¬ 
ing  with  his  overall  operational  and  tactical  aims  (Figure  4.6). 
T  he  one  exception  of  course  is  the  engagement  phase:  the 
engagement  is  a  consequence  of  the  decisions  taken  by  the 
Red  and  Blue  commanders  in  the  other  phases.^  T  herefore  at 
each  cycle,  the  decision  set  in  FAST H  EX,  D(t),  consists  of  4 
components,  D(f)={di(f),d2(f),d3(f),d4(f)},  where  t  indicates  the 
cycle  and  d, ft)  is  the  decision  taken  atthe/th  phase  in  theFAS- 
T  H  EX  model.  T  he  following  is  a  description  of  the  decisions 
taken  at  each  phase: 

1.  Reconnaissance,  di(f):  The  commander  must  decide 
which  of  the  reconnaissance  assets  at  his  disposal  to  allo¬ 
cate.  In  most  cases,  this  means  deciding  whereto  direct 
his  sensors  and  how  many  to  commit  to  the  process  (tak¬ 
ing  account  of  higher  level  assets  such  as  satellite 
surveillance).  It  should  denoted  here  that  reconnais¬ 
sance  assets  are  used  primarily  to  identify  valid  targets 
and  main  force  concentrations. 

2.  Strike,  d2(f):  The  strike  decisions  can  be  thought  of  as 
the  allocation  of  deep  fire  assets.  T  he  engagement  phase 
adjudicates  the  close  battle,  but  the  deep  fire  battle  is 
handled  separately. 


^Engagement  adjudication  is  done  using  look-up  tables  based  on  lower 
level  computer  modelling  of  the  various  types  of  engagement. 
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3.  M  OVEMENT,  dslf):  T  he  commander  must  decide  which 
of  his  units  to  move  and  how  far  they  are  to  move.  He  is 
constrained  by  terrain,  the  maximum  speed  of  his  units, 
and  the  degree  to  which  the  units  are  fit  (in  terms  of 
damage  inflicted)  to  accomplish  movement.  U  nits  move 
at  the  lowest  level  of  resolution:  the  battalion,  battle 
group,  or  squadron. 

4.  Post-Engagement  Moves,  d4(f):  After  each  engage¬ 
ment,  the  commander  assesses  the  damage  done  to  his 
units.  If  the  units  are  not  capable  of  continuing  as  an 
integral  force,  they  can  be  withdrawn,  consolidated  with 
other  forces,  or  both. 

Optimal  Control  Formulation 

Each  of  the  decisions  is  taken  from  among  a  discrete  set  of  the 
possible  choices  described  above  and  therefore  the  set  {D(t)} 
consisting  of  all  possible  decisions  at  cycle  t  has  cardinality 
equal  to  the  product  of  the  cardinality  of  the  4  phase-decision 
sets.  The  collection  of  decision  sets  at  each  cycle  in  the  game  is 
then  referred  to  as  the  decision  stream  for  that  game  and  there¬ 
fore,  we  denote 


the  decision  stream  for  an  m-cycle  game.  Clearly,  the  number 
of  possible  decision  streams  for  even  a  simple  wargame  can 
easily  become  unmanageably  large  and  thus  we  are  burdened 
by  the  "curse  of  dimensionality." 

T  he  consequences  of  the  commander's  decisions  at  each  cycle 
can  be  measured  in  several  ways.  As  mentioned  earlier,  we 
have  selected  the  status  of  friendly  and  enemy  forces.  Conse- 
quently,ifwelet 
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be  the  vector  of  combat  strength  of  the  b  friendly  force  compo¬ 
nents  and  the  r  enemy  force  components  at  cycle  t,  then  the 
decision  stream  can  be  viewed  as  a  memoryless  multistege  decision 
process  as  depicted  in  Figure  4.7.  At  each  cycle,  the  commander 
wishes  to  select  D  ft)  so  that  a  performance  function  dependent 
upon  the  vector  X  ft)  is  optimised  in  some  way. 


DIO)  D{I)  Dim-l) 


Figure 4. 7:  k  Wargameas  an  Open  Dynamical  Process 

In  this  formulation,  is  a  transition  function  so  that 
x{t+\)=eXx{t),D{t)\  The  status  of  both  Red  and  Blue  forces 
in  terms  of  combat  strength  in  cycle  t-Fl  depends  upon  the 
their  combat  status  in  cycle  t  and  the  decisions  made  in  cycle  t. 
T  he  initial  condition  X  (0)  is  the  total  combat  strength  of  the 
friendly  and  enemy  forces  at  the  beginning  of  the  campaign, 
and  X(m)  is  the  status  of  both  at  the  termination  of  the 
campaign. 

T  he  commander's  problem  then  is  to  select  the  decision  stream 
Dg  that  optimises  the  performance  (i.e.,  utility)  function: 

P  =  YCX,Wt),D(l)V‘t>[x{m)\ 

Assuming  that  we  are  able  to  find  a  reasonable  representation 
for  /,(■)  and  o[x(m)],  finding  the  optimal  decision  stream  is 
then  an  optimal  control  problem  of  the  form  discussed  in 
Chapter  1  and  Appendix  1.  In  this  case,  the  decision  variables 
D  (f)  at  each  cycle  are  the  control  variables  and  the  X  (f)  are  the 
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State  variables.  In  Chapter  1  and  Appendix  1,  we  discussed 
under  what  conditions  we  might  expect  such  a  problem  to 
have  a  unique  solution.  For  further  discussion  of  such  solu¬ 
tions,  see  also  Bryson  and  H  o  [14]. 

The  Two-Sided  Game 

The  problem  with  a  two-sided  game  is  the  development  of  a 
second  performance  function  for  the  opposing  side.  This  prob¬ 
lem  can  easily  be  solved  if  we  design  P  such  that  if  Blue 
chooses  to  maximise  P,  then  Red  will  choose  to  minimise  P.  I  n 
this  formulation,  we  assume  two  decision  streams,  one  for  the 

friendly  commander  Bg  =  {6(0),6(1) . 6(m-l)}  and  one  for 

the  enemy  commander  fig  =  {fi (0),fi (1) . fi(m-l)}.  Our 

assumption  implies  that  the  game  is  zero-sum,  that  is,  any 
increase  in  P  for  the  friendly  force  results  in  an  equal  decrease 
for  the  enemy  force  and  vice  versa.  Consequently,  we  wish  to 

maximise  minimise .  , 

(p) 

Bit)e  Rit)e 

subject  to  the  transition  constraint: 

x(t  +  l)=B,  [x(0,  B(t),  0  <  ^  <  m  - 1  . 

Solutionsto  problemsofthistypearefairly  complex  for  all  but 
very  trivial  examples.  For  example,  see  Hillestad  [15]  and 
Berkovitz  and  Dresher  [16].  However,  the  objective  here  is 
not  to  solve  the  wargame  using  one  of  these  techniques,  but 
rather  to  use  the  two-sided  memoryless  multistage  optimal 
control  formulation  as  a  convenient  construct  for  a  formal 
statement  of  the  problem. 
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Gameswith  Equivalent  Decision  Streams 

Given  the  complexity  of  the  two-sided  optimal  control  con¬ 
struct  for  the  wargames,  we  forego  attempts  to  apply  any 
closed  form  solution  and  rely  instead  on  the  replication  of 
instances  of  the  game  from  each  of  the  scenarios  and  for  the 
several  cases  to  be  examined.  Even  this  however  can  be 
extremely  time  consuming  and  therefore  we  wish  to  examine 
those  alternative  cases  in  which  the  decision  streams  are  essen¬ 
tially  different.  That  is,  if  the  introduction  of  alternative 
weapon  systems  in  a  game  does  not  significantly  alter  the  deci¬ 
sion  stream,  Dq,  then  the  two  games  are  considered 
equivalent.  In  this  way,  the  results  of  one  game  can  be  rerun 
under  the  conditions  of  the  second.  Although  the  results  may 
vary  with  each,  the  decision  stream  is  taken  to  be  constant. 

A  simple  example  will  serve  to  illustrate  the  process.  Consider 
a  conflict  in  which  Red  and  Blue  commanders  have  only  artil¬ 
lery  and  tanks  with  which  to  conduct  a  two-cycle  campaign, 
and  we  focus  on  the  use  of  artillery  resources.  0  ur  state  vector 
then  is  x{t)=[xXt\xMyXt\y2{t)y  where: 

Xi(r),ji(r)=the  number  of  BLUE  and  RED  artillery 
pieces  respectively,  and 

X2(r),j2(^)=the  number  of  BLUE  and  RED  tanks 
respectively. 

In  both  cycles,  we  assume  that  the  decision  on  both  sides  is  the 
allocation  of  artillery  fires.  Consequently,  5(f)=[/)j(f),/)2(r)]^  is 
the  Blue  commander's  decision  at  cycle  t  (f  =  0,l)  and 
(f)  =  [/'i(f)/2(f)r  is  the  Red  commander's  decision,  where 
is  the  fraction  of  Blue/  Red  artillery  allocated  to  attack 
Red/  Blue  artillery,  b2it}  =  l-bj(t}  is  the  fraction  of  Blue  artil¬ 
lery  to  be  allocated  against  Red  tanks,  and  r2(f)  =  l-ri(f)  is  the 
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fraction  of  Red  artillery  to  deallocated  against  Blue  tanks.  For 
simplicity,  we  restrict  the  domain  to  0,  .5,  and  1.  The  two- 
cycle  game  thus  described  isillustrated  in  Figure4.8: 


F / gure 4.8:  An  ExampleTwo-CydeGame 

The  transition  function,  ,  issimply  the  combat  adjudication 
model.  If  we  assume  a  simple  Lanchester  differential  model 
and  let  /  =  1  for  Blue  and  i  =  2  for  Red, 

6',  [x{t),  B{t),  i?(r)]  =  {t)y,  {t),  {t)y^  {t),  (r)x,  {t),  {t)x,  (?)]  ^ , 

where: 

0<cr,.i  <1  is  the  effect  of  Blue/  Red  artillery  against  Red/  Blue 
artillery;  and 

0<cr.2  <1  is  the  effect  of  Blue/  Red  artillery  against  Red/  Blue 
tanks. 

The  a^jS  can  be  thought  of  as  single  shot  kill  probabilities 
(SSK  Ps)  and  bj(t}Xj(t}  and  rj(f)yi(f)  represent  the  number  of  Blue/ 
Red  artillery  allocated  to  Red/  Blue  artillery  and  tanks.  T  here- 
fore  the  transition  equations  become: 

Xi(r  +  l)  =  Xi(r)-a2/i(^k(0 
X2  (r  + 1)  =  X2  {t)-a^2r^  {t)y,  [t] 
y,(t  +  \)  =  y,(t)-a,,bXt)x,{t) 

(^  + 1)  =  (^)-  «i2^2  (^) 
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We  must  also  have  that  a^jr.{t)y^{t)<Xj{t)  and  similarly, 
a^jbj{t)xXt)^yj{t)  ■  I  n  other  words,  the  number  of  kills  cannot 
exceed  tne  number  of  target  weapon  systems  available. 

What  remains  to  be  defined  is  a  utility  function  that  is  some 
measure  of  how  well  both  sides  accomplish  their  mission.  For 
this  simple  problem,  we  assume  that  both  sides  wish  to  maxi¬ 
mise  the  number  of  tanks  available  at  the  end  of  the  second 
cycle.  T  heir  reasoning  is  that  as  the  opposing  sides  come  into 
closer  contact,  tanks  are  more  effective  than  artillery.  A  utility 
function  that  does  this  is: 


^  =  ZLh(0->^2(0]  +  0-9[x2(2)-j2(2)]  +  0.l[xi(2)->;i(2)]. 

The  Blue  commander  therefore  wishes  to  maximise  P  and  the 
Red  commander  wishes  to  minimise  P.  Note  that  the  weights, 
0.9  and  0.1,  reflect  the  relative  importance  assigned  to  tanks 
and  artillery  by  the  two  commanders.  0  f  course,  this  objective 
function  is  not  unique.  T  here  are  several  other  possibilities. 

An  interesting  problem  arises  when  the  opposing  commanders 
do  not  agree  on  a  common  objective.  T  hat  is,  what  happens 
when  the  game  is  not  zero-sum?  I  n  this  case,  we  would  need  to 
evaluate  a  separate  objective  function  for  each  side.  In  the 
example  above,  this  would  correspond  to  having  a  utility 
P(Blue)  and  a  utility  P (Red)  that  might  correspond  (for  exam¬ 
ple)  to  different  weightings  for  the  value  of  tanks  versus 
artillery  due  to  different  perceptions  of  the  endstate  (Red  may 
simply  wish  to  survive  with  a  roughly  balanced  force,  for 
example).  We  would  then  require  some  higher  level  measure 
of  what  this  set  of  outcomes  implies.  In  some  cases  (such  as 
0  perations  0  ther  T  han  War),  we  might  be  seeking  to  maxi¬ 
mise  the  utility  of  both  sides  (i.e.,  a  win-win  situation  rather 
than  the  win-lose  assumption  we  made  above). 
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For  the  example  we  have  considered  so  far,  we  can  now  anal¬ 
yse  an  equivalent  decision  streams  case.  The  allocation  of 
artillery  to  enemy  tanks  and  artillery  over  the  two  cycles  is 
referred  to  as  the  allocation  strategy.  Figure  4.9  lists  all  the 
possible  strategies  by  describing  the  game  for  the  Blue  com¬ 
mander  in  extensive  form.  Note  that  the  diagram  reflects  only 
the  allocation  of  Blue  artillery  to  Red  tanks,  b2it}. 


Figure4.9:  BLUE  Commander's  Allocation  Strategy 

The  Blue  commander  reasons  that  during  the  first  cycle, 
because  the  opposing  forces  are  not  in  direct  contact,  the 
greatest  threat  to  his  forces  is  the  enemy  artillery.  Therefore, 
he  allocates  all  of  his  artillery  against  Red  artillery.  I  n  the  sec¬ 
ond  cycle,  as  the  forces  begin  to  close,  he  sees  Red  tanks  as  the 
more  serious  threat  and  therefore  allocates  all  of  his  resources 
against  Red  tanks.  Therefore,  his  decision  stream  (allocation 
strategy)  is: 


5<P={S(0),S(1)}  =  {[1,0],[0,1]}. 
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This  is  refiected  in  Figure  4.9  by  the  boid  path,  and  corre¬ 
sponds  to  the  beng-beng  soiution  discussed  in  Appendix  1, 
where  oniy  extreme  settings  of  the  controi  variabie  are  used. 
(This  soiution  wouid  resuit  from  a  Lanchester  square  iaw 
Biue/  Red  interaction  for  exampie,  which  wouid  give  rise  to  a 
iinear  H  amiitonian,  asexpiained  in  Appendix  1.) 

Now,  assume  that  a  new  game  is  to  be  piayed  that  differs  from 
the  current  game  oniy  in  that  the  Biue  force  has  been  aug¬ 
mented  by  more  tanks,  if  the  addition  of  these  tanks  does  not 
aiter  the  decision  stream,  then  we  consider  them  equivaient 
and  the  game  may  be  repiayed  with  the  same  decisions  on 
both  sides  and  the  outcomes  (vaiueofP)  compared. 

DECiSiON  U  NCERTAiNTY 

ASTOR's  primary  function  is  to  contribute  to  tacticai  situa¬ 
tion  assessment  by  observing  the  battiefieid,  detecting  and 
identifying  enemy  units,  and  reporting  on  its  findings.  Conse- 
quentiy,  a  metric  designed  to  measure  how  weii  situation 
assessment  isaccompiished  in  aii  cases  tested  was  seen  asuse- 
fui  to  this  study  of  the  ASTO  R  sensor  system.  Such  a  metric 
aiiows  us  to  measure  the  degree  of  confidence  the  commender  has  thet 
he  possesses  an  eccurete  picture  of  the  batti^ieid  in  his  area  of  interest. 
We  wouid  expect  that  the  greater  his  knowiedge  about  the 
iocation,  size,  and  composition  of  the  enemy  force,  the  greater 
his  confidence  in  making  decisions  concerning  the  aiiocation 
of  his  weapons  and  the  movement  of  his  forces.  Weaiso  recog¬ 
nise  that  information  of  this  type  is  not  aii  he  wouid  require, 
information  concerning  enemy  intent  gieaned  from  COM- 
iNT,  SiGiNT,  and  known  enemy  fighting  doctrine  wouid  aiso 
assist  in  compieting  the  picture. 
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We  developed  such  a  metric  that  reflects  the  amount  of  knowl¬ 
edge  the  commander  has  concerning  the  enemy  forces  arrayed 
against  him  in  his  area  of  interest.  The  measure  is  a  function  of 
the  size,  diversity,  and  effectiveness  of  the  sensor  suite,  and  the 
effectivenessof  the  command  and  control  system  used  to  pro¬ 
cess  the  reported  sensor  observations.  The  detailed  develop¬ 
ment  of  this  metric  is  covered  next. 

Probability  Distribution 

We  begin  by  letting  the  vector  U  represent  the  competing 
hypotheses  that  any  number  of  enemy  units  are  arrayed 
against  the  friendly  commander  at  time  cycle  t  so  that 

U  ={0,1,2 . n}.  Given  the  level  of  resolution  for  the  ASTOR 

games,  a  unit  was  taken  to  be  a  battalion.  We  omit  the  cycle 
index,  f,  for  now  focusing  instead  on  analysiswithin  a  timestep. 
T  he  term  arrayed  aga/nsf  is  taken  to  indicate  the  units  located  on 
the  battlefield  in  some  area  of  interest  to  the  friendly  com¬ 
mander.  T  his  may  mean  along  some  avenue  of  approach  in  a 
defensive  operation  or  blocking  a  route  of  advance  in  an  offen¬ 
sive  operation.  Figure  4.10  depicts  a  notional  defensive 
campaign  situation. 

We  assume  that  the  friendly  commander  knows  the  number  of 
enemy  units  that  might  be  brought  to  bear  against  him  during 
the  campaign.  T  hat  is,  we  assume  that  he  knows n.  T  his  is  a  rea¬ 
sonable  assumption  in  that  it  is  highly  likely  that  the  I  ntelligence 
Preparation  of  the  Battlefield  (IPB)  would  yield  this  informa¬ 
tion.  What  is  unknown  is  the  tactical  deployment  of  the  units  at 
each  timestep.  Tactical  situation  assessment  then  is  taken  to  be 
the  process  of  estimating  the  enemy's  tactical  deployment  at 
time  t  and  the  effectiveness  of  this  estimate  is  the  degree  of 
uncertainty  associated  with  his  current  state  of  knowledge. 


T  he  Decision  Probiem 


110  Complexity  Theory  and  Network  Centric  Warfare 


Friendly 

Defensive 

Positions 


u\+U2+Ui—n 


Figure4.10:  BLUE  C ommdnder's  Situdtion  Assessment  Problem 


Bayesian  Decisionmaking 

We  begin  by  analysing  the  intelligence  gathering  process  at 
each  timestep.  We  first  assume  that  a  Bayesian  update  meth¬ 
odology  for  tactical  situation  assessment  is  appropriate  within 
a  wargame  cycle,  but  not  between  wargame  cycles,  given  the 
assumptions  concerning  the  M  arkov  properties  (i.e.,  lack  of 
memory)  of  the  FASTHEX  game  with  2-hour  timesteps.'^ 
Consequently,  the  process  described  here  is  repeated  prior  to 
each  decision  to  commit  forces. 

1.  Input  Distribution :Thefriendly  commander  may 
or  may  not  have  some  idea  of  the  likely  disposition  of 
enemy  units.  If  he  does,  we  may  describe  it  using  an 
empirical  distribution.  H  owever,  for  this  analysis,  we 
assume  that  the  friendly  commander  is  completely  igno- 


%e  later  exploit  Bayesian  updating  by  assuming  multiple  sensor  sweeps 
within  a  single  decision  cycle. 
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rant  of  the  enemy  commander's  intentions.  T  his 
provides  us  with  a  worst  case  situation,  corresponding  to 
the  assumed  lack  of  memory  between  timesteps.  We  let 
P(U  =  u)  represent  the  probability  that  the  enemy  com¬ 
mander  will  commit  u  of  his n  units  in  a  specified  area  of 
interest  in  the  area  of  operations  AO  (avenue  of 
approach  in  Figure  4.10).  Assuming  that  the  enemy 
commander  is  equally  likely  to  deploy  any  number  of 
units  in  the  area  of  interest,  we  have  that 
PiU  =u}  =  l/  n+1.  T  he  friendly  commander  hopes  to 
refine  this  distribution  using  his  sensor  assets. 

2.  The  Sensor  Model:  We  next  let  1/  =  {0,1,2 . n}  rep¬ 

resent  the  number  of  units  detected  by  the  sensor  assets 
allocated  to  the  area  of  interest.^  T  herefore,  P  (1/  =  i/)  is 
the  probability  that  the  sensors  will  detect  i/ of  the  enemy 
units  arrayed  against  the  friendly  forces.  H  owever,  this 
number  is  conditioned  on  the  number  of  units  in  the 
area  of  interest  deployed  by  the  enemy  commander. 
Consequently,  we  focus  on  the  conditional  probability, 
p(v  =  V I  c/  =  m)  .  For  simplicity,  we  assume  a  single  sen¬ 
sor  is  cued  to  search  the  specified  area  of  interest.®  We 
further  assume  that  the  sensor  is  capable  of  detecting  a 
unit  in  the  area  of  interest  with  probability  q  and  that 
there  are  no  false  detections  from  the  sensor  or  else¬ 
where.^  Consequently,  the  conditional  probability 
distribution,  F(r  =  viu  =  u),  isthe binomial 
distribution: 


^By  "detect"  we  mean  that  sufficient  information  is  provided  to  allow  the 
unit  to  be  targeted  by  a  weapon. 

his  assumption  can  be  relaxed  to  allow  for  the  characterisation  of  a 
multisensor  suite,  provided  that  the  sensors  are  independent. 
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p{V  =  v\U  =  u)  =  \ 


for  v<  M 

0  otherwise. 


v^y 


3.  Sensor  0  perations:  0  ur  objective  isto  clarify  the 
enemy  force  deployment  picture  based  on  the  sensor 
observations  by  refining  the  friendly  commander's  ini¬ 
tial  and  subsequent  probability  distributions  on  d.That 
is,  we  wish  to  calculate  p{u  =  u\  v  =  v^) ,  where  i/^  is  the 
number  of  detections  reported  in  the  cycle,  and  thus 
assess  the  i  mpact  of  the  evidence  provided  by  the  sensor 
on  our  estimate  of  the  number  of  enemy  units  arrayed 
against  the  friendly  forces  in  the  area  of  interest.  Opera¬ 
tionally,  we  assume  that  the  sensor  sweeps  the  area  of 
interest  once  in  a  cycle.  Asa  detection  occurs,  it  is 
immediately  reported  so  that  there  are  i/d-i-1  reports 
from  the  sensor  per  cycle.  T  he  additional  report 
accounts  for  the  fact  that  a  report  of  0  detections  is  sent 
initially.  Since  it  is  impossible  to  control  the  time  when 
detections  occur  within  a  FAST  H  EX  game  cycle,  we 
assume  a  uniform  distribution  of  reports.  T  hat  is,  a 
report  of  no  detections  occurs  at  time  f/  (i/d+1),  a  report 
of  one  occurs  at  2f/  (1/^+1),  etc.  The  estimate  is  refined  at 
every  subinterval  using  Bayes'  formula  as  follows: 


p{U  =  u\V  =  v) 


p{u  =  u\  v  =  v-\)p{y  =  v\U  =  u) 

Y,‘_j,P(U=i\V  =  v-\)P(V  =  v\U=i) 


^We  later  relax  this  assumption  by  allowing  for  the  possibility  that  the 
sensor  detections/  identifications  are  false,  that  the  command  and  control 
system  used  to  transmit  the  sensor  information  may  report  a  false 
detection/  identification  as  real,  and  that  the  intelligence  processing  centre 
may  interpret  a  false  detection/  identification  as  real. 
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In  this  formulation,  p{u  =  u\  v  =  v-\)  isthe  prior  prob¬ 
ability,  p{y  =  v\u  =  u)  isthe  knowledge  contributed  by 
the  latest  report  (the  probability  that  one  more  unit  is 
detected),  and  p{u  =  u\  v  =  v)  isthe  posterior  probabil¬ 
ity  on  U  given  the  last  report.  N  ote  that 
p(u  =  u\v  =  -I}=p(u  =  u}=i/(n+i};  that  is,  the  prior 
distribution  before  sensors  a  re  deployed  is  flat,  as 
described  above.  T  his  process  is  repeated  for 

1/  =  0,1 . 1/^.  M  aking  the  appropriate  substitutions  in  (2), 

we  get: 


P(U  =  U  IF  =  v-l]l 


p(U  =  u\F  =  v}-- 


x;y/’(t/=;iF=v-i)  5I1-9)' 


(3) 


P{U  =  u\V  =  v-\\  {\-q) 


v^y 


Y,^P(U  =  i\V  =  v-i)  (l-q)‘ 


V’'/ 


where  1/  =  0,1 . i/^  isthe  number  of  units  detected  by  the 

sensor  and  u>v  at  each  iteration.  Figure  4.11  depictsthe 
process  diagrammatically.  Note  the  difference  between  no 
sensor  sweep  in  progress  and  a  report  of  no  detections 
Theformer  isdepicted  by  a  flat  probability  distribution  on 
U  whereas  the  latter  is  a  refinement  to  the  flat  distribution. 
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Initial  Estimate 


(t/=«)={0,l . n]  j/^(t/  =  H)=^. 


(U=ulV=Vj)={v^  Vrf+l, 


P{U  =  u\V  =  v,)=- 


'  ’  d‘ 

Refined  Estimate 

F igure 4.11:  D  e/d oping  a  R  efined  E stimdte 


A  Simple  Example 

Thefollowing  illustratesthe  process.  Table  4.1  summarises  the 
results  of  a  simple  situation  in  which  three  units  are  known  to 
be  available  to  the  enemy  commander.  T  he  sensor  system  has 
a  probability  of  detection/  identification  of  q  =  0.8.  The  entries 
in  the  rows  are  the  refined  probabilities  from  0,  1,  2,  and  3 
detections.  T  he  first  row  is  the  a  priori  probability  assessment 
on  U  assuming  initial  total  ignorance.  Figure  4.12  depicts  the 
results  graphically. 
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V 

f(u  =  oiv} 

F(U  =  11  v) 

F(U  =  2IV} 

F(U  =  3IV} 

- 

.0321 

1 

.0736 

.9634 

3 

0 

1 

Table  4.1:  Refined  Probability  Assessments:  Example  1 


nr=«iry 


Figure4.12:  R  efined  Probability  Assessments 

Multiple  S\a/eeps 

We  now  refine  the  analysis  to  show  the  value  that  multiple  sen¬ 
sor  sweeps  within  the  same  cycle  have  on  refining  the 
probability  estimates  for  U.  Suppose  that  we  assume  that  our 
sensors  are  capable  of  k  sweeps  of  the  area  of  interest  within 
the  commander's  decision  cycle.  That  is,  the  sensor  can  per¬ 
form  k  sweeps  of  the  area  of  interest  before  the  enemy 
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commander  can  move  his  units  in  any  significant  way.  i  n  each 
of  these  sweeps  (/),  1/^,  enemy  units  are  detected.  We  further 
assume  that  the  probabiiity  estimates  are  made  sequentiaiiy, 
and  that  the  sweep  time  issufficientiy  smaii  to  aiiow  for  a  sin- 
gie"end  of  sweep”  report.  Using  B  ayes  for  mu  i  a  we  get: 


p{U=u\V=v,,) 


P(l/=«|F=v,„.„)p(K  =  v„|t/=n) 

Tr-Av  =>!’'  =  =  vJ  t/  =  j) 


(4) 


where  /  =  1,2 . k,  and  F(u  =  ujF  =  v^J  =  F(u  =  u}=l/n+l 

M  aking  the  appropriate  substitutions,  (4)  becomes: 


p(u=„IF  =  vJ= - 

in  generai,  Bayesian  updating  has  a  tendency  to  converge 
rather  rapidiy-especiaiiy  in  cases  such  as  this  where  faise 
detections/  identifications  are  not  aiiowed:  that  is,  it  is  impossi- 
bie  to  overstate  the  number  of  units  actuaiiy  present.  The 
effect  is  that  subsequent  detections  that  report  fewer  units  than 
the  previous  are  totaiiy  ignored.  To  iiiustrate,  consider  a  sim- 
pie  case  in  which  n  =  3  units.  We  assume  that  three  sweeps 
were  conducted  resuiting  in  three  sequent! a i  detections  using  a 
sensor  with  probabiiity  of  detection:  q  =  0.8.  Tabie  4.2  sum¬ 
marises  the  resuits  of  appiying  equation  (5)  with  k  =  3.  The 
number  of  units  detected  each  time  is  iisted  in  the  tabie.  T  he 
number  of  units  in  the  area  of  interest  is  actuaiiy  three  and 
subsequent  observations  that  two  units  were  detected/  identi¬ 
fied  arecompieteiy  ignored. 
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i 

o 

II 

b 

II 

b 

b 

(N 

II 

b 

b 

F(U  =  31  v) 

- 

1 

2 

0 

0 

0.625 

0.375 

0 

0 

0 

1 

3 

2 

0 

0 

0 

1 

T able  4.2:  M  ultiple  Sweeps  C  ase  1 


Now  consider  a  second  case  with  a  somewhat  different  history 
as  depicted  in  T able  4.3.  In  this  case,  four  sweeps  were  con¬ 
ducted  resulting  in  the  sequential  detections  depicted  in  the 
Table.  The  detection  of  one  unit  persisted  for  three  reports. 
Note  the  rapid  convergence  of  p{u  =  \\v).  However,  the  sin¬ 
gle  detection  of  two  units  in  Sweep  4  shifts  the  mode  of  the 
distribution  to  U  =2.  Because  we  exclude  false  detections,  all 
reports  less  than  the  current  number  detected  will  be  ignored. 


/ 

b 

o 

II 

b 

II 

b 

b 

F(U  =  21  v) 

F(U  =  31  v) 

l!l 

Bi 

0:250 

0.250“ 

“0.250 

0 

0.658 

0.263 

0 

0.767 

0.122 

0 

0.925 

0.059 

4 

2 

0 

0 

0.855 

T able  4.3:  M  ultiple  Sweeps  C  ase  2 


False  Target  Detections/ 

Identifications 

U  p  to  this  point,  we  have  assumed  that  false  targets  were  not 
present.  It  is  possible  to  relax  this  assumption  and  recognise 
that  targets  can  be  misclassified  in  several  ways:  (1)  the  sensor 
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system  may  be  defective;  (2)  the  command  and  control  system 
used  to  transmit  the  detection  to  a  central  processing  centre 
may  have  erroneously  introduced  a  false  target;  (3)  the  pro¬ 
cessing  centre  equipment  or  personnel  may  have  misinter¬ 
preted  the  data  being  received;  and  (4)  the  enemy  may  be 
actively  engaging  in  deception  activities  (i.e.,  Information 
Operations).  The  way  in  which  this  is  done  is  described  in 
detail  in  reference  [11],  using  a  Poisson  process  to  represent 
the  "flow”  of  false  targets  through  the  sensing  process. 

Knowledge  Representation 

It  now  remains  to  ascertain  the  degree  of  uncertainty  existing 
in  the  mind  of  the  friendly  commander  at  the  time  he  must 
take  a  decision  on  the  employment  of  his  forces.  H  is  current 
knowledge  consists  of  two  components:  (1)  the  fact  that  his  sen¬ 
sor  suite  detected  a  number  of  enemy  units  in  his  area  of  inter¬ 
est;  and  (2)  the  refined  probability  distribution  over  the 
possible  number  of  enemy  units  that  might  be  in  his  area  of 
interest  based  on  his  most  recent  sensor  report.  The  value  of 
the  first  component  depends  upon  whether  false  detections  are 
possible.  The  second  depends  upon  the  number  of  enemy 
units  detected  and  the  reliability  of  the  sensor  system.  The  task 
isto  develop  a  knowledge  metric  that  incorporates  these  two  com¬ 
ponents,  thereby  quantifying  the  likelihood  that  the  com¬ 
mander  has  a  true  picture  of  the  number  of  enemy  units 
arrayed  against  him  in  his  area  of  interest. 

Information  Entropy 

We  draw  on  information  science  to  develop  a  knowledge  met¬ 
ric  that  is  a  function  of  the  average  information  present  in  the 
set  of  all  possible  uncertain  events.  T  his  quantity  is  referred  to 
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as  informetion  entropy^  and  it  measures  the  emount  of  uncerteinty  in 
a  probability  distribution. 

T  he  amount  of  information  available  from  the  known  occur¬ 
rence  of  the  event,  U  =u,  i.e.  that  u  enemy  units  are  indeed 
arrayed  against  the  friendly  force,  is  inversely  proportional  to 
the  likelihood  that  the  event  will  occur.  An  event  that  is  very 
likely  to  occur  provides  little  information  when  it  does  occur. 
On  the  other  hand,  an  unlikely  event  provides  considerable 
information  when  it  occurs.  M  athematically,  we  define  infor¬ 
mation  as  follows: 

/(U  =  u}  =  ln  =  -lnF(U  =  u).  ^ 

If  we  now  consider  all  of  the  events  in  the  refined  set 
c/|F  =  v^,  we  reason  that  each  occurs  with  probability 
p(u  =  uiv  =  vj.  Therefore,  the  information  available  from 
the  occurrence  of  each  event  is: 


/(l/  =  ulF  =  vJ  =  -lnF(U  =  M  I  F  =  vj  , 

and  the  expected  information  from  the  occurrence  of  each 
event  is: 


M  I  r  =  vj/(c/ =  M  I  F  =  vj  = -F(r/ =  M  I  F  =  vJlnF(C/ =  M  I  F  =  vj. 


he  term  entropy  is  used  because  the  information  entropy  function  isthe 
same  as  that  used  in  statistical  mechanicsfor  the  thermodynamic  quantity 
entropy.  For  a  more  complete  discussion  of  entropy,  see  Blahut  [17]  and 
Zurek,  ed.  [18]. 

^1  n  communication  theory  the  units  of  measurement  are  "bits"  if  base  2 
logarithms  are  used  and  "nits"  if  natural  logarithms  are  used  (See 
Kullback  [19]  p.  7).  For  our  purposes  we  will  assume  a  dimensionless 
quantity. 
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Consequently,  the  average  amount  of  information  in  the  prob¬ 
ability  distribution:  p{u\v  =  v^)  can  be  expressed  as: 

h[p{U\V  =  v,)]  =  H{U\V  =  v,)  =  P{U  =  i\V  =  v,)\n[p{U  =  i\V  =v, )]. 

The  entropy  quantity  h{u\v  =  Vj)  is  the  residudi  uncertdinty 
regarding  U  given  that  1/  is  instantiated  to  i/^.  The  average 
uncertainty  then  is  the  sum  of  the  residual  uncertainties 
weighted  by  the  probability  distribution  on  the  sensor  detec¬ 
tion/  identifications: 


H(uiv)= -X%p(v  =  J)Y,:^P(U = ;  I K  =  y)in[p(iy  = ;  |  f  =  y)J 


Propertiesof  Information  Entropy 

Information  entropy  has  properties  that  make  it  ideal  as  a  met¬ 
ric  for  measuring  the  commander's  uncertainty  prior  to 
making  a  decision  and  for  measuring  the  uncertainty  in  the 
entire  campaign: 

1.  Maximum  Entropy:  The  entropy  function  is  maxi¬ 
mised  when  the  uncertainty  in  the  distribution  is 
greatest.  M  aximum  uncertainty  occurs  when  the 
friendly  commander  has  no  sensor  assets  to  deploy.  I  n 
this  case,  any  number  of  units  might  be  arrayed  against 
him  with  equal  probability.  M  athematically  we  have 
that  p(u  =  u}=l/n+l .  The  entropy  in  this  case  is: 

n  +  l  n  +  l 

Thus  the  maximum  uncertainty  in  P(U)  isln(n-Fl).  Note 
that  asn  grows  larger,  the  entropy  increases  as  we  might 
expect;  the  more  units  available  to  the  enemy  com- 
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mander,  the  less  clear  we  are  about  their  deployment  in 
the  absence  of  sensor  outputs.  In  general,  a  probability 
distribution  with  a  wide  variance  exhibits  high  entropy. 

2.  Minimum  ENTROPY:Theentropyfunction  ismini- 
mised  at  0.  T  his  occurs  when  P{U  =  u,)  =  1.0  and 

P(U  =Uj)  =  0forall  This  represents  total  certainty 
or  minimum  uncertainty. 

3.  Campaign  Entropy:  The  total  campaign  entropy, 
denoted  H  (U  i,U  ml  where  m  is  the  total  number  of 
game  (campaign)  cycles  satisfies  the  relation: 

The  equality  condition  holds  when  the  process  is  memo¬ 
ryless  as  approximated  in  the  FAST  HEX  games  (for 
purposes  of  tactical  situation  assessment),  when  the  situ¬ 
ation  being  considered  is  rapidly  changing  across  the 
timespan  of  the  campaign. 

Combat  Cycle  Knowledge 

For  the  wargames  we  were  dealing  with,  we  found  it  impor¬ 
tant  to  develop  a  metric  that  was  capable  of  providing  an 
ordinal  ranking  of  the  wargame  cases  across  all  scenarios  in 
terms  of  the  knowledge  possessed  by  the  commander  prior  to 
making  a  decision  at  eac/i  cycle.  Although  entropy  is  a  conve¬ 
nient  measure  of  decision  uncertainty,  making  direct 
comparisons  among  the  cases  examined  could  be  misleading. 
We  need  only  recall  that  maximum  entropy  isH  (U}  =  ln(n-Fl) 
to  realise  that  the  varying  number  of  enemy  units  in  the  AO 
makes  a  direct  comparison  incorrect.  I  n  addition,  it  is  incom¬ 
plete  because  it  addresses  only  the  second  knowledge 
component,  namely  the  knowledge  gained  from  the  refined 
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probability  distribution.  What  is  needed  is  a  more  comprehen¬ 
sive  metric  incorporating  the  residual  uncertainty  in  the  refined 
distribution,  and  the  detection  information  gained  by  the  sensor 
report.  H  aving  said  this,  in  later  applications  of  this  idea  to 
quantifying  Information  Dominance  [6]  and  the  benefits  of  a 
network-centric  approach  to  collaboration  and  task  prosecu¬ 
tion  [9],  the  simpler  form  of  Residual  U  ncertainty  (and  hence. 
Residual  K  nowledge-which  is  a  measure  of  uncertainty 
removed)  has  been  found  to  be  adequate. 

For  the  wargaming  application  then,  we  let  =  i/^)  repre¬ 
sent  the  knowledge  gained  from  detecting  i/^  enemy  units  when 
there  are  U  enemy  units  in  the  area  of  interest.  Symbolically 
we  have: 


K(U,V  =  v,}  =  K(U\V  =  vMv  =  v,} 

where  k(u\v  =  vJ  isthe  knowledge  associated  with  the  resid¬ 
ual  uncertainty  in  the  refined  probability  distribution  given  a 
sensor  report  of  i/^  units,  and  /((!/=  1/^)  isthe  knowledge  gained 
by  detecting/  identifying  i/^  enemy  units.  If  we  can  ensure  that 
both  k(u\v  =  vJ  and  KiV  =  i/^)  are  confined  to  the  interval 
[0,1],  then  we  can  think  of  KiU,V  =  i/^)  as  a  probability.^®  As 
such,  it  represents  the  likelihood  that  the  commander  has  a 
complete  picture  of  the  battlefield  at  the  time  he  makes  a  deci¬ 
sion.  This  can  be  a  very  powerful  statistic  when  correlated 
with  the  Force  Loss  Exchange  Ratio,  enemy  attrition,  and 
friendly  survivability  as  discussed  below. 

1.  Residual  KNOWLEDGE:The maximum  uncertainty  in 
p(i/ \v  =  vj  is In(n-Fl).  T herefore,  maximum  certainty 


%([/,!/  =  l/d)satisfiestheprobability  axioms  (see  Stark  and  Woods  [13] 
p.  9  for  instance)  and  therefore  can  bethought  of  as  a  subjective 
probability. 
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can  bedefined  as  in(n+i)-//(c/|K  =  v^).^^ Normalising 
this  quantity  provides  us  with  the  following  definition  of 
residual  knowledge: 

Residual  knowledge  is  maximised  when  residual  entropy 
is  0  and  it  is  minimised  when  residual  entropy  is  ln(n+l). 
In  general,  residual  knowledge  reflects  the  amount  of 
uncertainty  in  the  refined  probability  distribution. 

2.  Detection  Knowledge:  Given  that  i/^  enemy  units 
were  detected,  we  are  now  concerned  with  the  addi¬ 
tional  information  this  provides  concerning  the 
likelihood  that  there  are  actually  i/^  or  more  enemy  units 
in  the  area.  T  his  is  clearly  a  function  of  the  reliability  of 
the  sensors  and  the  command  and  control  system  used 
to  process  the  sensor  data  it  receives  M  athematically, 
we  are  interested  in  the  information  content  for  the 
event:  c/  >  |  K  =  .  T  hat  is,  the  information  that  will 

be  provided  from  the  detection  reports  this  cycle,  or  the 
prior  information  content  of  the  event,  1/  =  This  is  cal¬ 
culated  to  be: 

l{U>vAV  =  v,)  =  -\v{p{U>vAV  =  v,)]  =  -\n^l^P{U  =  i\V  =  v,)\ 

If  1/d  =  0,  we  get  no  information  because  p(/7>o)=i . 
However,  if  i/d  =  n,  the  information  content  is  maxi¬ 
mised  at  -in[p(c/  =  n|K  =  v^)].  This  is  due  to  the  fact 
that  p{u>u\v  =  v^)  decreases  monotonically  with 


^^1  n  general,  the  change  in  information  resulting  from  detecting  1/  =  l/^ 
units  is  M{U  I V  =  v^)=H{u)-H{!J  |  F  =  vJ  . 
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increasing  u  and  therefore  is  smallest  for  u  =  n.  T  his  sug¬ 
gests  the  following  definition  for  /((!/  =  1/^): 

,  ,  -1)1 
ln[P{U  =  n\V  =  Vj-l)] 

(We  use  vl  to  ensure  that  the  denominator  never  goes  to 
zero). 

T  he  total  knowledge  gained  is  then  defined  to  be  the  product 
of  residual  and  detection  knowledge: 


K{U,V  =  v,)  = 


\n{n  +  \)-H{U\V  =  Vj) 
lii(«  +  l) 


In 


P(U  =  i\V  =  v,-l) 


\n[p{U  =  n\V  =  v,-l) 


(7) 


We  can  apply  equation  7  to  the  example  depicted  in  Table 
4.1.  T  he  first  5  columns  of  Table  4.4  repeat  the  information  in 
T able  4.1  for  convenience.  T  he  last  two  columns  contain  the 
entropy  and  knowledge  figures  based  on  the  refined  distribu¬ 
tions  at  each  iteration,  and  the  intermediate  values  of  1/. 
Figure  4.13  depicts  the  results  graphically. 
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T able  4.4:  T otal  K  nowledge:  Example  1 
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F  igure  4.13:  Know  ledge  and  E  ntropy  for  E  xample  1 

Cam paign  K nowledge 

A  similar  formulation  may  now  be  used  to  calculate  campaign 
knowledge,  given  that  the  FAST  HEX  games  are  taken  to  be 
memoryless  processes  for  tactical  situation  assessment.  Con¬ 
sider  a  campaign  consisting  of  m  cycles.  At  each  cycle,  t, 
enemy  units  are  detected  by  the  sensor.  At  each  cycle,  the 
number  of  possible  enemy  units  arrayed  against  the  friendly 
forces,  n,  is  likely  to  be  reduced  as  a  resultof  combat  during  the 
cycle  so  that  is  the  total  number  of  enemy  forces  that  might 
be  arrayed  against  the  friendly  forces  in  the  area  of  interest. 
H(i/^  I  F  =  Vjj  then  represents  the  residual  uncertainty  at  each 
cycle  and  the  total  campaign  uncertainty  is  expressed  as: 

//  (C/,  I  F  =  V,  „  C/,  I  F  =  , . . . ,  (7  J  F  =  v„„ )  =  £"  W(C/,  I  F  =  v,„ ) 

=  -E'.,E:o^W,  =.|F=v„)lnP(C/,  =.|F=v„) 
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By  analogue,  residual  knowledge  for  the  entire  campaign  can 
be  defined  as: 


,  ,  Y'\\n{n,+\)-H{U,\V  =  v,)] 

K{U,\V  =  vMV  =  v,„-,UJV^vJ=^-^ 


Detection  knowledge  can  be  calculated  in  a  similar  way.  In 
this  case,  we  are  interested  in  the  total  information  available 
from  having  detected/  identified  i/^t  enemy  units  at  each  of  the 
m  cycles,  f.  To  do  this,  we  rely  on  the  fact  that  the  total  infor¬ 
mation  available  from  the  occurrence  of  m  independent  events 
is  the  sum  of  the  information  available  from  the  occurrence  of 
each  of  them.  T  herefore  we  get  that: 


I(U,>v„\V  =  -1,C/,  >  V,,  I  F  =  >  V,.  I  F  =  -1) 

Detection  knowledge  for  the  entire  campaign  can  then  be 
expressed  as: 

im  >v„  I  F  =  v„ -l.U,>v„  I  F  =  v,, -l-.U,  >v„.  I  F  =  v^  -1) 

X:,ln[/>(t/,=«,|F  =  v„-l)] 
and  total  campaign  knowledge  is: 

K(U„V=v_„,U,.V  =  v,,,-.U„V=yJ 
[to(«, +!)-//(;/,  I  F=vJ] 

+1)  I  -1)] 
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An  Exam  ple 

C  onsider  the  example  summarised  in  T able  4.5.  T  he  campaign 
consists  of  5  cycles.  At  each  cycle  t,  the  detection  probability  q^, 
the  number  of  units  detected  i/^tand  the  maximum  size  of  the 
enemy  force  is  given.  T  he  last  three  columns  depict  the  residual 
uncertainty  in  the  refined  probability  distribution,  the  informa¬ 
tion  available  from  the  detection  of  i/^t  enemy  units,  and  the 
"probability"  that  the  commander  has  an  accurate  picture  of 
the  number  of  enemy  units  in  his  area  of  interest.  T  he  last  row 
reflects  his  total  campaign  knowledge. 
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T able  4.5:  T otal  K  nowledge 

The  detection  probability  Qf  is  assumed  to  change  as  a  func¬ 
tion  of  time  t  to  reflect  the  changing  sensor  mix.  The 
maximum  size  of  the  enemy  force  fij-  reduces  over  time  to 
reflect  attrition.  The  fluctuation  in  numbers  detected  % 
(including  a  complete  lack  of  detections  during  one  time  cycle) 
leads  to  a  reduced  value  of  overall  campaign  knowledge.  This 
is  scaled  to  vary  between  0  and  1,  with  0  representing  com¬ 
plete  ignorance,  and  1  representing  complete  knowledge  of 
the  number  of  enemy  units  in  the  commander's  area  of  inter¬ 
est  at  every  stage  of  the  campaign. 
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The  Effects  of  K  nowledge 

We  stated  at  the  outset  that  it  was  desirable  to  assess  the  effects 
of  increased  knowledge  on  the  outcome  of  the  campaign.  One 
way  to  do  this  is  to  compare  K{U,V  =  v^)  to  the  Force  Loss 
Exchange  Ratio  (FLER)  and  the  friendly  and  enemy  combat 
attrition.  Comparisons  with  the  FLER  measure  how  knowl¬ 
edge  influences  the  relative  losses  in  combat.  Comparisons 
with  friendly  and  enemy  attrition  measure  the  degree  to  which 
knowledge  enhances  the  survivability  of  friendly  forces  and  the 
destruction  of  the  enemy.  Statistically,  we  have  shown  a  posi¬ 
tive  linear  relationship  when  the  FLER  or  enemy  attrition  is 
compared  to  knowledge,  and  a  negative  linear  relationship 
when  friendly  casualties  are  compared  to  knowledge  [11].  The 
results  also  showed  that  the  approach  and  structure  (the  use  of 
the  epitomising  approach,  and  open  gaming)  adopted  in  the 
wargameswasableto  produce  a  set  of  coherent  and  quantified 
alternatives,  which  formed  the  basis  of  the  balance  of  invest¬ 
ment  analysis.  As  a  consequence  of  this  analysis,  information 
could  be  weighed  in  the  same  scale  as  weapon  effects,  and  the 
benefit  clearly  demonstrated. 

Figure  4.14  shows  the  relation  between  campaign  level  knowl¬ 
edge  and  attrition  of  enemy  forces,  as  assessed  using  the 
FAST  HEX  wargaming  experiments,  giving  a  correlation 
value  of  0.8. 

Figure  4.15  shows  the  effect  of  an  increase  in  campaign  level 
knowledge  on  own  force  casualties,  again  as  a  result  of  the 
FAST  HEX  wargaming  experiments,  giving  a  correlation 
value  of -0.4. 
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As  discussed  in  [11]  from  the  experimental  results,  there 
appears  to  be  a  point  where  the  knowledge  available  to  the 
commander  exceeds  his  capacity  to  act  on  it,  either  to  gener¬ 
ate  more  enemy  losses  or  to  prevent  further  friendly  casualties, 
and  we  need  to  represent  this  effect. 

Quantifying  the  Benefit  of 
Collaboration  Across  AN  Information 
Network 

In  further  exploitation  of  these  ideas.  Perry  [9]  has  shown  how 
this  approach  can  be  used  as  a  basisfor  quantifying  the  benefit 
of  collaborating  across  an  information  network.  An  example 
of  the  approach  is  described  below,  as  used  in  recent  work  by 
DstI  in  the  UK  in  the  context  of  modelling  a  time-critical 
operation.  Full  details  of  the  general  approach  and  other  areas 
of  application  are  in  [9]. 

We  assume  we  have  a  network  of  command  and  control  nodes 
that  are  involved  in  coordinating  a  time-critical  operation. 
Each  of  these  nodes  has  a  number  of  information  pro¬ 
cessing  tasks  to  perform.  If: 


X 

is  the  mean  time  for  node  i  to  complete  all  of  its  tasks,  we 
assume  that  this  completion  time  is  distributed  exponentially 
(an  exponential  distribution  is  used  to  model  the  time  between 
events  or  how  long  it  takes  to  complete  a  task),  so  that  if  /jff)  is 
the  probability  of  completing  all  tasks  at  node  /  by  timet,  then: 

In  general,  there  will  be  a  number  of  parallel  and  sequential 
nodes  in  the  network  sustaining  the  operation.  Let  this  total 
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number  be  r .  I  n  the  simplest  case,  there  is  a  critical  path  con¬ 
sisting  of  p  nodes  where  is  a  subset  of  r,  as  shown  in 
Figure  4.16. 


Figure4.16:  T  heCriticdl  Path 

We  define  the  total  latency  of  the  path  as  the  sum  of  the 
delays  (latencies)  at  each  of  the  nodes,  plus  the  time,  defined 
as  tm,  required  to  move  a  terminal  attack  system  (such  as  an 
aircraft)  to  the  terminal  attack  area.  In  this  sequential  case, 
we  thus  have  that  the  total  expected  latency  T  is  the  sum  of 
the  expected  latencies  at  each  node  on  the  critical  path,  plus 
the  time  fn,: 


If  there  are  sequential  and  parallel  nodes  on  the  critical  path, 
these  can  be  dealt  with  in  the  way  shown  by  the  example 
below: 


Xoclc  3 


Node  2 


Figure4.17:  Parallel  N  odes  on  theCrItIcal  Path 
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In  this  example: 


T  =  max(— ,  — )  +  —  +  —  +  f  ■ 

/t,  /I2  A,  Xt4 

Returning  now  to  the  case  of  a  serial  set  of  nodes  that  consti¬ 
tute  the  critical  path,  for  each  such  node  i  on  the  critical  path 
defi ne  the /ndegreed/ to  be  the  number  of  command  and  control 
(C2)  network  edges  having  /  as  a  terminal  link. 

For  each  node  j  in  the  C2  network,  we  assume  (based  on  our 
earlier  discussion  of  information  entropy  and  knowledge)  that 
the  amount  of  knowledge  available  at  node  j  concerning  its 
ability  to  process  the  information  and  provide  quality  collabo¬ 
ration  is  a  function  of  the  uncertainty  in  the  distribution  of 
information  processing  time  ffi)  at  node  j.  T  hus  the  more  we 
know  about  node  j  processes,  the  better  the  quality  of  collabo¬ 
ration  with  node). 

Let  H jit)  be  the  Shannon  entropy  of  the  function  ffi).  Then 
H  j(t}  is  a  measure  ofthis  (residual)  uncertainty  defined  in  terms 
of  a  lack  of  knowledge.  By  definition  of  the  Shannon  entropy, 
we  have: 


Ilj(t)  =  -jln(Aje  ^^‘dt 

0 

Since  the  differential  of  {xe'‘  -e'‘)  is 
it  follows  that  H j  (t)  =  ln(— ) 

If  Ajmin  corresponds  to  a  minimum  rateof  task  completions  at 
node),  then: 


1 

XjVcim 


Q  usntifying  the  B  enefit  of  C  olldboration  k  cross  an  I  nformation  N  etwork 


Chapter  4  133 


corresponds  to  a  maximum  expected  time  to  complete  all  tasks 
at  node  j.  In  order  to  provide  a  normalised  value  of  the  knowl¬ 
edge  Kjit)  available  at  node  j  in  terms  of  the  Shannon  entropy, 
we  define  this  as: 


K/t)  =  In 


f  \ 

e 

J 


-In 


f  \ 
e 


J 


In 


yA,jmvn  J 


if/l^min  <  X-  <  eX-vnm 


■  0  if  Xj  <  X^mm 

■  1  if  X^  >  cX-mm 


Suppose  now  that  node  /  is  on  the  critical  path,  and  node  j  is 
another  network  node  connected  to  node  /.  Let  qj  represent 
the  quality  of  collaboration  obtained  by  including  node  j.  If 
this  is  high,  reference  [9]  assumes  Kj{t)  will  be  close  to  1.  T  he 
effective  latency  at  node  i  is  thus  assumed  to  be  reduced  by  the 
factor  due  to  the  effect  of  this  high  quality  of  col¬ 

laboration.  T  he  factor  is  assumed  to  be  1  if  j  is  one  of  the 
nodes  directly  involved  in  the  time-critical  operation  (but  not 
on  the  critical  path).  It  is  assumed  to  be  0.5  if  node)  is  one  of 
the  other  network  nodes,  to  reflect  a  lower  level  of  collabora¬ 
tion  with  these  nodes. 


It  is  important  to  note  that  the  actual  latency  may  not  be 
reduced  by  this  collaboration,  but  the  ability  to  use  the  time 
more  wisely  through  collaboration  (to  fill  in  missing  parts  of 
the  operational  picture  that  are  available  from  other 
nodes,  etc.)  has  an  impact  that  can  be  expressed  equivelently  in 
terms  of  latency  reduction.  T  he  use  of  such  time  more  wisely 
implies  a  good  knowledge  of  expected  time  to  complete  tasks 
that  can  provide  such  information.  This  is  similar  to  the  use  of 
entropy  and  knowledge  in  the  FASTHEX  wargames  where 
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increased  knowledge  led  to  better  awareness  of  the  layout  of 
the  enemy  force,  and  hence  to  wiser  use  of  the  commander's 
own  forces.  Such  wiser  use,  due  to  information  superiority, 
could  be  quantified  in  that  case  by  an  equivalent  (linear) 
improvement  in  the  number  of  enemy  units  destroyed,  or  a 
reduction  in  own  force  casualties. 

The  balance  to  be  struck  is  that  between  such  enhanced  col¬ 
laboration  and  the  effects  of  information  overload  due  to 
increasing  network  complexity  (which  we  assess  separately  as  a 
function  of  the  number  of  elements  of  the  network  involved  in 
the  task). 

The  total  (equivalent)  reduction  in  latency  at  node  /  due  to  col¬ 
laboration  with  the  network  nodes  connected  to  node  /  is  then 
given  by: 


y=i 

y=i 

Thus  the  total  effective  latency  along  the  critical  path, 
accounting  for  the  positive  effects  of  collaboration,  isgiven  by: 

P 

T  =y^+t 

/-I  A- 

1=1  A  j=\ 

We  noted  in  the  experimental  data  from  the  ASTOR  study 
discussed  earlier  that  we  need  to  represent  the  effect  of  infor¬ 
mation  sdturdtion.  Reference  [9]  also  includes  a  comp/ex/'fypena/fy 
to  account  for  the  fact  that  taking  account  of  additional  net- 
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work  connectivity  leads  to  such  information  overload  effects. 
This  is  the  negative  effect  of  collaboration.  It  leads  to  an 
increase  in  effective  latency  on  the  critical  path.  Following  [9], 
we  define  C  to  be  the  total  number  of  network  connections 
accessed  by  nodes  on  the  critical  path.  For  each  node  i  on 
the  critical  path,  this  isthe  indegreed/.  Thus: 

c='Z‘‘.- 

i=l 

The  value  of  C  is  then  a  measureof  the  complexity  of  the  net¬ 
work.  We  assume  that  the  complexity  effect  associated  with  a 
particular  value  of  C  follows  a  nonlinear  S-shaped  curve  as 
shown  below. 


Figure4.18:  T  heLogisticsS-Shaped  Curve 

The  equation  used  to  describe  this  effect  is  a  Logistics 
equation: 

a+bC 

g(C)  =  T-^^- 
1  -I-  e 

T  he  penalty  for  information  overload  is  then  defined  as: 

1 

i-g(C)' 
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T  he  total  effective  latency,  taking  account  of  both  the  positive 
and  negative  effects  of  C  2  network  collaboration,  isthen: 


Network-Centric  Benefit 

This  network-enabled  approach  thus  allows  us  to  compute  the 
distribution  of  the  response  time  of  the  system  as  a  function  of 
the  network  assumptions.  As  we  increase  the  collaboration 
throughout  the  network  in  going  from  platform-centric  to  net- 
work-centric  to  futuristic  network-centric  (to  use  the  RAND 
categories  [9]),  so  the  positive  effects  of  enhanced  collabora¬ 
tion  have  to  balance  off  against  the  downside  effects  of 
information  overload  and  increasing  network  complexity. 
Going  back  to  the  discussion  in  Chapter  2  on  the  Conceptual 
Framework  of  Complexity,  we  can  call  this  overall  assessed 
performance  of  the  network  the  plecticity^^  of  the  network, 
since  it  characterises  the  combined  positive  and  negative 
effects  of  network  complexity  and  collaboration. 

Stochastic  Networksand  Network 
Vulnerability 

So  far,  we  have  shown  how  it  is  possible  to  calculate  the  posi¬ 
tive  and  negative  effects  of  network  complexity  and 
collaboration  based  on  the  use  of  Shannon  entropy  as  a  mea¬ 
sure  of  (lack  of)  knowledge.  We  can  extend  this  model 
potentially  in  a  number  of  ways.  The  length  of  the  critical 
path,  for  example,  if  the  network  is  adapting  over  the  course  of 
a  series  of  tasks,  will  be  a  stochastic  variable.  We  would  expect, 

term  proposed  by  Perry  (RAND  Corp.)  -  personal  communication. 
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from  the  theory  we  have  considered  so  far  that  the  size  of  the 
network  (i.e.,  the  number  of  nodes  on  the  critical  path)  should 
be  sampled  from  a  power-law  distribution  of  network  size.  The 
exponent  of  this  power-law  is  then  a  characteristic  measure  of 
the  ability  of  the  nodes  in  the  network  to  form  and  reform 
dynamically  over  time.  Similarly,  we  can  consider  the  indegree 
of  a  node  on  the  network  to  be  a  stochastic  variable.  If  the 
indegree  of  a  node  is  sampled  from  a  power-law  distribution  of 
the  number  of  links,  then  the  network  is  said  to  be  "scale-free” 
[20].  This  corresponds  to  a  network  with  a  small  number  of 
nodes  with  very  rich  connections,  and  many  nodes  with  sparse 
connections.  (The  Internet  is  an  example  of  a  scale-free  net¬ 
work.)  Conversely,  if  the  indegree  of  a  node  is  sampled  from  a 
normal  distribution  of  the  number  of  links,  then  the  network  is 
of  "random”  type.  Characterising  networks  in  this  way  allows 
us  to  investigate  the  vulnerability  of  such  networks  to  attack,  as 
discussed  in  [20]. 
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An  Extended 
Example  of  the 
Dynamics  OF  Local 
Collaboration  and 
Clustering,  and 
Some  Final 
Thoughts^ 


Towards  the  end  of  Chapter  4,  we  discussed 
the  way  in  which  a  particular  network  could 
be  analysed  using  what  we  called  the  plecticity  of  a 
network,  which  includes  both  the  positive  effects 
of  collaboration  and  the  downside  effects  of  infor¬ 
mation  overload.  This  assumes  a  network  of  a 
particular  size  and  configuration,  and  thus  raises 
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the  question  of  what  sizes  and  configurations  of  such  networks 
of  collaboration  are  likely  to  emerge.  I  n  some  models  of  natural 
systems*  we  have  already  seen  that  the  networks  of  interaction 
that  form  (i.e.,  the  dynamics  of  cluster  formation)  can  be  pre¬ 
dicted  ahead  of  time  to  have  a  particular  form.  For  example  the 
Bak-Sneppen  model  of  the  coevolution  of  the  species  within  an 
evolving  ecosystem  (described  in  Chapter  1)  gives  rise  to  clusters 
of  coevolution  that  tend  to  a  power  law  distribution  of  cluster 
size.  C  lusters  of  burning  trees  (forest  fires)  also  show  such  power 
law  effects*  which  turn  out  to  be  similar  to  the  distribution  of 
casualties  in  war  (Chapters  2  and  3).  Clustering  of  force  units  in 
the  ISAAC  "distillation”  model  of  manoeuvre  warfare  produces 
fractal  clustering  (Chapter  4). 

We  wish  to  finish  with  an  extended  example  analysis  of  the 
ISAAC  distillation;  recall  that  this  is  a  simple  agent-based  model 
of  land  warfare*  incorporating  small  rule  sets  that  govern  agent 
decisionmaking,  movement,  and  engagement.  M  ore  detail  is 
availablein  [1].  Figure  5.1  is  a  screenshot  of  the  start  of  atypical 
ISAAC  simulation  run. 
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Figures. 1:  Screenshot  of  the  Start  of  a  Typical  ISAAC  Simulation  Run 
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ISAAC  is  based  on  cellular  automata,  which  uses  simple  local 
rules  to  describe  interactions  between  units.  Different  scenar¬ 
ios  can  be  set  up  by  changing  the  parameters  of  these  rules. 
The  initial  laydown  of  forces  is  carried  out  stochastically,  by 
the  model  within  user-defined  limits,  and  so  different  replica¬ 
tions  of  the  same  basic  scenario  are  possible. 

T  he  particular  ISAAC  scenario  we  use  in  this  chapter  was  sup¬ 
plied  to  us  by  Dr.  Gary  Horne  of  the  U.S.  Marine  Corps 
Warfighting  Lab  and  has  three  phases.  In  the  first  phase,  the 
Red  forces  move  to  meet  the  Blue;  the  second  phase  is  the 
engagement  between  Red  and  Blue  forces.  It  is  in  the  third 
phase  that  either  theRed  forces  take  control  of  the  Blue  flag  (in 
the  top  right  corner)  or  the  Blue  forces  retain  control  of  their 
flag.  This  scenario  is  interesting  as  in  all  stochastic  replications 
but  one  (replication  40),  theRed  forces  are  successful  in  achiev¬ 
ing  their  goal  of  taking  control  of  the  Blue  flag.  The  analysis 
that  follows  explores  the  clustering  and  "swarming"  of  the 
agents  and  the  similarities  and  differences  in  the  replications. 

Clustering  and  Sia/arming 

In  order  to  analyse  the  swarming  dynamics  of  cluster  forma¬ 
tion  and  dissolution  in  ISAAC,  we  need  to  consider  firstly 
what  this  means.  There  are  two  ways  to  define  a  cluster  of 
agents.  The  first,  and  most  usual,  is  to  define  neighbouring 
agents  only  by  those  that  are  north,  south,  east,  or  west  adja¬ 
cent  to  the  agent  in  question,  known  as  nedrest  ndghbour 
clustering.  The  second  (and  although  most  intuitive,  less  used) 
definition  isto  include  all  eight  neighbours  of  the  central  agent 
as  part  of  a  cluster,  as  shown  in  Figure  5.2,  which  is  known  as 
next  nearest  neighbour  clustering. 
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m  H 


Figure 5.2:  N &rest  and  NextN &rest N eighbour  C lustering 

For  a  single  timestep  in  any  run  of  a  cellular  automata-based 
model,  the  number  and  size  of  clusters  of  agents  can  be  deter¬ 
mined  using  a  simple  algorithm-the  H  oshen-Kopelman 
algorithm-for  the  nearest  neighbour  case.  T  his  algorithm  can 
be  modified  for  use  with  eight  neighbours.  For  details  of  the 
algorithm,  see  references  [2,3]. 

Cluster  Distribution 

0  nee  the  cluster  numbers  and  sizes  can  be  determined,  there 
are  a  number  of  ways  to  analyse  the  data.  The  first  that  we 
look  at  is  the  size  of  the  largest  cluster.  This  gives  an  indication 
of  the  ability  of  the  agents  to  cluster  or  the  amount  of  dispersal 
of  the  agents.  For  example,  if  the  largest  cluster  size  is  near  to 
the  total  number  of  agents,  we  know  that  that  is  the  only  clus¬ 
ter.  FI  owever,  if  the  largest  cluster  is  small,  then  we  know  that 
the  agents  are  dispersed  in  many  small  clusters. 

The  following  plots  are  of  the  largest  cluster  size  against  the 
timestep  for  several  different  replications  with  different  ran¬ 
dom  seeds  for  the  same  basic  run  of  the  ISAAC  model.  The 
clustering  algorithm  used  is  that  of  nearest  neighbours  and  dif¬ 
ferent  plots  are  graphed  for  Red  and  Blue  agents.  T  he  agents 
can  be  ordered  by  state-alive,  injured,  or  dead.  T  he  plots  that 
follow  are  for  only  those  agents  that  are  alive.  For  each 
timestep,  we  plot  the  largest  cluster  size  for  Blue  and  the  larg¬ 
est  cluster  size  for  Red.  For  the  first  iteration  of  our  example 
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run  of  ISAAC,  Figure  5.3  shows  the  evolution  of  the  largest 
cluster  size  for  Red  as  a  function  of  simulated  time.  Figure  5.4 
shows  the  same  thing  for  Blue. 


Time 

Figure  5.3:  L  argest  Cluster  Size  as  a  Function  of  Simulated  T  ime 
(First  I  teration,  Red  A  gents) 


Figures. 4:  L  argest  Cluster  Size  as  a  Function  of  Simulated  T  ime 
(First  Iteration,  Blue  Agents) 


In  these  two  plots  (Figures  5.3  and  5.4),  it  is  possible  to  see  a 
smoothly  changing  pattern  in  the  largest  cluster  size  for  the 
Red  forces,  but  not  for  the  Blue.  Looking  at  the  clustering  of 
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Red,  it  is  clear  that  the  dynamic  behaviour  can  be  split  into 
distinct  areas  in  time  by  the  rate  of  change  of  cluster  size- the 
slope  of  the  plot.  T  he  distinct  areas  in  time  of  the  changes  in 
the  slope  of  the  plot  of  largest  cluster  size  of  Red  agents  corre¬ 
spond  to  the  three  phasesof  the  ISAAC  run. 

In  further  replications  of  the  ISAAC  run,  a  similar  pattern 
emerges-each  time  Red  succeeding  in  his  objective.  H  owever, 
in  the  40th  replication.  Red  fails  to  secure  the  Blue  flag. 
Dr.  H  orne  estimates  that  the  Red  force  is  successful  in  excess 
of  100  replications,  with  only  this  one  failure.  Figures  5.5  and 
5.6  show  the  evolution  of  the  largest  cluster  size  for  this 
replication. 

We  can  see  that  there  is  now  no  clear  evidenceof  a  third  phase 
of  operation  in  the  plot  of  the  largest  Red  cluster  size.  I  n  fact, 
the  plot  of  the  largest  Blue  cluster  size  is  now  more  structured 
and  shows  evidence  of  a  third  phase  of  the  force's  operation. 
We  suggest  that  Red's  failure  is  due  to  the  increased  clustering 
ability  of  the  Blue  forces,  thereby  reaching  a  greater  largest 
cluster  size  than  in  other  replications.  Such  behaviour  is  con¬ 
sistent  with  the  mathematical  metamodel  of  ISAAC  discussed 
in  Chapter  4.  This  metamodel  indicates  that  agile  clustering 
and  reclustering  should  lead  to  better  local  force  ratios  and 
hence  improved  ability  to  cause  attrition  to  the  enemy  (locally) 
and  to  thus  move  freely. 

To  gain  some  additional  insight,  let  us  now  consider  the  larg¬ 
est  cluster  size  at  each  timestep  of  the  simulation,  and  plot 
this  as  a  frequency  distribution  of  cluster  size.  We  have  done 
this  for  a  number  of  replications.  Figure  5.7  shows  the  fre¬ 
quency  of  largest  cluster  size  for  Red  agents  for  typical  and 
exceptional  replications. 
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Figures. 5:  Largest  Cluster  Size  as  a  Function  ofTime(40th  Iteration,  Red) 
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Figures. 6:  L  argest  Cluster  Size  as  a  F  unction  ofTime(40th  Iteration,  Blu^ 
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F  igure  S.7:  F  requency  D  istribution  of  the  L  argest  C  luster  Size  for  Red  A  gents 
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We  can  see  from  this  that  Red  is  abie  to  generate  a  wide 
spectrum  of  ciuster  sizes.  Figure  5.8  shows  the  same  piot  for 
the  Biue  agents  with  the  number  of  each  repiication  shown 
on  the  piot. 


Largest  Cluster  Size 

F igure  5.8:  F requency  D  istribution  of  the L  argest  C luster  Size  for  B  lueAgents 

We  can  see  from  Figure  5.8  that  the  spread  of  ciusters  for 
Biue  is  much  smaiier  in  generai.  Fi  owever,  for  the  40th  repii¬ 
cation,  Biue  is  abie  to  generate  a  wider  spread  of  ciusters,  and 
thus  succeed. 
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The  Distribution  of  Cluster  Size 

Let  us  move  on  from  largest  cluster  size  now  and  look  at  all 
the  clusters  formed  over  time  in  the  simulation.  From  our 
mathematical  metamodel  of  ISAAC  discussed  in  Chapter  4, 
and  from  the  general  emergent  behaviour  of  natural  systems 
we  have  discussed  in  this  book,  we  anticipate  that  the  distribu¬ 
tion  of  cluster  size  should  approximate  to  a  power  law 
distribution.  T  huson  a  log-log  scale,  the  distribution  of  cluster 
size  should  be  a  straight  line,  with  end  effects  where  the 
assumptions  break  down. 

Firstly  let  us  just  look  at  one  replication  of  the  simulation.  Fig¬ 
ure  5.9  shows  the  distribution  of  cluster  size  for  Red  agents  for 
the  2nd  replication  of  I SAAC ,  plotted  on  a  log-log  scale. 


Logarithm  of  cluster  size 

Figures. 9:  Distribution  of  duster  Sizes  (2nd  R  epiicdtion,  R  ed  Agents) 

We  can  see  that  in  the  intermediate  regime,  the  plot  forms  a 
straight  line,  confirming  the  theoretical  expectation.  Figure 
5.10  shows  a  similar  plot  with  the  other  replications 
superimposed. 
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F  igure  5.10:  D  istribution  of  C  luster  Size  for  Red  A  gents 


Finally,  from  theory  and  from  the  analysis  in  Chapter  3,  we 
expect  that  the  time  series  of  casualties  produced  by  a  model 
such  as  ISAAC  should  show  evidence  of  fractal  clustering  in 
time.  In  precise  terms,  this  implies  that  the  power  spectrum  of 
the  time  series  of  casualties  should  be  related  to  the  casualty 
size  as  a  power  law  relationship.  Using  a  related  distillation 
model  called  MANA  [4],  evidence  of  this  effect  has  been 
found  by  Lauren  [5],  as  we  have  already  discussed  in  C  hapter 
3.  T  his  is  an  area  that  we  intend  to  investigate  further  in  the 
context  of  building  metamodelling  equations. 


Final  Thoughts 

We  started  by  considering  what  we  can  learn  from  natural 
systems:  an  ecosystem  in  which  species  coevolve  locally;  a 
fluid  forming  an  interface  when  it  is  pinned;  the  effect  of  forest 
fires.  All  of  these  show  regularities  and  emergent  behaviours 
of  the  whole  system  that  can  be  captured  and  deduced  using 
mathematical  models.  We  have  also  shown  how  the  same 
ideas  of  local  coevolution  within  such  "open”  systems  are  very 
relevant  to  thinking  about  the  consequences  of  a  network-cen¬ 
tric  form  of  warfare,  where  units  coevolve  (self-synch ronise) 
across  an  information  grid.  By  exploiting  this  linkage,  it  ispos- 


Final  T  houghts 


Chapters  149 


sibleto  build  quantitative  models  that  help  us  to  understand 
the  likely  emergent  behaviour  of  such  coevolving  networks  of 
force  interaction. 

This  is  work  in  progress  that  we  hope  will  contribute  to  the 
new  science  of  understanding,  analysing,  and  modelling  the 
effects  of  Information  Age  warfare.  In  doing  so,  we  aim,  as 
remarked  in  the  preface  to  a  previous  contribution,^  to  "gain  a 
deeper  understanding  not  only  of  conflict,  but  also  of  the 
avoidance  of  conflict,  which  is  the  ultimate  aim  of  the  politi¬ 
cal/  military  art.” 
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Optimal  Control 
WITH  A  Unique 
Control  Solution 


In  this  Appendix,  we  investigate  the  case  of  a 
unique  optimai  controi  soiution  to  the  prob- 
iem  of  system  controi,  and  show  that  this  unique 
soiution  is  of  the  form  of  beng-beng  controi  when 
the  system  isiinear  in  nature. 

As  in  C  hapter  1,  we  assume  that  our  system  can 
be  described  by  thefunctionai  reiationship: 


where  f,  are  the  rate  iaws,  and  are  the  con¬ 
troi  variabies.  in  matrix/ vector  notation,  we 
write  this  as: 


x{t)  =  F{X{t\m) 
with  initiai  conditions  xxt,)  =  xl 
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We  assume  fairly  weak  conditions  on  the  continuity  of  f,  X, 
and  X  sufficient  to  make  the  equations  "well  behaved.” 

We  consider  here  analysis  of  this  relationship  as  time  varies 
over  a  fixed  time  interval  [fo,fi].  At  time  fi,  the  state  variables 
will  have  values/, (fi),  and  the  objective  is  to  maximise  or  min¬ 
imise  a  linear  combination  of  these  endstate  values.  The 
problem  can  thus  be  written  as: 

n 

Optimise 

/-I 

(where  the  c,  are  constant  coefficients  or  "weights")  subject  to 
the  constraints: 


V/. 

We  now  define  a  process  [1]  that  yields  a  necessary  condition 
for  a  vector  of  control  variables  X{t)  to  optimise  the  objective 
function.  In  other  words,  any  control  vector  that  gives  rise  to 
an  optimal  value  of  the  objective  function  must  satisfy  this  con¬ 
dition.  Although  this  does  not  guarantee  that  a  solution  X{t) 
satisfying  this  condition  is  optimal,  other  information  (such  as 
the  uniqueness  of  such  a  solution)  can  be  used  in  particular 
cases  to  prove  that  X{t)  is  indeed  an  optimal  control  vector. 

T  he  first  step  in  this  procedure  is  to  introduce  a  set  of  "dual" 
variables: 


¥iitX . 

that  are  defined  by  the  relationships: 


¥Xt)  =  -Y¥jit) 

7=1 


^F.{X,X) 

ax, 
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with  final  conditions: 

¥i(ti)  =  -c,  V/ 

where  the  values  c,  are  the  same  as  the  coefficients  appearing 
in  the  objective  function. 

A  H  amiltonian  function  H  is  now  defined  by: 

X,  A.)  =<ii/,X> 

where  <  >  denotes  the  inner  product  of  the  two  vectors 

y/  and  X . 

T  hus: 


H{y/,X,X)  =  Y,¥iit)XXt) 

/-I 

=  t^¥imiX,A) 

/-I 

From  the  definition  ofH  ,  we  have  the  dual  relationship: 


X:=^  Vz 


dH 


¥i=- 


dx.. 


Vi 


T  he  corresponding  "boundary  conditions”  are: 


X,(^o)  =  X“  Vz 
¥,{t,)  =  -c,  Vi 
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PONTRYAGIN'S  MAXIMUM  PRINCIPLE 

Using  the  notation  we  have  now  developed,  Pontryagin's 
M  aximum  Principle  [2]  states  that  if  JL\t)  is  a  control  vector 
that  maximises  (resp.  minimises)  the  objective  function: 

/—I 

then  the  Hamiltonian  H{y/,x,X)  achieves  a  minimum  (resp. 
maximum)  at  the  point  Xit)  for  any  value  of  X  or  y/. 

It  is  worth  noting  here  that  if  the  set  U  of  admissible  control 
vectors  is  a  topologically  compact  set,  then  the  continuous 
function  H  will  realise  its  minimum  or  maximum  value  on  U. 

T  he  control  vectors  X  that  minimise  or  maximise  the  H  amil- 
tonian  H  are  called  0(tremal  controls.  If  we  denote  this  subset  of 
U  byU*,  then  we  know  from  Pontryagin's  M  aximum  Principle 
that  if  X*  is  an  optimal  control  vector,  then  X*eU\  Thus,  if 
we  know  that: 

1.  An  optimal  control  vector  exists;  and 

2.  There  is  only  one  extremal  control  (i.e.,  U*  is  a  single 
point) 

then  the  single  element  of  U  *  must  be  the  optimal  control  vec¬ 
tor.  Whether  (1)  and  (2)  apply  depends  on  the  particular 
problem  under  study. 

Determining  THE  Extremal  Controls 

T  he  process  for  determining  these  extremal  controls  is  as  fol¬ 
lows  (for  definiteness,  we  assume  that  we  are  maximising  the 
objective  function,  and  hence  minimising  the  H  amiltonian): 


Pontrydgin's  M  aximum  Principle 
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1.  Compute  the  form  of  the  H  amiltonian  function 

2.  Compute  the  values  of  X  which  minimise  H{\i/,x,X) 
for  /le  u . 

T  hese  are  the  extremal  controls.  T  hey  are  expressed  as 
functions  of  X  and  y/ 

X  =  X(X,  y/) . 

3.  Substitute  this  extremal  control  X  into  the  relations: 


dH 


in  order  to  solve  for  the  extremal  system  trajectory  X* 
and  extremal  dual  function  y/ . 

4.  Substitute  these  values  into  the  expression  for  the  extre¬ 
mal  control: 


X  =  X{X\y/) 

to  give  an  explicit  formulation  of  this  extremal  control 
vector. 

Knowledge  of  the  extremal  trajectory  X*  allows  the 
objective  function: 


/=! 


to  be  evaluated. 
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Linear  M odels 

A  linear  system  model  isdefined  as  having  a  relationship  of  the 
form: 


X{t)  =  A{t)X{t)  +  Bm{t)  +  g{t) 

in  a  matrix  and  vector  representation,  with  initial  conditions 

V/. 

When  the  system  behaviour  is  of  this  form,  it  is  possible  to 
characterise  the  nature  of  the  optimal  controls  under  fairly 
general  conditions.  T  his  is  particularly  the  case  if  the  objective 
function  itself  is  linear,  i.e.  it  is  of  the  form: 

h 

\<s,x>  +  <w,x> 

where,  as  before,  <  >  denotes  the  inner  product  of  two  vec¬ 
tors,  and  S  and  1/1/  are  time-dependent  vectors  of  known  value. 

Assume  then  that  the  objective  function  is  of  this  form,  and 
that  the  system  model  is  linear  in  the  way  that  we  have 
described.  Without  loss  of  generality,  we  can  set  g(t}=0  and 
write  the  system  behaviour  model  in  the  form: 

X  =  AX  +  BA 

where  X  is  the  vector  of  state  variables  (e.g.,  force  levels)  and  X 
is  the  vector  of  control  variables. 

M  ake  the  transformation: 

t 

x^,,{t)  =  \<s,x>^<w,X>. 
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The  objective  function  now  becomes: 

Optimise 

and  the  reiationship  x^^,{t)=<S{t),x{t)>  +  <w{t),X{t)>  is 

added  to  the  set  of  equations  describing  the  system  behaviour. 
(This  can  be  done  since  the  above  equation  isiinear  and  so  has 
the  same  form  as  the  others.) 

Consider  now  the  form  of  the  H  amiitonian  for  such  a  iinear 
system.  We  have: 


H{y/,  X,  X)  =<y/,X> 

=<  y/,  AX  +  BX> 

=<  y/,  AX  >  +  <y/,BX> 

Since  we  are  interested  in  the  extremai  controis  X  that  maxi¬ 
mise  or  minimise  H  ,  oniy  the  second  term  is  of  interest,  the 
first  not  being  a  function  of  A . 

Letusiookat  this  second  term  in  moredetaii.  We  have: 

<  y/,  BA  >=  X  ¥i  E  BijAj  =  E  (Z 

i  j  J  i 

Let:  (pj=Y.y^iBy. 

i 

T  hen :  <y/,BA>='^  cp-A.  =<(p,A>. 

J 

Let  us  assume  that  the  objective  function  is  to  be  maximised. 
By  Pontryagin's  M  aximum  Principie  [2],  we  thus  wish  to  con¬ 
sider  controi  vectors  A  that  minimise  the  H  amiitonian  H  .  T  his 
is  then  equivaient,  as  we  have  seen,  to  minimising  the  expres¬ 
sion  Y,(p.{t)A.{t). 
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We  can  define  such  an  extremal  control  vector  /I*  as  follows, 
(provided  that  the  set  U  of  all  possible  control  vectors  is  topo¬ 
logically  compact): 

If  ^j(0  ^  0  let  Aj(t)  =  U} 

and  if  ^j(t)  >  0  let  Aj(t)  =  U} 

T  hese  are  well  defined  since  a  continuous  function  will  attain 
its  max  or  min  on  a  compact  set. 

For  any  fin  the  interval  [fo,  fj  we  then  have: 

J  J 

For  any  control  vector  X  in  the  admissible  set  U  of  control 
vectors. 

Uniqueness  OF  THE  Extremal  Control 
FOR  A  Linear  System 

If  y  is  any  other  extremal  control  vector,  then  since  it  mini¬ 
mises  the  FI  amiltonian  H  ,  it  must  satisfy: 

J  j 

for  any  X  in  the  admissible  set  d  of  control  vectors. 

FI  owever,  if  (p^{t)  >  0 ,  then  it  is  clear  that: 


Vj(t)  =  mm{Xj{t),X^  U)  . 

Otherwise  it  would  be  possible  to  define  a  vector  giving  a 
smaller  value  of  the  FI  amiltonian,  contradicting  the  extremal 
nature  of  1/.  It  follows  that  every  extremal  vector  must  be  of 
the  form  X\ 
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At  the  points  where  (p.{t)  changes  sign,  X]{t)  changes  from 
an  extreme  minimum  value  to  an  extreme  maximum,  or 
vice  versa.  Such  a  form  of  control  is  known  as  beng-beng 
since  the  value  "bangs"  from  one  extreme  possible  value  to 
another,  and  never  assumes  any  intermediate  values.  What 
we  have  shown  is  that  for  a  linear  system  with  a  linear  objec¬ 
tive  function,  every  extremal  control  (including  therefore  the 
optimal  control)  must  be  in  bang- bang  form.  In  this  sense,  the 
optimal  control  vector  always  lies  on  the  boundary  of  the 
admissible  set  d. 

If  it  can  be  shown  that  the  are  unique,  then  the  above  con¬ 
struction  yields  a  unique  control  vector  that  must  then  be  the 
optimal  control.^  Now,  we  have  that: 

(Pj=Y.y'Aj. 

i 

Thus  the  uniqueness  of  (pj  depends  on  the  uniqueness  of 

¥i 

We  have  that: 


7  I  m 

J 

with  final  conditions: 


Wih )  =  -(ci , . ,  )  =  -(0, . ,0,1). 


^Itcan  be  shown  that  for  this  type  of  system,  an  optimal  control  mustexist[l]  [2]. 
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Since  the  matrixA  isknown,  thisequation  has  a  unique  solution 
[1]  and  thus  the  vector  y/  is  unique.  Hence,  the  optimal  control 
for  such  a  linear  system  that  optimises  the  objective  function  is 
precisely  defined  by  the  ibang-ibang  control  function  JL*. 
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Coalition  Connnnand  and  Control* 

(M  aurer,  1994) 

Peace  operations  differ  in  significant  ways  from  tra¬ 
ditional  combat  missions  Asa  result  of  these  unique 
characteristics  command  arrangements  become  far 
more  complex.  T  he  stress  on  command  and  control 
arrangements  and  systems  is  further  exacerbated  by 
the  mission's  increased  political  sensitivity. 


The  M  esh  and  the  N  et 

(Libicki,  1994) 

Considers  the  continuous  revolution  in  information 
technology  as  it  can  be  applied  to  warfare  in  terms 
of  capturing  more  information  (mesh)  and  how  peo¬ 
ple  and  their  machines  can  be  connected  (net). 


Connnnand  Arrangennentsfor 
Peace  Operations 
(AlbertsSi  Hayes,  1995) 

By  almost  any  measure,  the  U  .5.  experience  shows 
that  traditional  C  2  concepts  approaches  and  doc¬ 
trine  are  not  particularly  well  suited  for  peace 
operations  This  book  (1)  explores  the  reasons  for 
this  (2) examinesalternative command  arrangement 
approaches,  and  (3)  describes  the  attributes  of  effec¬ 
tive  command  arrangements 
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■  Standards:  The  Rough  Road  to  the 
Connnnon  Byte 
(Libicki,  1995) 

T  he  inability  of  computers  to  "talk"  to  one  another  is  a 
major  problem,  especially  for  today's  high  technology 
military  forces  T  his  study  by  the  C  enter  for  Advanced 
C  ommand  C  oncepts  and  T echnology  looks  at  the 
growing  but  confusing  body  of  information  technology 
standards.  Among  other  problems  it  discovers  a  persis¬ 
tent  divergence  between  the  perspectives  of  the 
commercial  user  and  those  of  the  government. 


What  Is  Infornnation  Warfare?* 

(Libicki,  1995) 

Is  Information  Warfare  a  nascent,  perhaps  embryonic 
art,  or  simply  the  newest  version  of  a  time-honored  fea¬ 
ture  of  warfare?  I  s  it  a  new  form  of  conflict  that  owes  its 
existence  to  the  burgeoning  global  information  infra¬ 
structure,  or  an  old  one  whose  origin  lies  in  the  wetware 
of  the  human  brain  but  has  been  given  new  life  by  the 
Information  Age?  I  sit  a  unified  field  or  opportunistic 
assemblage? 


INr()KM\fI()N 
WXW  AKI  ^ 


0  perations  0  ther  Than  War* 

(AlbertsSi  Hayes,  1995) 

This  report  documents  the  fourth  in  a  series  of  work¬ 
shops  and  roundtables  organized  by  the  I N  SS  C  enter 
for  Advanced  C  oncepts  and  T echnology  (ACT ).  T  he 
workshop  sought  insights  into  the  process  of  determin¬ 
ing  what  technologies  are  required  for  GOT  W.  The 
group  also  examined  the  complexities  of  introducing 
relevant  technologies  and  discussed  general  and  specific 
OOTW  technologies  and  devices. 
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Donninant  Battlespace  Knowledge* 

(Johnson  &  Libicki,  1996) 

T  he  papers  collected  here  address  the  most  critical 
aspects  of  that  problem—  to  wit:  I  f  the  U  nited  States 
develops  the  means  to  acquire  dominant  battlespace 
knowledge,  how  might  that  affect  the  way  it  goes  to  war, 
the  circumstances  under  which  force  can  and  will  be 
used,  the  purposes  for  its  employment,  and  the  resulting 
alterations  of  the  global  geomilitary  environment? 


Interagency  and  Political-M  ilitary 
Dinnensions  of  Peace  0  perations: 

Haiti  -  A  Case  Study 

(HayesSi  Wheatley,  1996) 

This  report  documents  the  fifth  in  a  seriesof  workshops 
and  roundtables  organized  by  the  I N  SS  C  enter  for 
Advanced  Conceptsand  Technology  (ACT).  Widely 
regarded  as  an  operation  that  "went  right,"  H  aiti 
offered  an  opportunity  to  explore  interagency  relations 
in  an  operation  close  to  home  that  had  high  visibility 
and  a  greater  degree  of  interagency  civilian-military 
coordination  and  planning  than  the  other  operations 
examined  to  date. 


OpcfMioci*: 

Mili* 


The  U  nintended  Consequences  of  the 
Infornnation  Age* 

(Alberts,  1996) 

T  he  purpose  of  this  analysis  is  to  identify  a  strategy  for 
introducing  and  using  Information  Age  technologies 
that  accomplishes  two  things:  first,  the  identification 
and  avoidance  of  adverse  unintended  consequences 
associated  with  the  introduction  and  utilization  of  infor- 


CAT-3 


CCRP  Publications 


mation  technologies;  and  second,  the  ability  to 
recognize  and  capitalize  on  unexpected  opportunities 


J  oint  Training  for  Infornnation  M  anagers* 

(M  axwell,  1996) 

T  his  book  proposes  new  ideas  about  joint  training  for 
information  managers  over  Command,  Control,  Com¬ 
munications  C omputers  and  I ntelligence  (C 41 )  tactical 
and  strategic  levels  It  suggests  a  substantially  new  way 
to  approach  the  training  of  future  communicators, 
grounding  its  argument  in  the  realities  of  the  fast-mov¬ 
ing  C  41  technology. 


iD  EFE-NSIVE 
INFORMATION 
WARFARE 


Defensive  Infornnation  Warfare* 

(Alberts,  1996) 

T  his  overview  of  defensive  information  warfare  is  the 
result  of  an  effort,  undertaken  at  the  request  of  the  D  ep- 
uty  Secretary  of  Defense,  to  provide  background 
material  to  participants  in  a  series  of  interagency  meet¬ 
ings  to  explore  the  nature  of  the  problem  and  to  identify 
areas  of  potential  collaboration. 


ConnnnancI,  Control,  and  the  Connnnon 

Defense 

(Allard,  1996) 

The  author  provides  an  unparalleled  basis  for  assessing 
where  we  are  and  were  we  must  go  if  we  are  to  solve  the 
joint  and  combined  command  and  control  challenges 
facing  the  U .5.  military  as  it  transitions  into  the  21st 
century. 
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Shock  &  Awe: 

Achieving  Rapid  Donninance* 

(U liman  &  Wade,  1996) 

T  he  purpose  of  this  book  is  to  explore  alternative  con¬ 
cepts  for  structuring  mission  capability  packages 
around  which  future  U .  S.  military  forces  might  be 
configured. 


Infornnation  A  ge  A  nthology: 

Volunne  I* 

(Alberts^  Papp,  1997) 

In  this  first  volume,  we  will  examine  some  of  the 
broader  issuesofthe  Information  Age:  what  the  I  nfor- 
mation  Age  is;  how  it  affects  commerce,  business,  and 
service;  what  it  means  for  the  government  and  the  mili¬ 
tary;  and  how  it  affects  international  actors  and  the 
international  system. 


IVFOKMAIION  AGL 
WIHOUKA: 


Connplexity,  Global  Politics, 
and  N  ational  Security* 

(Alberts^  Czerwinski,  1997) 

T  he  charge  given  by  the  President  of  the  N  ational 
Defense  University  and  RAND  leadership  was  three¬ 
fold:  (1)  push  the  envelope;  (2)  emphasize  the  policy  and 
strategic  dimensions  of  national  defense  with  the  impli¬ 
cations  for  C  omplexity  T  heory;  and  (3)  get  the  best 
talent  available  In  academe. 
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Target  Bosnia:  Integrating  Infornnation 
Activities  in  Peace  0  perations* 

(Siegel,  1998) 

This  book  examines  the  place  of  PI  and  PSYOP  in 
peace  operations  through  the  prism  of  NAT  0  opera- 
tionsin  Bosnia-H erzegovina. 


Coping  Mtth 
th«  Bounds 


Coping  with  the  Bounds 

(Czerwinski,  1998) 

T  he  theme  of  this  work  is  that  conventional,  or  linear, 
analysis  alone  is  not  sufficient  to  cope  with  today's  and 
tomorrow's  problems,  just  as  it  was  not  capable  of  solv¬ 
ing  yesterday's.  Its  aim  isto  convince  usto  augment  our 
efforts  with  nonlinear  insights,  and  its  hope  isto  provide 
a  basic  understanding  of  what  that  involves 


Infornnation  Warfare  and 
International  Law* 

(Greenberg,  Goodman,  &  Soo  Hoo,  1998) 

The  authors  membersof  the  Project  on  Information 
T echnology  and  I  nternational  Security  at  Stanford 
U  niversity'sC  enter  for  International  Security  and  Arms 
Control,  have  surfaced  and  explored  some  profound 
issues  that  will  shape  the  legal  context  within  which 
information  warfare  may  be  waged  and  national  infor¬ 
mation  power  exerted  in  the  coming  years 
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Lessons  Fronn  Bosnia: 

The  IFO  R  Experience* 

(Wentz,  1998) 

This  book  tellsthe  story  of  the  challenges  faced  and 
innovative  actionstaken  by  NAT 0  and  U.S.  personnel 
to  ensure  that  IFOR  and  0  peration  Joint  Endeavor 
were  military  successes  A  coherent  C  41  SR  lessons 
learned  story  has  been  pieced  together  from  firsthand 
experiences  interviews  of  key  personnel,  focused 
research,  and  analysis  of  lessons  learned  reports  pro¬ 
vided  to  the  National  Defense  U  niversity  team. 


Doing  Windows:  N  on-Traditional 
M  ilitary  Responses  to  Connplex 
Ennergencies 
(Hayes  &  Sands  1999) 

This  book  provides  the  final  results  of  a  project  spon¬ 
sored  by  the  joint  Warfare  Analysis  Center.  Our 
primary  objective  in  this  project  was  to  examine  how 
military  operations  can  support  the  long-term  objective 
of  achieving  civil  stability  and  durable  peace  in  states 
embroiled  in  complex  emergencies 


N  etwork  Centric  Warfare 

(Alberts,  Garstka,  &  Stein,  1999) 

1 1  is  hoped  that  this  book  will  contribute  to  the  prepara- 
tionsforNCW  in  two  ways  First,  by  articulating  the 
nature  of  the  characteristics  of  N  etwork  C  entric  War¬ 
fare.  Second,  by  suggesting  a  process  for  developing 
mission  capability  packages  designed  to  transform 
NCW  concepts  into  operational  capabilities. 


Doing 

Windows 


Network 

Centric 

Warfare 


CAT -7 


CCRP  Publications 


wONrRO\TA’'iC'\ 


I\I0IIM\T10\ 
C  WlPMC.V 
k'K 

rtwi 

Otermions 


CW\ 


Behind  the  Wizard's  Curtain 

(Krygiel,  1999) 

There  is  still  much  to  do  and  more  to  learn  and  under¬ 
stand  about  developing  and  fielding  an  effective  and 
durable  infostructure  as  a  foundation  for  the  21st  cen¬ 
tury.  Without  successfully  fielding  systems  of  systems 
we  will  not  be  able  to  implement  emerging  concepts  in 
adaptive  and  agile  command  and  control,  nor  will  we 
reap  the  potential  benefits  of  Network  Centric  Warfare. 


Confrontation  Analysis:  How  to  Win 
0  perations  0  ther  Than  War 
(Howard,  1999) 

A  peace  operations  campaign  (or  operation  other  than 
war)  should  be  seen  as  a  linked  sequence  of  confronta¬ 
tions  in  contrast  to  a  traditional,  warfighting  campaign, 
which  Isa  linked  sequence  of  battles  The  objective  in 
each  confrontation  is  to  bring  about  certain  "compli¬ 
ant"  behavior  on  the  part  of  other  parties  until  in  the 
end  the  campaign  objective  is  reached.  T  his  is  a  state  of 
sufficient  compliance  to  enable  the  military  to  leave  the 
theater. 


Information  Campaigns  for 
Peace  Operations 
(Avruch,  Narel,  &  Siegel,  2000) 

In  its  broadest  sense,  this  report  asks  whether  the  notion 
of  struggles  for  control  over  Information  identifiable  in 
situations  of  conflict  also  has  relevance  for  situations  of 
third-party  conflict  management— for  peace 
operations 
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Infornnation  A  ge  A  nthology: 

Volunne  II* 

(Alberts &  Papp,  2000) 

I  sthe  I  nformation  Age  bringing  with  it  new  challenges 
and  threats,  and  if  so,  what  are  they?  W  hat  sorts  of  dan¬ 
gers  will  these  challenges  and  threats  present?  From 
where  will  they  (and  do  they)  come?  I  s  information  war¬ 
fare  a  reality?  T  his  publication.  Volume  1 1  of  the 
Information  Age  Anthology,  explores  these  questions 
and  provides  preliminary  answers  to  some  of  them. 


Infornnation  A  ge  A  nthology: 

Volunne  III* 

(Alberts &  Papp,  2001) 

I  n  what  ways  will  wars  and  the  military  that  fight  them 
be  different  in  the  I  nformation  Age  than  in  earlier  ages? 
What  will  this  mean  fortheU.S.  military?  In  this  third 
volumeof  the  Information  Age  Anthology,  we  turn 
finally  to  the  task  of  exploring  answers  to  these  simply 
stated,  but  vexing  questions  that  provided  the  impetus 
for  the  first  two  volumes  of  the  I  nformation  Age 
Anthology. 


MOmUTHK  4GE 
ANTHOLOGY: 


U  nderstanding  Infornnation  A  ge  Warfare 

(Alberts,  Garstka,  Hayes,  &  Signori,  2001) 

T  his  book  presents  an  alternative  to  the  deterministic 
and  linear  strategiesofthe  planning  modernization  that 
are  now  an  artifact  of  the  I  ndustrial  Age.  T  he  approach 
being  advocated  here  begins  with  the  premise  that 
adaptation  to  the  Information  Age  centers  around  the 
ability  of  an  organization  or  an  individual  to  utilize 
information. 


Understanding 
Information  Age 
Warfare 


CCRf 
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Infornnation  A  ge  Transfornnation 

(Alberts,  2002) 

This  book  is  the  first  in  a  new  seriesof  CCRP  books 
that  will  focus  on  the  I  nformation  Age  transformation 
of  the  D epartment  of  Defense.  Accordingly,  it  deals 
with  the  issues  associated  with  a  very  large  governmen¬ 
tal  institution,  a  set  of  formidable  impediments,  both 
internal  and  external,  and  the  nature  of  the  changes 
being  brought  about  by  I  nformation  Age  concepts  and 
technologies. 


E)(PERImentat;:;\ 


Code  of  Best  Practice  for 
Experinnentation 
(CCRP,  2002) 

Experimentation  isthe  lynch  pin  in  the  DoD's strategy 
for  transformation.  W ithout  a  properly  focused,  well- 
balanced,  rigorously  designed,  and  expertly  conducted 
program  of  experimentation,  the  DoD  will  not  be  able 
to  take  full  advantage  of  the  opportunities  that  Informa¬ 
tion  Age  concepts  and  technologies  offer. 


Lessons  Fronn  Kosovo: 

The  KFO  R  Experience 

(Wentz,  2002) 

Kosovo  offered  another  uniqueopportunity  for  CCRP 
to  conduct  additional  coalition  C  41  SR  -focused  research 
in  the  areas  of  coalition  command  and  control,  civil- 
military  cooperation,  information  assurance,  C4ISR 
interoperability,  and  information  operations. 
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N  ATO  Code  of  Best  Practice  for 
C2  Assessnnent 
(2002) 

T 0  the  extent  that  they  can  be  achieved,  significantly 
reduced  levels  of  fog  and  friction  offer  an  opportunity 
for  the  military  to  develop  new  concepts  of  operations, 
new  organisational  forms,  and  new  approaches  to  com¬ 
mand  and  control,  as  well  as  to  the  processes  that 
support  it.  Analysts  will  be  increasingly  called  upon  to 
work  in  this  new  conceptual  dimension  in  order  to 
examine  the  impact  of  new  information-related  capa¬ 
bilities  coupled  with  new  waysof  organising  and 
operating. 
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Effects  Based  0  perations 

(Smith,  2003) 

T  histhird  book  of  the  I  nformation  Age  Transformation 
Series  speaks  directly  to  what  we  are  trying  to  accom¬ 
plish  on  the  "fields  of  battle"  and  argues  for  changes  in 
the  way  we  decide  what  effects  we  want  to  achieve  and 
what  means  we  will  use  to  achieve  them. 


The  Big  Issue 

(Potts,  2003) 

T  hisO  ccasional  considers  command  and  combat  in  the 
Information  Age.  It  is  an  issue  that  takes  us  into  the 
realms  of  the  unknown.  Defence  thinkers  everywhere 
are  searching  forward  for  the  science  and  alchemy  that 
will  deliver  operational  success 
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Power  to  the  Edge: 

Connnnand.. .Control...  in  the 
Infornnation  Age 
(Alberts &  Hayes,  2003) 

Powarfof/ieEdgearticulatestheprinciplesbeing  used  to 
provide  the  ubiquitous,  secure,  wideband  network  that 
people  will  trust  and  use,  populate  with  high  quality 
information,  and  use  to  develop  shared  awareness,  col¬ 
laborate  effectively,  and  synchronize  their  actions 
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